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Abstract 


The  AFIT  linear  Turbine  Cascade  Test  Facility  (ICTF)  was  used  to  study  the  effect 
of  the  small  changes  In  the  angle  of  incidence  on  turbine  blade  convective  heat  transfer. 

-5 

Other  parameters^in  the  study  were  the  model  Reynolds  number  and  the  freestream 

-r 

turbulence  level.  Characterization  test%were"^performed  to  determine  the  feasibility  of 
the  study,,followed  by  a  series  of  tests  designed  to  separate  the  effects  of  the  angle  of 
incidence,  Reynolds  number,  and  freestream  turbulence  level  on  the  convective  heat 
transfer. 

For  any  given  freestream  turbulence  level  or  angle  of  Incidence,  there  Is  an 
increase  in  the  heat  transfer  coefficient  for  an  increase  In  the  Reynolds  number.  For  the 
low  freestream  turbulence  (0.5%)  configuration  at  a  given  Reynolds  number,  there  Is  a 
decrease  in  the  convective  heat  transfer  coefficient  v/ith  an  increase  in  the  angle  of 
incidence,  partially  a  result  of  the  decrease  in  the  cascade  passage  velocity  and 
partially  a  result  of  the  variation  in  the  boundary  layer  behavior.  For  the  high  freestream 
turbulence  configuration  (10%)  at  any  Reynolds  number,  there  is  an  obvious  increase  in 
the  convective  heat  transfer  coefficient  over  that  for  the  low  turbulence  configuration; 
hov/ever,  the  transitional  or  fully  turbulent  boundary  layer  makes  it  difficult  to  observe 
any  relation  between  the  convective  heat  transfer  coefficient  and  the  angle  of  incidence. 


X 


/  Introduction 


1.1  General 

The  objective  of  this  thesis  was  to  determine  the  effect  of  the  angle  of  incidence  to 
a  linear  turbine  cascade  blade  row  on  the  turbine  blade  surface  convective  heat 
transfer.  Included  were  the  effects  due  to  the  variation  in  the  model  Reynolds  number 
and  the  freestream  turbulence  level. 

The  linear  turbine  cascade  has  been  used  as  an  experimental  device  for  predicting 
and  analyzing  the  behavior  of  operational  turbines  for  almost  1 00  years.  Because  of 
their  simplicity,  flexibility,  relative  low  cost,  and  over  100  years  of  validation,  the  linear 
turbine  cascade  allows  the  user  to  perform  a  variety  of  experiments  with  confidence  in 
the  results.  The  Air  Force  Institute  of  Technology  (AFIT)  linear  Turbine  Cascade  Test 
Facility  (TCTF)  was  the  testbed  for  this  series  of  experiments.  Described  in  detail  in 
Chapter  III,  the  TCTF  is  a  flexible  device,  allowing  variation  in  a  number  of  turbine 
cascade  parameters. 

Heat  transfer  from  a  turbine  blade  in  a  turbine  cascade  has  been  shown  to  be  a 
funotion  of  many  parameters,  including  the  freestream  turbulence  intensity,  the  model 
Reynolds  number,  the  blade  geometry,  and  the  cascade  geometry.  Many  experimental 
efforts,  including  those  of  Galassi  (Galassi,  1 989)  and  Acree  (Acree,  1 990)  using  the 
TCTF,  have  explored  the  effects  of  the  aforementioned  parameters  on  the  turbine  blade 
heat  transfer.  Each  effort  has  built  onto  the  wealth  of  k.nowledge  concerning  turbine 
blade  heat  transfer,  but  no  effort  has  supplied  all  of  the  answers.  Every  turbine  cascade 
is  different,  and  the  infiuence  of  the  individual  cascade  arrangement,  secondary 
end-wall  effects,  initial  conditions,  and  other  specifics  prevent  "an  experiment  to  end  all 
experiments";  in  addition,  advances  in  instrumentation  and  analysis  techniques  contin¬ 
ually  warrant  the  reexamination  of  many  issues. 
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1.2  Objective 

The  objective  of  this  effort  is  to  build  onto  that  wealth  of  turbine  blade  convective 
heat  transfer  knowledge  by  examining  the  turbine  blade  heat  transfer  behavior  at 
off-design  conditions.  The  principal  off-design  parameter  chosen  was  the  turbine  blade 
cascade  angle  of  incidence.  In  order  to  separate  the  effects  of  turbulence  and  model 
Reynolds  number,  a  test  plan  was  developed  to  isolate  and  identify  the  influence  of  all 
turbine  cascade  parameters. 

Unlike  an  axial  compressor  with  a  variable  inlet  guide  vane  and  variable  stators,  an 
axial  turbine  is  designed  to  operate  effectively  and  efficiently  at  one  operating  condition 
which  may  be  characterized  by  a  number  of  parameters,  including  the  model  Reynolds 
number,  freestream  turbulence,  and  angle  of  incidence;  however,  the  turbine  is  not 
always  operating  at  the  design  condition.  For  this  reason,  the  off-design  performance  is 
just  as  important  as  the  on-design  performance  and  deserves  the  same  experimental 
and  analytical  considerations. 

By  determining  the  behavior  of  the  turbine  blade  heat  transfer  at  specified  and 
controlled  design  and  off-design  conditions,  the  complicated  heat  transfer  mechanisms 
may  be  more  easily  understood,  predicted,  and  managed. 

1.3  Method 

An  opportunity  was  available  to  build  on  the  works  of  Galassi  (Galassi,  1 989)  and 
Acree  (Acree,  1 990)  using  the  TCTF  where  the  majority  of  the  hardware,  software,  and 
theory  was  in  place.  In  order  to  accommodate  a  variable  angle  of  incidence,  the  TCTF 
was  modified  upstream  of  the  cascade  test  section;  in  addition,  extensive  changes 
were  incorporated  into  the  calibration,  data  acquisition,  and  data  reduction  software. 
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After  completing  the  characterization  of  the  TCTF,  a  test  plan  was  developed 
which  altowed  variations  in  the  model  Reynolds  number,  freestream  turbulence  level, 
and  angle  of  incidence.  No  attempt  was  made  to  vary  the  amplitude  or  frequency 
nature  of  the  turbulence.  The  data  was  then  analyzed  and  compared  to  the  efforts  of 
others,  where  possible  and  relevant. 


I -3 


2.1  General 


The  theory  required  for  the  performance  of  this  thesis  was  that  of  iinear  turbine 
cascade  flow,  heat  transfer,  and  turbulence.  Since  the  objective  of  this  thesis  was 
limited  to  the  heat  transfer  behavior,  the  cascade  theory  and  turbulence  theory  will  be 
limited  to  the  topics  which  aid  the  heat  transfer  analysis. 

2.2  Cascade  Theory 

The  bulk  of  the  linear  turbine  cascade  theory  deals  with  the  performance  of  the 
cascade  as  determined  by  the  behavior  of  the  flow  through  the  cascade.  No  cascade 
"performance"  theory  is  required  for  this  thesis  other  than  that  needed  for  character¬ 
ization  of  the  cascade  flow.  A  well  behaved  cascade  will  exhibit  certain  flow  properties 
upstream,  downstream,  and  in  the  cascade  passages.  Without  these  certain  flow 
properties,  it  is  difficult  to  extrapolate  the  heat  transfer  results  to  real  turbine  behavior. 

Periodicity.  In  general,  a  well  behaved  cascade  will  exhibit  periodic  static  pressure 
and  velocity  (magnitude  and  direction)  profiles  both  upstream  and  downstream  of  the 
cascade.  The  amplitude  and  shape  of  the  periodicity  will  vary  with  model  Reynolds 
number,  the  angle  of  incidence,  and  the  turbulence  level.  It  is  beyond  the  scope  of  this 
work  to  predict  the  periodic  profiles  for  the  TCTF;  however,  periodicity  is  a  key  a 
characterization  requirement. 

Pressure  Coefficient.  The  only  other  aspect  of  cascade  theory  required  deals  with 
the  behavior  of  the  flow  in  the  cascade  as  determined  by  the  nondimensional  blade 
pressure  coefficient,  C  p ,  defined  as: 


where  p ,  is  the  local  static  pressure  on  the  blade  surface 


p„  is  the  freestream  static  pressure 
p  „  is  the  freesteam  density 
U  „  is  the  freestream  velocity 

The  term  freestream  here  refers  to  the  upstream  cascade  conditions  and  the  quantity 
is  also  known  as  the  freestream  <7  or  ,  commonly  called  the  tunnel  "head". 

The  behavior  of  the  blade  surface  C  p  distribution  is  a  good  indication  of  the 

overall  cascade  behavior.  When  plotted  as  a  function  of  position  along  the  chord,  a 
turbine  blade  pressure  distribution  is  formed.  As  the  cascade  parameters  are  varied, 
the  pressure  distribution  changes  and  can  be  observed  and  analyzed  for  characterizing 
the  cascade  flow. 

In  addition  to  providing  information  about  the  general  cascade  behavior,  manipu¬ 
lation  of  C  p  yields  the  blade  local  velocity,  an  important  parameter  in  evaluating  the 
heat  transfer.  According  to  Bernoulli’s  equation  for  incompressible  flow  along  a 
streamline, 

=  (2.2) 

where  U ,  is  the  local  velocity.  Combining  the  definition  of  C  p  and  Bernoulli’s 
equation,  the  local  velocity  U  t  can  be  determined: 

1 

2 

(2.3) 

The  local  density  p ,  is  determined  based  on  the  film  temperature  7/ 


2q„ 
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Tf  = 


(2.4) 


where  T  i  is  the  local  blade  surface  temperature 
T„  is  the  freestream  temperature 

The  values  of  U  i  and  p  ,  combine  with  the  viscosity,  p ,  and  a  geometiic  parameter 

such  as  the  surface  distance,  s ,  on  an  airfoil  to  form  the  local  Reynolds  number  Re 
defined  as: 


pf/s 
Re  =  - — 


(2.5) 


a  nondimensional  quantity  important  for  both  cascade  flow  analysis  and  heat  transfer 
analysis. 

2,3  Heat  Transfer  Theory 

Turbine  blade  heat  transfer  requires  the  understanding  of  the  conductive, 
convective,  and  radiative  heat  transfer  processes.  In  an  actual  turbine,  heat  is 
convected,  conducted,  and  radiated  from  the  hot  combustion  gases  into  ihe  turbine 
blade  surface.  In  the  tCTF,  heat  transfer  is  studied  by  heating  the  foil  surface  of  a 
specially  fabricated  turbine  blade;  in  this  case  the  heat  generated  in  the  thin  foil  surface 
is  convected  and  radiated  into  the  freestream  and  conducted  into  the  interior  of  the 
blade,  as  illustrated  in  Figure  2.1 .  Whatever  the  situation,  heat  transfer  into  or  out  of  the 
blade,  nondimensional  analysis  allows  valid  comparison. 

In  the  general  case,  the  heat  transfer  to  or  from  an  object  is  the  sum  of  the 
conductive,  convective,  and  radiative  heat  fluxes: 

(2.6) 


<T  =  Q 


rati 
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where  q  is  the  total  heat  flux  ( 1/  /  m  ^ ) 

Q  "  cond  is  the  conductive  heat  flux 
Q  "  conv  is  the  convective  heat  flux 
q"  rad  is  the  radiative  heat  flux 

Total  Heat  Flux.  For  the  TCTF,  the  heat  flux  q  is  the  heat  generated  in  the  blade 

foil  surface  as  a  result  of  current  passing  through  it  and  can  be  calculated  based  on  the 
resistance  of  the  foil: 

q toii~  ‘  fou^ foil  (2.7) 

where  /  /<,,(  is  the  current  passing  through  the  foil  surface 
R  foil  is  the  resistance  of  the  foil  per  unit  area 
The  resistivity  of  the  foil  p  is  a  function  of  the  local  foil  surface  temperature  and  has 
been  determined  (Acree,  1990)  to  be 

P/„„(n-m)  =  2.223-  10''°7(degf  )  +  8.446-  10‘7  (2.8) 

and 

=  ^  (2.8a) 

where  w  =  width  of  foil  =  G.0508  m 

t  =  thickness  of  foil  =  0.0000508  m 

Conductive  Heat  Transfer.  The  conductive  he.'^t  flux  q  "  cond  is  determined  from 
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Fourier’s  Law  of  cooling: 


(2.9) 


q' 


cond 


dx 


Where  A:  is  the  thermal  conductivity  and  dT/dx  is  the  temperature  gradient  in  the 
direction  of  the  heat  flow. 

Conductive  heat  transfer  is  a  prccej-  ..'xo-  v.-o  .lolecular  collisions  in  gases  and 
liquids  and  lattice  vibration  and  free  electron  u  nsport  in  solids;  therefore,  the  value  of 
k  \s  a  function  of  the  heat  transfer  medium  ar.<' uependent  on  the  matrjrial  compo¬ 
sition  and  temperature.  For  gases  in  particuia',  /  varies  widely  with  temperature. 

When  combined  with  the  material’s  density,  p ,  and  specific  heat,  c,  the  material’s 
temperature  dependent  thermal  diffusiviiy,  a,  is  obtained: 

a  =  —  (2,10) 

pc 

The  magnitude  of  a  is  a  good  indicator  of  the  conductive  heat  transfer  ability  of  the 
material  and  a  large  a  will  conduct  heat  faster  than  a  small  a . 

Convective  Heat  Transfer  The  convective  heat  flux  q"  is  determined  from 

Newton’s  law  of  cooling: 

Q"coa.  =  /i(7'/-7‘»)  (2  11) 

where  h  is  the  convection  heat  transfer  coefficient,  an  experimentally  determined  value. 

The  convection  coefficient  h  and  the  thermal  conductivity  k  can  be  combined  with  a 
geometrical  parameter  such  as  the  diameter  of  a  cylinder  D  to  form  the  nondimen- 
sional  Nusselt  number  Na  defined  as: 
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(2.12) 


Widely  used  in  convective  heat  transfer  analysis,  the  N  a  ha?  been  determined  both 

analytically  and  empirically  for  many  convective  heat  transfer  situations,  including  the 
turbine  blade.  Calculation  of  the  Na  will  be  explored  in  detail  later. 

Radiation  Heat  Transfer.  The  radiation  heat  flux  q"  ,ad 's  based  on  the  Stefan- 
Boltzmann  law  of  thermal  radiation: 

(2.13) 

where  e  is  the  emissivity  of  the  surface 

a  is  the  Stefan-Boltzmann  constant  =  5.669,v  lO'S-^ 

In  summary,  the  overall  heat  transfer  equation  for  the  TCTF  heat  transfer  turbine 
blade  is 

(2.14) 

The  objective  of  the  heat  transfer  experiment  Is  to  obtain  h,  the  convective  heat 
transfer  coefficient: 

^  foil  ^  foil'*'  ^  foam  ~  ^  foil^  I  ~  ”') 

h  = - - - (2.15) 

i  I  /  ee 

The  i\  n  can  then  be  calculated,  along  with  another  nondimensional  heat  transfer 
parameter,  the  St,  defined  as: 

o.-  h  _  Na 
pcU„  RePr 


(2.16) 


where  Pr  is  the  Prandtl  number  =  v/a 


V  is  the  kinematic  viscosity 

Several  empirical  means  exist  for  calculating  the  Na  and  Sf  numbers  for  a  wide 

variety  of  convective  heat  transfer  situations.  Convective  heat  transfer  is  a  boundary 
layer  mechanism,  highly  dependent  on  the  behavior  of  the  boundary  layer.  For  that 
reason,  empirical  relationships  are  often  presented  in  pairs:  one  relationship  for  the 
laminar  boundary  layer  case  and  one  relationship  for  the  turbulent  boundary  layer  case. 

Turbine  Blade  Heat  Transfer.  The  most  simplified  analysis  of  turbine  blade  heat 
transfer  approximates  the  curved  turbine  blade  surface  with  a  flat  plate.  For  the  TCTF 
configuration,  the  heated  turbine  blJ-d3  woi'!:^  be  approximated  by  a  constant  heat  flux 
flat  plate.  A  laminar  analytical  solutior.,  based  on  boundary  layer  theory,  is  given  as  a 
function  of  the  Re  and  the  Pr ,  a  common  convention: 

hx  -  - 

=  =  0.453  (2.17) 

Experimental  studies  have  verified  the  analytical  solution  above.  A  turbulent  analytical 
solution,  based  on  superposition  theory,  is  given  by, 

hx  - 

Na^  =  --  =  0,03078RefPr^  (2.18) 

K 

A  more  accurate  approximation  of  the  curved  turbine  blade  is  that  of  a  body  of 
arbitrary  shape.  These  analytical  solutions  are  more  complex  than  the  flat  plate, 
requiring  simplifying  assumptions,  yet  yield  very  good  results.  For  laminar  flow  over  an 
arbitrary  surface,  a  solution  is  given  by  Eckert  (Kays  and  Crawford,  1980): 


S(  =  0.418^; -  (2.19) 

[j  I 

0  ; 

For  the  turbulent  case,  a  procedure  by  Ambrok  (Kays  and  Crawford,  1 980)  yields  i 

j 

(J-  -J-  -^0.25  0.2 

St  =  0.0287  - — - — - - - (2.20)  1 

X  \0.2  ^  : 

|(T,-7„)'-2Sp,f/„  ^ 

_  0  J 

Regardless  of  which  approximation  is  used,  the  flat  plate  or  body  of  arbitrary  | 

- 

shape,  the  flow  across  the  leading  edge  of  the  turbine  blade  can  be  represented  as  flow  j 

across  a  cylinder,  an  important  case  of  heat  transfer  used  in  heat  exchanger  theory.  i 

For  laminar  flow  across  a  cylinder,  the  N u  takes  a  familiar  form: 

y 

I  ^ 

/Vus  =  0.81 (2.21)  j 

- 

where  Re  ^  is  calculated  using  the  cylinder  radius  R  and  the  upstream  velocity,  ^ 

2.4  Turbulence  Theory 

'Turbulent  fluid  motion  is  an  irregular  condition  of  flow  in  which  the  various 

fi 

quantities  show  a  random  variation  with  time  and  space  coordinates,  in  that  statistically  I 

distinct  average  values  can  be  discerned"  (Hinze,  1 959).  The  ability  to  calculate  distinct  | 

average  values  of  the  flow  properties  allows  for  comparative  analysis  based  on  these  ] 

J 

values.  I 

An  instantaneous  turbulent  property  X  is  represented  as  | 

1 

J 

X  =  X  +  x'  (2.22)  1 


where 


(2.23) 


f^Xdt 

T  Jo 

or,  for  discrete  data, 

X  =  IjIx,  (2.24) 

"i-i 


and  X '  represents  the  instantaneous  turbulent  fiuctuation  such  that 

F  =  0  (2.25) 

The  primary  measure  of  turbulence  used  in  this  thesis  is  the  relative  intensity  T a , 
defined  (for  the  u  velocity  component)  as 

(2.26) 

U 

where,  for  discrete  data, 

I 
Z 

(2.27) 

In  a  coordinate  system  where  K  =  1/  =  0,  the  relative  turbulence  intensity  Tu  is 

I 

Tu  =  ^ - = - ^  (2.28) 

U 

When  u’ fjns  -  rms  =  ^'rms’  ^^e  turbulence  is  said  to  be  isotropic,  where  the 
turbulent  properties  are  independent  of  coordinate  axes  rotation. 


u 
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other  measures  used  in  this  thesis  for  comparative  purposes  only  are  the  turbu¬ 
lence  microscale  and  integral  scale.  These  scales  represent  the  sizes  of  the  largest  and 
smallest  eddies  in  the  turbulent  flow.  The  methods  of  calculating  these  scales  are 
detailed  in  Appendix  A. 


Ill  Apparatus 


3. 1  TurSine  Cascade  Test  Facility 

This  investigation  was  performed  in  the  AFIT  Turbine  Cascade  Test  Facility  (TCTF) 
located  in  Building  1 9  at  Wright-Patterson  Air  Force  Base.  Designed,  built,  and  origi¬ 
nally  located  at  the  United  States  Air  Force  Academy,  the  TCTF  is  an  open  circuit, 
draw-down  tunnel  powered  by  a  20  Hp,  18-inch  diameter  centrifugal  blower  permitting 
test  section  speeds  up  to  90  m/s.  The  tunnel  speed  is  controlled  by  a  variable  vane 
arrangement  at  the  blower  inlet. 

The  TCTF  layout  is  illustrated  in  Figure  3.1 .  The  bell-mouth  inlet  directs  the  flow 
through  a  21  -inch  square,  8-inch  deep  honeycombed-type  flow  straightener.  The  flow 
continues  on  through  a  7:1  contraction  ratio  curved-wall  converging  nozzle  to  the 
4.5-inch  X  16-inch  stilling  chamber.  The  12-inch  deep  stilling  chamber,  a  new  addition 
to  the  TCTF  for  this  series  of  experiments,  has  been  outfitted  with  hand  cranks  allowing 
the  user  to  manually  adjust  the  inlet  sideboard  for  a  variable  angle  of  incidence  to  the 
linear  cascade,  illustrated  in  Figure  3.2.  The  installation  of  the  chamber  required  the 
fabrication  of  longer  Inlet  sideboards,  potentially  improving  the  cascade  inlet  flow. 

The  stilling  chamber  is  followed  by  the  24-inch  diameter  linear  turbine  cascade  test 
section.  Although  not  utilized  for  this  series  of  experiments,  the  test  section  can  be 
rotated,  also  allowing  a  variable  angle  of  incidence  to  the  linear  blade  array.  Both  the 
stilling  chamber  and  the  test  section  have  plexiglass  top  covers  allowing  visual  obser¬ 
vation  of  the  experiments.  For  turbulence  generation,  a  removeable  tube  grid  can  be 
inserted  Into  the  flow  upstream  to  the  test  section,  parallel  to  the  row  of  blades.  The 
turbulence  grid  consists  of  eight  1 /4-inch  diameter  tubes  that  span  the  height  of  the 
cascade  assembly  placed  1-inch  apart  on  center. 


After  passing  through  the  cascade,  the  flow  is  diffused  and  continues  on  through 
the  blower,  exiting  to  atmosphere.  The  open-circuit  design  results  in  test  section 
pressures  below  atmospheric  but  with  acceptable  Reynolds  numbers.  The  drawbacks 
of  the  open-circuit  design  are  high  noise  levels  and  the  ingestion  of  dust  and  insects, 
harmful  to  the  delicate  test  section  instrumentation. 

3.2  Linear  Turbine  Cascade 

The  4-blade  linear  turbine  cascade  arrangement  is  illustrated  in  Figure  3.3.  The 
reaction  stage  blade  profile  is  identical  to  that  used  in  several  studies  (Langston  et  al., 

1 977;  Moore  and  Ransmayr,  1 984)  and  is  listed  in  Appendix  B.  With  a  chord  of  4.5 
inches,  a  span  of  4.5  inches,  and  an  aspect  ratio  of  1 ,  secondary  flow  effects  limit  the 
instrumentation  area  to  the  middle  two-thirds  of  the  blade;  however,  the  large  blades 
allow  conventional  pressure  and  temperature  instrumentation  techniques.  Blades  #1 
and  #4  were  machined  from  solid  aluminuiri  at  the  USAFA.  Blades  #2  and  #3  are 
eurethane  foam,  fabricated  at  AFIT  from  molds  of  the  original  blades.  Blades  #1 ,  #3, 
and  #4  are  uninstrumented. 

Blade  #2  is  the  heart  of  the  cascade,  instrumented  for  both  surface  static  pressure 
and  temperature,  and  equipped  with  a  constant-current  heat  transfer  surface  as  illus¬ 
trated  in  Figure  3.4.  A  2-inch  wide  strip  of  very  thin  (0.002  inches)  stainless  steel  shim 
foil  was  fitted  to  the  surface  contour  of  the  blade,  beginning  1 /4-inch  from  the  trailing 
edge  of  the  suction  surface,  continuing  forward  and  around  the  leading  edge,  and 
ending  1 /4-inch  from  the  trailing  edge  on  the  pressure  side.  The  electrical  taps  are 
connected  to  the  underside  of  the  foil  and  pass  through  the  center  of  the  blade  and  exit 
through  the  bottom  endwall.  Power  for  the  heat  transfer  surface  is  supplied  by  a 
Hewlett  Packard  HP6456B  direct  current  power  supply. 


Heat  transfer  temperature  measurements  are  obtained  with  23  type  J  thermo- 
couples-soldereb  to  the  underside  of  the  foil,  and  23  static  pressure  taps  are  located 
along  a  1 /2-inch  wide  band  adjacent  to  the  foil  at  the  same  thermocouple  axial 
locations.  Both  the  thermocouple  leads  and  the  pressure  lines  pass  through  the  center 
of  the  blade  and  exit  through  the  bottom  endwall.  The  pressure  tap/thermocouple 
locations  are  listed  in  Appendix  C. 

3.3  Instrumentation 

Requirements.  The  instrumentation  for  this  series  of  experiments  consisted  of 
temperature,  pressure,  velocity,  electrical  power  measuring  devices.  The  measurement 
requirements  were  as  follows: 

temperature  -blade  surface,  23  locations 
-biade  interior,  5  locations 
-tunnel  air,  1  location 
■room  air,  1  location 

pressure  -blade  surface,  23  locations 

-tunnel  static,  1  location 
-tunnel  total,  various  locations 
-room  static,  1  location 

velocity  -tunnel,  various  locations 

power  -voltage  across  the  hf'^l  transfer  surface 

-current  through  the  heat  transfer  surface 

Temperature  Measurements.  The  blade  #2  heat  transfer  surface  was  outfitted 
with  28  type  J  thermocouples:  23  along  the  surface  and  5  embedded  in  the  blade 
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interior,  along  the  biade  centeriine  as  shown  in  Figure  3.5.  For  tunnei  temperature,  a 
singie  t^'e  J'thermocoupie  was  piaced  in  the  freestream,  and  a  singie  type  J  thermo- 
coupie  was  piace  in  the  room  alongside  the  TCTF  for  room  temperature. 

Pressure  Measurements.  For  blade  surface  pressure  measurements,  23  static 
pressure  taps  were  installed  on  blade  #2  at  the  same  axial  locations  as  the  thermo¬ 
couples.  The  tunnel  freestream  total  and  static  pressures  were  obtained  with  a  pitot- 
static  pressure  tube  located  upstream  of  the  cascade  at  mid  span  z  =  2.25  inches.  A 
pitot  tube  was  used  to  obtain  the  total  pressure  at  any  location  in  the  cascade.  Access 
to  the  cascade  test  section  was  obtained  through  the  circular  plexiglass  cover.  For  this 
experiment,  the  cover  was  modified  for  accurate  and  simple  total  pressure  scanning. 
Two  grooves  were  cut  in  the  cover  along  the  leading  and  trailing  edges  of  the  cascade, 
as  well  as  4  access  ports  to  the  cascade  passages.  Machined  aluminum  blocks  held 
the  total  pressure  tube  in  place  while  data  acquisition  took  place.  The  configuration  is 
illustrated  in  Figure  3.6. 

The  23  blade  static  pressures,  the  tunnel  total/static  pressures,  and  atmospheric 
pressure  were  measured  with  a  48  port  Scanivalve  pressure  transducer  equipped  with  a 
+/-  2  psi  pressure  module,  referenced  to  atmosphere.  The  Scanivalve  pressure  trans¬ 
ducer  received  excitation  and  signal  amplification  from  an  Endevco  model  1 09  power 
supply/model  1 06  conditioner  unit.  A  Scanivalve  CTRL2P/S2  solenoid  controller 
allowed  sequential  port  stepping  or  immediate  return  to  home  port  0  and  a  Scanivalve 
OED2  odd-even  decoder  provided  port  location  digital  readout.  The  complete  Scani¬ 
valve  port  assignments  are  listed  in  Appendix  D. 


For  tunnel  speed  setting  and  monitoring,  the  freestream  pitot-static  tube  also  was 
connected  to  a  Trimount  Instrument  U-tube  water  manometer  providing  1 .5  digits  of 
accuracy.  This  configuration  allowed  setting  the  tunnel  speed  to  within  1  m/s  of  the 
desired  speed.  The  pressure  system  schematic  is  detailed  in  Figure  3.7. 

Velocity  Measurements.  Velocity  measurements  were  obtained  with  a  2-channel 
constant  temperature  hot  wire  anemometer.  An  x-wire  of  the  arrangement  in  Figure  3.8 
was  used  to  gather  both  streamwise  and  normal  velocity  components  at  any  location  in 
the  cascade.  A  TSi  IFA-100  Intelligent  Flow  Analyzer  provided  the  controlling  circuits  for 
the  x-wire.  Depicted  Schematically  in  Figure  3.9,  the  IFA-100  combines  the  constant 
temperature  x-wire  bridge  and  controlling  circuits,  the  amplifier,  and  the  signal  condi¬ 
tioner.  The  entire  system  is  controlled  from  the  front  panel  of  the  IFA-1 00,  providing 
bridge  compensation,  gain  and  filtering,  and  fault  (broken-wire)  indication. 

Power  Measurements.  Power  to  the  heat  transfer  surface  was  estimated  based  on 
the  local  resistivity  of  the  foil  and  the  current  through  the  foil.  The  local  resistivity  of  the 
foil  Is  a  function  of  the  local  temperature  as  measured  by  the  thermocouples  and  was 
determined  by  Acree  (Acree,  1 990)  with  a  foil  resistivity  experiment.  The  current  through 
the  foil  was  measured  with  a  Weston  Electric  model  430  DC  ammeter  and  the  voltage 
across  the  foil  was  measured  with  a  Hewlett  Packard  HP3438A  digital  multimeter. 
Although  not  required  for  power  estimation,  the  voltage  across  the  foil  was  measured  to 
detect  any  current  leaks.  The  blade  power  schematic  is  illustrated  in  Figure  3.10. 

3.4  Data  Acquisition  Equipment 

The  heart  of  the  data  acquisition  system  was  the  Hewlett  Packard  HP3852A  Data 
Acquisition  and  Control  Unit  equipped  with  three  24-channel  high-speed  FET  multi¬ 
plexers  (slots  0,  2,  and  3),  an  8  channel  relay  module  (slot  1),  a  high-speed  voltmeter 
(slots  4  and  5),  and  an  integrating  voltmeter  (slot  6)  as  illustrated  in  Figure  3.1 1 .  The 


HP3852A  acquired  all  analog  transducer  signals  and  performed  the  analog  to  digital 
conversion.  The  HP3852A  was  remotely  programmed  and  controlled  by  a  Zenith  386 
PC  via  a  National  Instruments  General  Purpose  Interface  Bus  (GPIB)  as  illustrated  in 
Figure  3.12. 

All  thermocouples  were  connected  to  the  24-channel  multiplexers  in  slots  0  and  2 
and  the  signals  were  measured  by  the  integrating  voltmeter.  The  multiplexers  allowed 
direct  temperature  measurements  in  degrees  Celsius  through  temperature  compen¬ 
sation.  A  secondary  room  temperature  was  monitored  through  a  type  J  thermocouple 
and  Omega  Trendicator.  The  output  from  the  model  106  conditioner  was  connected  to 
channei  221  in  the  siot  2  multiplexer  and  the  signal  was  measured  with  the  integrating 
voltmeter.  The  Scanivalve  scanning  solenoid  was  controlled  by  the  relay  module, 
allowing  incremental  scanning  by  relay  channel  1 01  or  instant  return  to  port  0  by  relay 
channel  100.  The  high-speed  voltmeter  (100,000  measurements  per  second)  was  used 
for  the  hot-wire  anemometer  velocity  measurements.  The  output  from  the  TSi  IFA-1 00 
signal  conditioners  was  connected  to  channels  321  (wirel)  and  322  (wire2)  in  the  slot  3 
multiplexer.  The  IFA-1 00/high-speed  voltmeter  combination  allowed  hot-wire  scanning 
rates  up  to  lOOKHz.  A  complete  HP3852A  channel  assignment  listing  is  contained  in 
Appendix  E. 

3.5  Software 
3.5.1  Requirements 

In  order  to  acquire,  catalogue,  and  reduce  the  large  amounts  of  data  anticipated, 
the  TCTF  was  automated  in  the  areas  of  data  acquisition  and  data  reduction  and,  to  a 
smali  degree,  experiment  controi.  Ali  software  was  written  in  Microsoft  QuickBASIC 
version  4.5  (QB).  The  software  package,  listed  in  its  entirety  in  Appendix  F,  is  divided 
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into  3  task  oriented  subsections:  calibration,  data  acquisition,  and  data  reduction,  illus¬ 
trated  in’Figure  3.1 3.  The  tasks  (calibration,  data  acquisition,  and  data  reduction)  are 
described  in  detail  in  Chapter  IV. 

The  calibration  and  data  acquisition  software  primarily  instructs  the  HP3852A  data 
acquisition  unit  through  the  GPIB.  In  fact,  the  HP3852A  instructions  are  basic-like 
programs  within  the  QB  calibration  and  data  acquisition  programs.  For  more  infor¬ 
mation  on  programming  the  HP3852.A  through  the  PC  and  GPIB,  consult  the  HP3852A 
and  GPIB  user's  manuals.  Unlike  the  calibration  and  data  acquisition  programs,  the 
data  reduction  programs  operate  independently  of  the  HP3852A:  therefore,  all  data 
reduction  can  be  performed  with  any  QB  equipped  PC. 

3.5.2  Calibration  Software 

The  function  of  a  calibration  program  is  to  aid  the  user  in  calibrating  a 
measurement  device  and  produce  the  calibration  resuits  in  the  form  of  calibration 
constants.  Since  no  thermocouple  calibrations  were  required  (see  Section  4.2),  the 
only  calibration  programs  required  and  developed  were  PRESSCALbas  for  the  Scani- 
valve  pressure  transducer  calibration  and  XWIRECAL.bas  for  the  hot-wire  anemometer 
calibration.  Both  calibration  programs  were  initially  written  by  Galassi  (Galassi,  1 989) 
and  then  modified  by  Acree  (Acree,  1990),  but  further  refinements  were  incorporated  to 
improve  and  simplify  the  calibration  process.  Exact  details  about  the  program  steps 
can  be  found  in  the  program  comments. 

Pressure  Transducer  Calibration  Software.  The  PRESSCALbas  program  directs 
the  user  through  the  Scanivalve  pressure  transducer  calibration  process.  The  program 
originally  only  acquired  the  calibration  data  in  the  form  of  pressures  and  associated 
Scanivalve  pressure  transducer  voltages;  however,  the  program  was  modified  and  the 
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current  menu-driven  PRESSCAL.bas  program  reduces  the  calibration  data  and  writes 
the  calibration  data  and  calibration  constants  to  files  which  are  accessed  by  the  data 
acquisition  and  reduction  programs. 

Hot-Wire  Anemometer  Calibration  Software.  The  XWIRECALbas  program  directs 
the  user  through  the  hot-wire  calibration  process.  As  with  the  PRESSCALbas  program, 
the  original  XWIRECALbas  program  only  acquired  the  hot-wire  calibration  data  in  the 
form  of  velocities  and  associated  voltages  for  the  2  wires;  the  program  also  was 
modified  to  reduce  the  calibration  data  and  write  the  calibration  data  and  calibration 
constants  to  files;  additionally,  an  x-wire  alignment  subroutine  was  added  to  the  menu- 
driven  XWIRECALbas  program.  The  subroutine  allows  the  user  to  somewhat  simply 
align  the  x-wire  in  the  calibration  stand,  ensuring  a  repeatable  and  accurate  calibration. 

3.5.3  Data  Acquisition  Software 

Originally,  the  TCTF  data  acquisition  program  both  acquired  and  reduced  the 
experimental  data.  Although  very  useful  for  real-time  confirmation  of  data  validity,  the 
automated  data  reduction  aspect  of  the  program  was  time  consuming;  furthermore,  a 
well  thought  out  test  plan  is  preferred  over  a  "let’s  see  what  happens".  For  this 
experiment,  data  acquisition  was  confined  to  one  user-friendly  program  and  data 
reduction  to  other  user-friendly  programs. 

The  majority  of  the  data  acquisition  is  performed  by  the  menu-driven  data  acqui¬ 
sition  program  ACQUIRE.bas.  Functions  of  the  menu  driven  program  include  acquiring 
the  blade  pressure  and  temperature  data  and  the  cascade  flow  turbulence  data  and 
displaying  real-time  blade  temperature  and  tunnel  velocity  information. 

Temperature  Data  Acquisition  Software.  The  blade  temperature  data  is  acquired 
with  the  subroutine  TEMPACQ.  The  TEMPACQ  subroutine  instructs  the  HP3852A  to 
momentarily  interrupt  the  power  to  the  heat  transfer  surface  using  a  relay  channel,  and 


then  scan  the  30  thermocouples  with  high  speed  FET  multiplexers.  The  HP3852A 
automatically  converts  the  thermocouple  voltages  to  temperatures  in  degrees  Celsius 
and  the  30  thermocouple  data  points  are  transferred  to  the  PC  in  the  form  of  ASCII  data 
files.  Once  stored,  the  temperatures  are  then  listed  on  the  screen  for  a  validity  check. 

Pressure  Data  Acquisition  Software.  The  blade  and  tunnel  pressure  data  is 
acquired  with  the  subroutine  PRESSACQ.  In  addition  to  instructing  the  HP3852A  to 
acquire  the  Scanivalve  pressure  transducer  voltage  data  through  one  multiplexer 
channel,  the  program  also  instructs  the  HP3852A  "call  the  Scanivalve  home"  and  to  step 
the  Scanivalve  through  ports  0  -  28  through  relay  channels.  Unlike  the  thermocouple 
readings,  the  Scanivalve  pressure  transucer  voltages  are  not  converted  to  pressures; 
instead,  the  voltages  are  transferred  to  the  PC  in  an  ASCII  data  file.  Once  stored,  the 
voltages  are  then  plotted  as  a  function  of  their  pressure  tap  location,  allowing  the  user 
to  check  for  an  obvious  error.  The  pressure  data  acquisition  and  storage  process  takes 
approximately  2  minutes,  mainly  due  to  the  port  stepping  speed. 

Hot-Wire  Anemometer  Data  Acquisition  Software.  The  x-wire  data  acquisition  is 
the  most  complex  and  time  consuming  data  acquisition  process,  requiring  extensive 
revisions  to  the  original  x-wire  data  acquisition  subroutine.  Since  the  x-wire  voltages  are 
acquired  with  a  "dedicated"  high-speed  multiplexer/high-speed  voltmeter  combination, 
special  data  acquisition,  transfer,  and  storage  considerations  were  required.  The 
VELACQ  subroutine  instructs  the  HP3852A  to  acquire  a  specified  number  of  x-wire 
voltages  at  a  specified  frequency  with  the  multiplexer/voltmeter  combination.  The 
unusually  large  amount  of  sequential  (almost  simultaneous)  wire1/wire2  data  is  trans¬ 
ferred  to  the  PC  for  temporary  storage. 

After  the  data  transfer  is  complete,  the  subroutine  then  reads  in  the  data  from  the 
temporary  storage  files,  sorts  the  wirel  and  wire2  readings,  removes  the  IFA-100  signal 
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processor  gain  and  offset,  and  then  rewrites  the  voltage  data  to  the  binary  data  files, 
ready  for  data  reduction.  Binary  data  files  were  used  since  each  data  file  would  require 
over  500kb  of  storage  in  ASCII  format.  As  an  error  check,  the  first  5  and  last  5  voltages 
from  each  data  file  are  printed  to  the  screen.  Depending  upon  the  number  of  data 
points  acquired,  (16,384  per  channel  maximum),  the  x-wire  data  acquisition,  transfer, 
and  storage  process  can  take  up  to  5  minutes. 

Real-Time  Data  Display  Software.  To  aid  the  user  in  the  data  acquisition  process, 
the  real-time  display  subroutines  TSET  and  VSET  were  included  in  the  data  acquisition 
program  ACQUIRE.bas.  The  subroutine  TSET  displays  a  bar  graph  of  the  room  temper¬ 
ature  from  \iiermocouple  29  and  the  5  blade  internal  temperatures  from  thermocouples 
24  -  28.  This  display  is  useful  during  the  heat  transfer  data  acquisition  process  when 
the  tunnel  speed  and  turbulence  is  changed,  altering  the  steady-state  conditions.  The 
VSET  subroutine  is  similar  to  the  x-wire  alignment  routine  in  XWIRECALbas,  allowing 
the  user  to  properly  align  the  x-wire  in  the  flow  through  a  real-time  display  of  the  x-wire 
voltages  and  the  streamwise  and  normal  velocity  components. 

3.5.4  Data  Reduction  Software 

There  are  2  levels  of  data  reduction:  low  level,  consisting  of  the  conversion  of  raw 
data,  usually  in  the  form  of  voltages,  to  useful  data  by  means  of  the  calibration 
constants:  and  high  level,  where  the  useful  data  is  manipulated  using  the  theory  of 
Chapter  II  to  tell  the  user  something  about  what  is  going  on  in  the  experiment. 

Temperature  Data  Reduction  Software.  The  low  level  temperature  data  reduction 
is  simplified  by  the  HP3852A,  which  presents  thermocouple  data  in  the  form  of  tempera¬ 
tures  automatically.  Further  low  level  temperature  data  reduction  is  performed  by 
TREDUCE.bas  which,  using  the  thermocouple  location  file  CHORD.dat,  produces  a 
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screen  plot  of  the  temperatures  and  writes  the  thermocouple  locations  and  associated 
temperatures  to  files  for  plotting  purposes.  High  level  temperature  data  reduction  is 
performed  by  HTXFER.bas,  described  below. 

Pressure  Data  Reduction  Software.  The  low  level  and  some  high  levei  pressure 
data  reduction  is  performed  by  PREDUCE.bas.  This  program  reads  in  the  raw  pressure 
data  files  and  the  calibration  data  file  produced  by  PRESSCAL.bas  and  produces  a  file 
of  static  pressures  listing  the  static  pressure  in  kPa  for  each  of  the  28  ports; 
additionally,  the  program  calculates  the  pressure  coefficient  for  each  point  and,  using 
the  CHORD.bas  file,  produces  a  screen  plot  of  the  pressure  coefficient  profile  and  writes 
the  pressure  tap  locations  and  associated  pressure  coefficients  to  a  file  for  plotting 
purposes. 

Hot-Wire  Anemometer  Data  Reduction  Software.  Like  the  hot-wire  data  acquisition 
subroutine  VELACQ,  the  hot-wire  data  reduction  program  VELREDUC.bas  is  the  most 
complex  of  the  data  reduction  programs,  with  5  subroutines  for  performing  both  the  low 
level  and  the  high  level  data  reduction.  The  low  level  subroutines  REDUCE  and  STATS 
convert  the  hot-wire  voltages  to  velocities  and  convert  the  velocities  to  a  mean  value,  a 
standard  deviation,  and  a  turbulence  level.  For  the  high  level  data  reduction,  the 
subroutines  RealFT,  Fourt,  and  SCALES  utilize  Fast  Fourier  Transform  (FFT)  analysis 
and  the  theory  of  Chapter  II  to  calculate  the  Energy  Spectrum,  the  integral  scale,  and 
the  microcscale  of  the  turbulence.  As  with  the  other  data  reduction  programs,  the 
processed  data  is  printed  to  the  screen.  The  x-wire  data  reduction  process  is  time 
consuming  and  can  take  up  to  5  minutes. 

Heat  Transfer  Data  Reduction  Software.  Another  high  level  of  data  reduction  is 
performed  by  HTXFER.bas,  which  uses  the  temperature  data  files  and  the  pressure  data 
files  to  calculate  the  heat  transfer  properties  of  the  instrumented  blade. 
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Duplicating  the  procedures  of  TREDUCE.bas  and  PREDUCE.bas,  HTXFER.bas 
uses  the'pressure  calibration  data  and  the  thermocouple  position  data  in  files 
CHORD.dat,  S.dat,  and  DELTAL.dat,  the  temperature  and  pressure  data  to  obtain 
temperature  and  pressure  values  for  each  of  the  23  heat  transfer  surface  data  locations. 
The  data  is  further  reduced,  using  the  heat  transfer  theory  of  Chapter  II,  to  obtain  the 
Reynolds  numbers,  heat  transfer  coefficients,  Nusselt  numbers,  and  Stanton  numbers 
for  each  of  the  23  thermocouple  locations.  The  program  prints  the  raw  and  reduced 
data  to  the  screen  and  writes  the  reduced  data  to  a  data  reduction  file  which  can  be 
accessed  for  plotting  purposes. 
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IV  Procedures 


4.1  General 

This  section  describes  the  operational  procedures  followed  for  this  series  of 
experiments.  The  procedures  are  intended  to  be  a  guide  for  performing  experiments 
using  the  TCTF  with  the  aforementioned  instrumentation  and  data  acquisition 
equipment.  Recommendations  for  changes  to  these  procedures  are  listed  in  Chapter 
VI.  As  with  the  classification  of  the  software,  the  experiment  itself  is  divided  into  3 
distinct  elements  or  tasks:  calibration  of  instrumentation,  data  acquisition,  and  data 
reduction. 

4.2  Instrumentation  Calibration 

The  most  important  element  of  any  experiment  is  properly  calibrated  instrumen¬ 
tation.  Proper  calibration  will  verify  that  the  instrument  is  operating  correctly  and  will 
insure  that  the  experiment  will  yield  repeatable  results. 

Thermocouple  Calibration.  A  thermocouple  is  a  junction  of  2  dissimilar  metals 
which  generates  a  temperature-dependent  thermoelectric  potential,  as  a  result  of 
the  Seebeck  effect,  illustrated  in  Figure  4.1a.  For  small  changes  in  temperature,  E  ab 's 
linearly  proportional  to  the  junction  temperature: 

£^,  =  a7  (4.1) 

where  a  is  the  Seebeck  coefficient 
T  is  the  junction  temperature 

For  a  given  pair  of  metals  at  a  given  junction  temperature,  the  thermoelectric 
potential  will  always  be  the  same,  eliminating  a  calibration  requirement;  however, 
measuring  E^b  results  in  the  introduction  of  several  more  junctions,  illustrated  in 
Figure  4.1b,  which  now  must  be  compensated  for.  The  HP3852A  multiplexer/voltmeter 


combination  is  configured  for  this  compensation  by  measuring  the  temperature  of  the 
terminal-block,  schematically  illustrated  in  Figure  4.1c,  and  factoring  in  the  appropriate 
correction  electronically  in  real-time. 

Inexpensive  type  J  thermocouples  were  chosen  for  this  experiment.  The 
iron(-f)/constantan(-)  combination  has  a  high  Seebeck  coefficient  and  excellent 
performance  over  the  expected  temperature  range.  The  only  calibration  procedural 
requirement  for  thermocouples  is  configuring  the  HP3852A  to  measure  temperature  with 
type  J  thermocouples  (see  HP3852A  handbook). 

Pressure  Transducer  Calibration.  Unlike  the  thermocouples,  the  pressure  trans¬ 
ducer  must  be  calibrated  over  its  expected  range  of  operation.  The  Scanivalve  system 
employs  a  diaphram-type  pressure  transducer  which  uses  a  piezoresistive  element  to 
measure  the  "pressure"  applied  to  the  diaphram  to  which  it  Is  bonded.  This  diaphram/ 
piezoresistive  element  must  be  calibrated,  yielding  a  relationship  between  the  pressure 
and  the  bridge  voltage.  Since  the  bridge  voltage  is  not  referenced  to  ground,  a 
requirement  for  the  HP3852A,  the  signal  passes  through  the  Endevco  operational 
amplifier  where  it  is  amplified  and  referenced  to  ground. 

Calibration  of  the  pressure  transducer  requires  applying  a  known  pressure  to  the 
transducer  and  recording  the  associated  voltage.  A  properly  designed  pressure  trans¬ 
ducer  will  exhibit  a  linear  pressure/voltage  relationship  over  the  designed  range; 
therefore,  the  calibration  data  reduces  to  a  linear  equation: 

P  =  V„  +  mV  (4.2) 

where  Vg  is  the  transducer  voltage  as  atmospheric  pressure 
m  is  the  pressure/voltage  curve  slope 

A  detailed  pressure  transducer  calibration  procedure  is  included  in  Appendix  G. 


Hot-Wire  Anemometer  Calibration.  A  hot-wire  probe  operates  on  the  principles  of 


convective  heat  transfer.  When  placed  in  a  flowing  fluid,  the  rate  of  convective  heat 
transfer  from  a  heated  cylinder  is  proportional  to  the  velocity  of  the  fluid.  This 
relationship  is  described  by  King’s  Law, 

h  =  A  +  BV''^  (4.3) 

where  h  is  the  heat  transfer  coefficient 

A  and  B  are  constants 

V  is  the  fluid  velocity 

A  constant  temperature  hot-wire  anemometer  operates  by  maintaining  a  heated 
probe  at  a  constant  temperature.  The  probe  is  constructed  of  a  very  fine  piece  of 
material  which  has  a  high  temperature  coefficient  of  resistance.  When  used  in  the 
bridge  circuit  of  Figure  3.9,  the  temperature  maintenance  circuits  generate  a  fluctuating 
bridge  voltage  which  can  be  measured  and  correlated  to  the  heat  transfer.  This  corre¬ 
lation  can  be  calibrated,  using  known  velocities,  yielding  the  relationship 

+  (4.4) 

where  E  is  the  bridge  voltage 

A  and  B  are  calibration  constants 

V  is  the  fluid  velocity 

n  is  the  power  which  gives  the  best  linearization 

The  exponent  n  Is  a  function  of  the  fluid  and  its  value  is  determined  when  the  data  is 

linearized,  usually  using  a  least-squared  error  reduction  technique.  Several  references 
have  determined  that  for  air  n  =  0.45. 
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When  the  x-wire  of  Figure  3.8  is  placed  in  a  flowing  fluid,  each  wire  is  subjected  to 
both  a  normai  and  tangentiai  flow  component;  in  this  case,  V  in  equation  4.4  is 
repiaced  by  the  effective  veiocity  Ko// . 

=  (4-5) 

where  K  ^  is  the  veiocity  component  normai  to  the  wire 

K  7-  is  the  velocity  component  tangential  to  the  wire 
fc  is  a  yaw  correction  factor 

The  values  o\  V  ^  and  Vt  are  determined  from  the  freestream  V  and  the  angle 

between  V  and  the  wire.  The  correction  factor  A:  is  a  function  of  the  wire  geometry 
and  varies  linearly  with  the  length/diameter  ratio  of  the  probe  wire.  Details  of  the 
hot-wire  calibration  are  included  in  Appendix  G. 

4.3  Data  Acquisition 

General.  The  data  acquisition  procedure  involves  2  steps:  configuring  the  TCTF 
for  a  specific  experiment  and  acquiring  the  required  data.  In  order  to  perform  each  of 
the  2  steps,  the  user  must  develop  a  test  plan  which  states  the  desired  TCTF  configu¬ 
ration  and  the  desired  data  requirements. 

TCTF  Configuration  For  a  fixed  cascade,  the  variables  in  the  TCTF  configurations 
are  the  model  Reynolds  number,  the  freestream  turbulence  level,  and  the  cascade 
angle  of  incidence.  Changing  a  given  configuration  usually  requires  shutting  down  the 
TCTF  for  a  brief  period;  starting  the  TCTF  back  up  again  requires  a  settling  period, 
allowing  the  TCTF  to  arrive  at  steady-state  conditions,  ready  for  data  acquisition.  For 
these  reasons,  a  well  thought  out  test  matrix  prepared  in  advance  is  necessary  for 
accomplishing  the  desired  TCTF  test  in  any  reasonable  period  of  time.  The  test  matrix 
chosen  for  this  heat  transfer  experiment  is  illustrated  in  Table  4.1 . 
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Table  4.1  TCTF  Test  Configurations 


TEST 

ANGLE 

MODEL 

FREESTREAM 

NO. 

of 

REYNOLDS 

TURBULENOE 

INOIDENCE 

NUMBER 

INTENSITY 

1 

-2.5 

low 

low 

2 

-2.5 

med 

low 

3 

-2.5 

high 

low 

4 

-2.5 

low 

high 

5 

-2.5 

med 

high 

6 

-2.5 

high 

high 

7 

0.0 

low 

low 

8 

0.0 

med 

low 

1 

9 

0.0 

high 

low 

i' 

10 

0.0 

low 

high 

11 

0.0 

med 

high 

12 

0.0 

high 

high 

13 

2.5 

low 

low 

y 

14 

2.5 

med 

low 

15 

2.5 

high 

low 

16 

2.5 

low 

high 

17 

2.5 

med 

IKliHl 

18 

2.5 

high 

high 

19 

5.0 

low 

low 

20 

5.0 

med 

low 

21 

5.0 

high 

low 

> 

22 

5.0 

low 

high 

23 

5.0 

med 

high 

24 

5.0 

high 

high 

Steady  State  Conditions.  Prior  to  acquiring  any  data,  the  TCTF  steady-state 
conditions  must  be  insured.  The  steady-state  requirement  also  applies  to  the  data 
acquisition  equipment.  Of  primary  concern  is  the  period  required  for  data  acquisition 
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equipment  warm-up.  Figure  4.2  iiiustrates  the  warm-up  behavior  of  the  Scanivai- 
ve/HP3852A  combination  and  cleariy  indicates  a  minimum  2  hour  warm-up  period  for 
the  data  acquisition  equipment. 

Aithough  not  iiiustrated  here,  the  time  required  for  the  TCTF  to  arrive  at  a  steady- 
state  flow  condition  is  approximately  5  minutes;  however,  when  using  the  heat-transfer 
blade,  the  time  required  for  the  heat-transfer  surface  to  arrive  at  a  steady-state  condition 
increases  the  waiting  period  dramatically. 

Figure  4.3  illustrates  the  blade  behavior  from  power-on  at  time  0  to  a  reasonable  steady- 
state  condition  at  time  =  30  minutes.  A  change  in  the  heat-transfer  surface  due  to 
altering  the  TCTF  flow  conditions,  such  as  the  freestream  velocity  or  turbulence  level, 
exhibits  a  similar  stabilizing  behavior;  therefore,  a  reasonable  steady-state  waiting 
period  between  heat  transfer  data  acquisitions  was  determined  to  be  approximately  30 
minutes. 

The  experimental  data  is  acquired  in  three  sets:  TCTF  flow  characterizaton  data, 
heat  transfer  temperature  and  pressure  data,  and  TCTF  flow  turbulence  data.  Acquiring 
the  data  in  this  order  allows  the  user  to  identify  and  correct  any  instrumentation,  data 
acquisition,  or  flow  condition  discrepancies. 

Cascade  Flow  Characterization  Data.  Characterizing  the  flow  consisted  of 
obtaining  turbine  blade  static  pressure  scans  and  performing  velocity  magnitude/di¬ 
rection  and  total  pressure  scans  along  the  leading  and  trailing  edges  of  the  cascade. 
The  blade  static  pressure  scans  were  recorded  automatically  using  ACQUIRE.bas; 
however,  since  the  leading  and  trailing  edge  scanning  process  required  manually 
moving  the  probes  to  the  desired  locations  and  aligning  the  probes  with  the  flow,  that 
data  was  recorded  manually  on  specially  prepared  data  sheets.  For  the  velocity  magni¬ 
tude/direction  data,  the  x-wire  probe  was  aligned  with  the  streamwise  velocity  vector. 


effectively  nulling  any  normal  components;  the  program  VSET.bas  displayed  the 
real-time- velocity  magnitude  and  the  angle  was  recorded  from  the  fixture  mounted 
protractor.  For  the  total  pressure  data,  the  total  pressure  probe  was  aligned  with  the 
flow  direction  as  determined  from  the  inlet  and  exit  sideboards;  the  program 
TOTPRESS.bas  displayed  the  real  time  Scanivalve  voltages  for  conversion  to  pressures 
from  the  calibration  coefficients. 

Heat  Transfer  Data.  Unlike  the  flow  characterization  data  acquisition,  the  heat 
transfer  temperature  and  pressure  data  acquisition  is  fully  automated,  initiated  by 
executing  the  menu-driven  ACQUIRE.bas  data  acquisition  program  and  selecting  the 
appropriate  menu  choice.  The  raw  data  is  automatically  recorded,  stored,  and 
displayed,  ready  for  data  reduction. 

Turbulence  Data.  The  turbulence  data  was  collected  after  all  characterization  and 
heat  transfer  data  was  collected.  Since  the  hot-wire  probe  operates  on  the  principle  of 
correlating  a  probe  temperature/resistance  to  the  flow  velocity,  the  turbulence  data 
could  not  be  acquired  while  the  blade  was  heated.  A  more  advanced  probe  calibration 
procedure  may  be  able  to  compensate  for  the  changing  freestream  temperatures. 

The  turbulence  data  acquisition  is  the  most  time  consuming  process,  due  to  both 
the  large  number  of  data  points  acquired  for  each  measurement  and  the  number  of 
measurements.  This  time  and  data  constraint  limited  the  turbulence  data  acquisition  to 
the  low  and  high  blade  Reynolds  number  TCTF  configurations.  The  fully  automated 
data  acquisition  process  was  initiated  by  the  appropriate  menu  selection  in  the  ACQUI¬ 
RE.bas  program;  the  data  was  automatically  recorded,  stored,  and  displayed,  ready  for 
data  reduction. 


4.4  Data  Reduction 


Although  very  little  real-time  data  reduction  was  performed,  the  bulk  of  the  data 
collected  was  evaluated  at  the  end  of  the  day  it  was  acquired,  aliowing  adjustments  to 
the  test  schedule,  if  necessary.  For  the  most  part,  raw  data  reduction  was  simpie,  using 
the  software  described  in  Section  3.5. 

The  manually  recorded  velocity  and  total  pressure  characterization  data  was 
manuaily  entered  in  a  QUATTRO  PRO  spreadsheet  for  evaluation  while  the  blade  static 
pressure  characterization  data  was  processed  with  PREDUCE.bas.  The  end  products 
were  various  plots  of  the  data,  discussed  in  Chapter  V,  Resuits. 

The  raw  heat  transfer  data,  consisting  of  temperature  and  blade  static  pressure 
scans,  was  processed  with  the  low  level  TREDUCE.bas  and  PREDUCE.bas  data 
reduction  programs.  The  products  of  these  programs,  blade  surface  temperature  and 
static  pressure  plots,  were  not  intended  to  reveal  anything  about  the  heat  transfer; 
instead,  they  were  used  to  show  behaviors  and  indicate  faulty  data.  High  level  data 
reduction  was  performed  by  the  HTXFER.bas  program  which  provided  the  heat  transfer 
results  in  the  form  of  the  nondimensionai  Nusselt  number  and  Stanton  number. 

The  time  consuming  turbulence  data  reduction  was  performed  with  VELRE- 
DUC.bas.  Unlike  the  temperature  and  pressure  scans,  reducing  the  data  of  one  turbu¬ 
lence  scan  did  not  provide  enough  information  for  plotting.  Ten  turbulence  scans  were 
required  for  each  TCTF  configuration.  VELREDUC.bas  was  written  to  batch  process  the 
turbulence  data,  providing  results  that  could  be  plotted  and  further  evaluated. 


V  Results  and  Discussion 


5. 1  Tunnel  and  Cascade  Flow  Characterization 

Prior  to  conducting  any  heat  transfer  data  experiments,  it  was  necessary  to 
characterize  the  flow  through  the  cascade.  Although  both  Galassi  (Galassi,  1 989)  and 
Acree  (Acree,  1 990)  had  characterized  the  flow  for  their  expe'imental  configurations,  the 
installation  of  the  longer  variable-angle  inlet  sideboards  called  for  another  look  at  the 
flow  behavior  to  ensure  proper  cascade  flow  for  all  experimental  configurations. 

5.1.1  Cascade  Inlet  Velocity  Profiles 

General.  A  basic  indicator  of  the  cascade  behavior  is  the  cascade  inlet  velocity 
(magnitude  and  direction)  profile.  Since,  for  this  experiment,  the  TCTF  was  limited  to  a 
four-blade  cascade,  it  was  important  to  verify  the  proper  periodic  profile  for  inlet 
velocities.  Unlike  previous  efforts  with  the  TCTF,  the  velocity  profile  measurement  was 
limited  to  a  scan  across  the  cascade  inlet;  it  was  assumed  that  any  inconsistency  in  the 
velocity  profile  further  upstream  would  show  up  at  the  cascade  entrance. 

Low  and  High  Freestream  Velocity  Profiles.  The  inlet  velocity  profiles  for  both  low 
and  high  freestream  turbulence  levels  at  a  single  model  Reynolds  number  are  illustrated 
in  Figures  5.1a  and  5.1b.  Although  the  magnitude  and  direction  of  the  freestream 
velocity  were  measured,  only  the  normalized  velocity  component  normal  to  the  cascade 
inlet  is  displayed;  this  convention  combines  the  direction  and  magnitude  information. 
Normalizing  by  the  velocity  measured  at  the  leading  edge  of  the  instrumented  blade 
eliminates  any  variations  in  the  profile  caused  by  variations  in  the  freestream  velocity. 

The  velocity  profiles  indicate  that  the  flow  upstream  of  the  cascade  is  well- 
behaved.  It  appears  that  the  profile  for  the  high  turbulence  configuration  may  be 
slightly  higher  than  the  low  configuration,  probably  caused  by  the  local  high  and  low 
velocity  areas  caused  by  the  tubes  of  the  turbulence  grid.  The  odd  behavior  on  the 
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ends  of  the  turbulent  profiles  in  Figure  5.1b  is  due  to  the  changing  space  between  the 
turbulence  grid  and  the  inlet  sideboards.  As  the  angle  of  incidence  (AOI)  to  the  blade 
row  was  either  increased  or  decreased  by  adjusting  the  sideboards,  the  spacing 
between  the  sideboards  and  the  fixed,  non-adjustable  grid  either  increased  or 
decreased,  causing  a  region  of  local  acceleration  or  deceleration.  Since  this  irregularity 
did  not  seem  to  effect  the  region  of  interest,  it  was  ignored. 

5.1.2  Pressure  Coefficient  Profiles 

General.  A  blade  static  pressure  scan  was  performed  for  each  of  the  24  configura¬ 
tions  listed  in  Table  4.1  and  a  typical  low  freestream  turbulence  pressure  coefficient 
(C  p)  plot  is  illustrated  in  Figure  5.2.  Since  Cp  is  nondimensionalized  by  the  same 
upstream  quantities  p„  and  q„,the  Cp  plots  for  the  various  configurations  listed  in 
table  4.1  can  be  compared  to  each  other  for  determining  flow  differences  due  to  the 
changing  configurations. 

Most  importantly,  a  well-instrumented  blade  will  Immediately  reveal  the  location  of 
the  stagnation  point.  By  Bernoulli’s  equation  for  incompressible  flow: 

(5.1) 

along  a  stream  line.  At  the  stagnation  point  i/ ,  =  0  and 

Pi  =  P.*q«  (5.2) 

or,  substituting  for  C  p 


C 


p 


(5.3) 
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By  Figure  5.2,  the  stagnation  point  occurs  at  about  point  2.  This  point  determines  the 
dividing^ine  between  the  suction  surface  and  the  pressure  surface,  important  for  calcu¬ 
lating  the  local  heat  transfer  quantities. 

The  shape  of  the  C  p  curve  also  illustrates  the  behavior  of  the  local  velocity  of  the 

flow  around  the  blade.  Again  by  Bernoulli’s  equation 

i/,  =  /(l-Cp)  (5.4) 


For  increasing  Cp  there  is  an  adverse  pressure  gradient  and  decreasing  local  velocity 

and  for  decreasing  C  p  there  is  a  favorable  pressure  gradient  and  increasing  local 
velocity. 

On  the  pressure  surface  of  Figure  5.2,  the  overall  general  trend  is  decreasing  Cp 


corresponding  to  a  favorable  pressure  gradient  and  acceleration  of  the  flow  along  the 
surface  to  the  “throat"  of  the  cascade  passage  and  the  trailing  edge  of  the  turbine  blade, 
illustrated  in  Figure  5.3.  On  the  suction  surface  of  Figure  5.2,  the  favorable  pressure 
gradient  accelerates  the  flow  as  it  approaches  the  throat  and  then  an  adverse  pressure 
gradient  decelerates  the  flow  in  the  diffuser  region  near  the  trailing  edge  of  the  blade, 
illustrated  in  Figure  5.3. 

Acceleration  Parameter.  The  acceleration  and  deceleration  of  the  flow  along  the 
turbine  blade  surface  can  be  described  by  the  acceleration  parameter  K : 


K  = 


V  dU„ 
Li  I  ds 


(5.5) 


It  is  has  been  experimentally  shown  (Kays  et  al.,  1 970)  that  a  turbulent  boundary  layer 
tends  to  relaminarize  for  /C  >  2.5  x  1 0‘*  as  the  "viscous  sublayer  thickens  so  as  to 
engulf  the  entire  boundary  layer"  (Kays  and  Crawford,  1980);  in  this  situation,  'Ihe  local 
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heat  transfer  coefficients  are  suppressed"  (Priddy  and  Bayley,  1 985).  Figure  5.4  illus¬ 
trates  the  behavior  of  K  for  the  0  degree  AO!  TCTF  configuration.  Starting  from  zero  at 
the  stagnation  point,  the  flow  is  highly  accelerated  on  either  side  of  the  leading  edge, 
representative  of  flow  around  a  cylinder.  On  the  suction  surface,  K  drops  below  the 
critical  value  and  remains  nearly  constant;  on  the  pressure  surface,  however,  K 
exceeds  the  critical  value  from  tap  5-8,  where  the  boundary  layer  was  observed  to  be 
reattached.  Suppressed  heat  transfer  coefficients  are  expected  in  this  region. 

Flow  Separation.  Combined  with  flow  visualization  experiments,  Cp  analysis  may 

also  indicate  the  presence  of  deviations  from  the  expected  flow  behavior,  such  as  local 
regions  of  flow  separation.  On  the  pressure  surface  of  Figure  5.3,  there  is  a  small 
adverse  pressure  gradient  and  an  area  of  flow  acceleration  and  deceleration  between 
taps  2  and  4,  followed  by  a  gradual  acceleration  to  the  trailing  edge.  Extensive  oil  drop 
flow  visualization  experiments  have  verified  flow  separation  and  reattachment  in  this 
region  for  low  freestream  turbulence;  in  addition,  others  with  similarly-shaped  turbine 
blades  (Hippensteele  et  al.,  1 985;  Consigny  and  Richards,  1 982)  have  observed  compa¬ 
rable  flow  behavior  in  this  region.  On  the  suction  surface  of  Figure  5.3,  the  oil  drop 
experiments  have  indicated  flow  separation  and  reattachment  from  taps  1 8  to  1 6  for  low 
freestream  turbulence,  also  confirmed  by  the  works  of  others  (Priddy  and  Bayley,  1 985; 
Consigny  and  Richards,  1 982). 

Effect  of  Model  Reynolds  Number.  Based  on  the  turbine  blade  model  true  chord, 
the  model  Reynolds  numbers  are  2x  10®,  4x  1 0®,  and  5x10®  using  the  tunnel 
freestream  velocities  and  4x10®,  8x10®,  and  10x10®  using  the  cascade  exit 
velocities. 
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Since  C  p  is  nondimensionalized,  C  p  curves  for  the  same  cascade  configuration 

at  different  Reynolds  numbers  are  expected  to  be  the  same,  assuming  no  flow 
anomalies.  This  is  confirmed  in  Figure  5.5a,  where  the  Cp  curves  for  three  different 
Reynolds  numbers  are  nearly  identical.  Although  the  C  p  are  identical,  the  q„  are  not 
and  the  local  velocities  are  proportional  to  the  model  Reynolds  number,  illustrated  in 
Figure  5.5b. 

Effect  of  Freestream  Turbulence.  When  the  turbulence  grid  is  inserted,  increasing 
the  turbulence  level  from  0.5%  to  10%,  the  flow  along  the  suction  surface  appears  to 
gain  a  little  speed,  illustrated  in  Figures  5.6a  and  5.6b.  The  oil  drop  experiments  show 
that  the  region  of  separation  on  the  suction  surface  from  taps  18-16  has  disappeared 
and  the  flow  transitions  smoothly  from  accelerating  flow  near  the  leading  edge  to  decel¬ 
erating  flow  near  the  trailing  edge.  The  freestream  turbulence  has  "energized"  the 
boundary  layer,  making  it  adhere  to  the  blade  surface  and  delay  or  eliminate  separation. 
The  oil  drop  experiments  have  also  indicated  that  the  increased  freestream  turbulence 
has  not  eliminated  the  region  of  separation  between  taps  2  and  4  on  the  pressure 
surface. 

Effect  of  API.  As  the  API  is  increased  from  -2.5  degrees  to  5  degrees,  several 
phenomena  occur,  illustrated  in  Figuies  5.7a  and  5.7b.  Tne  stagnation  point  is 
expected  to  move  towards  the  pressure  side  as  the  API  is  increased.  Since  the 
pressure  taps  are  at  fixed  locations,  the  location  of  the  stagnation  point  could  not  be 
tracked:  however,  both  Figures  5.7a  and  5.7b  iliustrate  that  at  5  degrees  API  the 
stagnation  point  is  at  tap  2  but  moves  slightly  toward  tap  1  and  the  suction  side  as  the 
API  decreases,  as  indicated  by  the  behavior  of  taps  1  and  3.  As  the  API  increases. 
Figure  5.8  illustrates  that  the  velocity  at  tap  1  increases  and  the  velocity  of  tap  3 


decreases,  indicative  of  the  movement  of  the  stagnation  point  of  flow  around  a  cylinder. 
As  mentioned,  the  oil  drop  experiments  confirm  the  flow  separation  regions  on  both 
surfaces  at  all  AOIs  for  the  low  freestream  turbulence  configurations. 

Cascade  Passage  Streamwise  Velocity.  A  more  obvious  behavior  is  the  apparent 
"shift"  of  the  pressure  field  as  the  AOI  changes,  indicating  a  decrease  in  the  local 
velocity  on  the  pressure  and  suction  surfaces  with  an  increase  in  the  AOI,  illustrated  in 
Figure  5.8.  As  the  AOI  increases,  the  velocity  decreases  due  to  the  movement  of  the 
stagnation  point  and  the  requirement  of  the  flow  to  turn  through  a  greater  angle, 
incurring  more  losses  and  losing  energy  and  speed.  This  effect  is  further  illustrated  in 
Figures  5.9a  and  5.9b  which  show  the  behavior  of  the  midpassage  streamwise  velocity 
as  the  AOI  increases  for  both  the  low  and  high  freestream  turbulence  configurations. 
After  an  initial  velocity  rise  as  the  flow  enters  the  cascade,  the  velocity  decreases  near 
the  regions  of  maximum  curvature,  perhaps  due  to  the  suspected  regions  of  separation 
and  turbulence.  The  overall  favorable  pressure  gradient  then  accelerates  the  flow  to  a 
maximum  velocity  at  the  cascade  exit.  For  both  the  low  and  high  Reynolds  number 
configurations,  the  streamwise  velocity  at  a  given  point  increases  for  a  decrease  in  the 
AOI;  in  addition,  the  small  increase  in  the  surface  velocities  for  the  high  turbulence 
configuration  observed  in  Figure  5.6b  is  further  manifested  as  a  small  increase  in  the 
cascade  passage  streamwise  velocities  of  Figures  5.9b. 

Similar  to  the  variation  in  the  streamwise  velocity  through  the  cascade  passage, 
there  is  a  streamwise  variation  in  the  turbulence  levels  through  the  cascade  passage, 
illustrated  in  Figures  5.10a  and  5.10b.  For  the  high  turbulence  configurations,  the 
variation  in  the  turbulence  level  is  opposite  that  of  the  streamwise  velocity,  decreasing 
when  the  velocity  is  increasing  and  vice-versa.  For  the  low  turbulence  configurations, 
the  variation  in  the  turbulence  levels  are  not  as  obvious,  generally  increasing  from  0.5% 
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to  1.5%. 


Th^  vetocity  information  of  Figures  5.9a  and  5.9b  and  the  turbulence  information 
of  Figures  5.1  Oa  and  5.1  Ob  can  be  combined  to  determine  the  behavior  of  the  fluctu¬ 
ating  component  of  velocity.  By  combining  equation  2.28  and  the  assumption  of 
isotropic  turbulence  (confirmed  by  data),  the  product  of  the  mean  velocity  and 
turbulence  intensity  is 

UTu  =  u\^,  (5.6) 

where  represents  the  fluctuating  component  of  the  velocity.  Priddy  and  Bayley 

(Priddy  and  Bayley,  1 988)  have  shown  that  “the  absolute  amplitude  of  the  turbulent 
fluctuations  is  little  changed  through  the  passages,  certainly  when  compared  with  the 
observed  changes  in  the  local  turbulence  intensity."  Although  the  behavior  of  the  turbu¬ 
lence  intensities  in  this  experiment  is  not  comparable  to  the  observed  "monotonically" 
decreasing  turbulence  intensities  of  the  Priddy  and  Bayley  experiment.  Figures  5.1  Oc 
and  5.1  Od  support  the  observations  that  the  fluctuating  components  remains  nearly 
constant  through  the  cascade  passage.  In  all  cases,  low  and  high  freestream  turbu¬ 
lence,  low  and  high  Reynolds  numbers,  the  value  of  u  Vmi 's  nearly  constant  from 
location  2  through  location  4,  which  represents  the  cascade  passage.  Location  1  is  the 
freestream  and  location  5  is  beyond  the  trailing  edges  of  the  turbine  blades.  Figures 
5.10c  and  5.1  Od  further  illustrate  that  the  values  of  at  a  point  appear  nearly 
constant  over  a  range  of  AOIs. 

Although  the  C  p  plots  are  very  useful,  they  do  not  necessarily  reveal  particulars 

about  the  boundary  layer,  such  as  whether  it  is  laminar,  transitional,  or  turbulent; 
however,  the  C  p  and  heat  transfer  data  combination  provides  information  not  available 
from  either  set  of  data  alone. 
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5.2  Convective  Heat  Transfer 


Since  the  inlet  velocity  profiles  and  the  turbine  blade  C  p  plots  demonstrated 

acceptable  cascade  behavior  for  all  TCTF  configurations,  the  heat  transfer  experiment 
was  carried  out  according  to  the  test  plan  of  Table  4.1 .  A  typical  "raw"  temperature  scan 
is  illustrated  in  Figure  5.1 1  (see  Figure  3.5  for  thermocouple  locations).  Since  thermo¬ 
couples  9-12  were  located  very  close  to  the  copper  bus  bars  feeding  the  power  to  the 
foil  surface,  it  was  determined  that  those  temperature  readings  were  not  representative 
of  "normal"  surface  temperatures  as  indicated  by  the  remaining  surface  thermocouples; 
therefore,  the  temperature  readings  from  thermocouples  9  - 12  were  omitted  from  the 
heat  transfer  analysis. 

The  temperature  scan  of  Figure  5.1 1  illustrates  the  variation  in  temperatures  over 
the  surface  and  in  the  interior  of  the  blade;  however,  the  temperatures  alone  are  not 
indicative  of  the  local  convective  heat  transfer  rates  since  the  heat  transfer  rates  are  the 
result  of  a  combination  of  heat  transfer  effects.  By  including  the  conductive  and 
radiative  heat  transfer  effects  as  outlined  in  Chapter  II,  the  local  convective  heat  transfer 
coefficients  can  be  determined. 

5.2.1  Low  Turbulence  Convective  Heat  Transfer 

General.  A  typical  low  turbulence  convective  heat  transfer  coefficient  h  curve, 

generated  from  the  temperature  data  of  Figure  5.1 1 ,  is  illustrated  in  Figure  5.12.  This 
plot  illustrates  the  behavior  of  the  heat  transfer  and,  combined  with  the  C  p  plot  for  the 
configuration,  reveals  the  cascade  flow  /  heat  transfer  mechanisms. 
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Pressure  Surface.  As  expected,  the  heat  transfer  is  very  high  at  the  stagnation 
point  (tap  2)  due  to  the  very  thin  boundary  layer.  On  the  pressure  surface  (negative  s/c) 
the  heat  transfer  rate  drops  rapidly  due  to  the  observed  flow  separation  and  rises  to 
small  peak  as  the  flow  reattaches. 

In  Figure  5.12,  following  the  reattachment  peak  at  s/c  =  -0.3,  the  heat  transfer 
coefficient  appears  to  slowly  rise  as  the  flow  proceeds  to  the  trailing  edge.  The  analysis 
of  the  acceleration  parameter  K  predicted  boundary  layer  relaminarization  and  a 
suppressed  heat  transfer  coefficient  in  this  region;  however,  the  heat  transfer  coefficient 
is  Increasing,  indicating  a  transitional  or  fully  turbulent  boundary  layer.  Figure  5.1 3  illus¬ 
trates  the  experimental  behavior  of  the  local  Nusselt  number  on  the  pressure  surface  in 
comparison  to  that  predicted  for  a  flat  plate  and  an  arbitrary  surface  using  equations 
2.1 9  and  2.20.  In  the  region  from  tap  5  to  tap  8  the  experimental  Nusselt  number  values 
more  closely  match  the  turbulent  predictions  indicating  a  turbulent  boundary  layer  and 
not  a  relaminarized  boundary  layer  as  expected. 

Leading  Edge.  In  Figure  5.12,  a  peculiarity  exists  In  the  leading  edge  region  of  the 
suction  surface  where  the  heat  transfer  rate  decreases  at  tap  1  and  then  increases 
again  for  taps  23  and  22,  the  most  densely  instrumented  region  of  the  turbine  blade. 
Neither  the  Cp  plots  nor  oil-drop  flow  visualization  experiments  indicate  the  existence 
of  a  small  region  of  flow  separation  and  nothing  could  be  found  in  the  literature  on 
turbine  blade  heat  transfer  resembling  this  behavior;  however,  a  check  and  recheck  of 
the  blade  surface,  instrumentation,  and  data  acquisition  and  reduction  systems  with  the 
unheated  blade  indicated  that  the  temperature  reading  of  thermocouple  #1  may  be  2  - 
3  degrees  C  higher  than  it  should  be,  yielding  an  unusually  low  heat  transfer  coefficient. 
In  reality,  the  temperature  of  thermocouple  #1  probably  lies  in  between  the  tempera¬ 
tures  of  thermocouple  #2  and  thermocouple  #23.  In  this  case,  the  behavior  of  the  heat 
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transfer  coefficient  at  the  leading  edge  would  resemble  that  observed  in  similarly 
shaped-blades  (Consigny  and  Richards,  1982;  Daniels  and  Browne,  1981 ;  Priddy  and 
Bayley,  1985)  and  would  appear  as  in  Figure  5.14.  The  loss  of  this  data  point  limits  the 
ability  to  confirm  the  location  of  the  stagnation  point  and  to  study  in  deta;  .he  leading 
edge  heat  transfer. 

Suction  Surface.  On  the  suction  surface  of  Figure  5.1 2,  the  heat  transfer  rate 
drops  appreciably  after  tap  22,  typical  of  laminar  boundary  layer  heat  transfer  behavior. 
The  heat  transfer  rate  continues  to  steadily  drop  to  tap  17  where  it  rapidly  rises  from  tap 
1 6  to  tap  1 3.  The  rapid  rise  in  the  heat  transfer  coefficient  is  characteristic  of  a 
reattaching  boundary  layer;  a  smooth  transition  to  a  turbulent  boundary  layer  exhibits  a 
more  gradual  rise.  Figure  5.15  illustrates  the  experimental  behavior  of  the  local  Nusselt 
number  on  the  suction  surface  in  comparison  to  that  predicted  for  a  flat  plate  and  an 
arbitrary  surface  using  equations  2.19  and  2.20.  The  experimental  Nusselt  number 
values  follow  the  laminar  predictions  to  the  region  of  separation  and  reattachment, 
where  it  rises  to  a  level  between  laminar  and  fully  turbulent,  almost  meeting  the  Ambrok 
turbulent  prediction.  It  is  possible  the  boundary  layer  is  in  a  transitional  state  in  that 
region. 

Effect  of  Model  Reynolds  Number.  By  Figure  5.1 6,  the  heat  transfer  coefficients 
increase  with  increasing  Reynolds  number  and  the  suspected  regions  of  separation  and 
reattachment  remain  as  indicated  by  the  jumps  in  the  heat  transfer  coefficient  on  the 
suction  and  pressure  surfaces.  The  influence  of  the  model  Reynolds  number  on 
laminar  boundary  layer  heat  transfer  is  demonstrated  in  the  Nusselt  number  plots  of 
Figures  5.17a  and  5.17b.  As  with  the  heat  transfer  coefficient,  an  increase  in  the  local 
Reynolds  number  yields  an  increase  in  the  Nusselt  number  since,  by  equation  2.17: 

A/u  =  /(/?e°'®)  (5.7) 
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For  a  change  in  the  laminar  boundary  layer  Nusselt  number  that  is  due  solely  to  a 
changein  the  local  Reynolds  number,  the  values  of  Nu/Re°‘^  will  be  identical  for  a 
given  location.  On  the  pressure  surface  in  Figure  5.17b,  the  curves  do  not  coalesce 
further  reinforcing  the  absence  of  a  laminar  boundary  layer;  however,  also  in  Figure 
5.17b  along  the  suction  surface  from  0.2  <  s/c  <  0.7  the  curves  for  all  three  model 
Reynolds  numbers  completely  coalesce,  confirming  an  extensive  region  of  laminar 
boundary  layer  heat  transfer. 

Another  effect  of  the  variation  in  the  model  Reynolds  number  appears  in  the 
regions  of  separation  and  reattachment.  For  all  three  model  Reynolds  numbers,  the 
pressure  and  suction  surface  regions  of  flow  separation  and  reattachment  remain; 
however,  for  the  low  model  Reynolds  number  configuration  in  Figure  5.16,  the 
reattachment  on  the  suction  surface  at  s/c  =  0.7  appears  more  gradual,  indicating  more 
of  a  transition  than  a  reattachment  jump. 

Effect  of  AOI.  As  the  AOI  is  increased  from  -2.5  degrees  to  5.0  degrees,  the 
primary  observation  of  Figures  5.1 8a  is  the  decrease  in  the  heat  transfer  coefficient  with 
an  increase  in  AOI.  At  first  glance,  it  would  appear  that  variation  in  the  heat  transfer 
coefficient  is  simply  a  result  of  the  variation  in  the  local  velocity  or  local  Reynolds 
number,  observed  in  Figure  5.8  and  illustrated  in  Figure  5.18b  and  5.18c  where  the 
influences  of  the  local  Reynolds  number  are  visible.  On  both  the  pressure  and  suction 
surfaces,  the  decrease  in  the  local  Nusselt  number  is  associated  with  a  decrease  in  the 
local  Reynolds  number;  however,  by  the  same  analysis  of  Figures  5.1 7a  and  5.1 7b, 
Figures  5.19a  and  5.19b  illustrate  that  the  variation  in  the  Nusselt  number  is  not  simply 
a  result  of  the  variation  in  the  local  Reynolds  number.  Unlike  the  illustration  of  Figure 
5.17b,  the  Nusselt  number  curves  for  the  suction  surface  in  Figure  5.1 9a  do  not 
coalesce  in  Figure  5.19b,  indicating  some  effects  in  addition  to  the  variation  in  the  local 
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Reynolds  number. 

FoT  convenience,  the  heat  transfer  coefficients  for  the  -2.5  degree  AOI  and  the 
+5.0  degree  AOI  configurations  are  illustrated  in  Figure  5.20.  On  the  pressure  surfaces 
at  s/c  =  -0.1  and  on  the  suction  surfaces  at  s/c  =  0.7,  the  effect  of  the  flow  separation 
appears  to  have  a  greater  influence  for  the  5.0  degree  AOI  configuration.  The  jump  in 
the  heat  transfer  associated  with  the  reattaching  flow  decreases  with  increasing  AOI 
despite  the  nearly  equal  heat  transfer  at  the  leading  and  trailing  edges. 

5.2.2  High  Turbulence  Convective  Heat  Transfer 

General.  A  typical  high  freestream  turbulence  convective  heat  transfer  coefficient 
curve  is  illustrated  with  the  low  turbulence  curve  for  comparison  in  Figure  5.21 .  The 
primary  observation  is  the  obvious  across  the  board  increase  in  the  heat  transfer  coeffi¬ 
cient.  It  was  noted  that  increasing  the  freestream  turbulence  level  apparently  affected 
the  boundary  layer  and  the  distribution  of  the  velocity  around  the  blade.  These 
changes  are  further  manifested  in  the  distribution  of  the  heat  transfer  coefficient.  For 
the  high  turbulence  configuration,  the  heat  transfer  coefficient  appears  almost  constant 
along  the  aft  pressure  and  suction  surfaces,  very  much  unlike  the  dips  and  rises 
observed  in  the  low  turbulence  configuration. 

Pressure  and  Suction  Surfaces.  In  comparison  to  the  low  turbulence  configu¬ 
ration,  the  nearly  constant  heat  transfer  coefficients  for  the  high  turbulence  configu¬ 
ration  of  Figure  5.21  appear  to  correspond  to  turbulent  boundary  layer  throughout  both 
surfaces.  On  the  pressure  surfaces  of  Figure  5.21 ,  the  values  of  the  heat  transfer 
coefficients  at  the  trailing  edges  are  nearly  equal,  supporting  the  belief  that  the  pressure 
surface  boundary  layer  for  the  low  turbulence  configuration  is  transitional  and 
approaching  fully  turbulent  near  the  trailing  edge.  Again,  Figure  5.22a  illustrates  that 
the  experimental  Nusselt  numbers  are  best  approximated  by  the  turbulent  predictions. 
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Figure  5.22b  indicates  that  the  suction  surface  boundary  layer  is  somewhere  between 
transition  and  fully  turbulent  with  the  experimental  Nusselt  number  values  rising  steadily 
in  between  the  laminar  and  turbulent  predictions. 

Figure  5.21  illustrates  that  the  increased  freestream  turbulence  may  have  slightly 
delayed  the  flow  separation  and  reattachment  process  on  the  pressure  surface  at  -0.3 
<  s/c  <0.0.  On  suction  surfaces  of  Figure  5.21,  it  is  difficult  to  locate  any  obvious 
regions  of  separation  and  reattachment.  The  oil  drop  flow  visualization  experiments 
indicated  that  the  high  freestream  turbulence  eliminated  the  separation  observed  on  the 
suction  surface  for  low  freestream  turbulence;  however,  the  experiments  also  indicated 
turbulence  did  not  eliminate  the  separation  observed  near  the  leading  edge  of  the 
pressure  surface. 

Effect  of  Model  Reynolds  Number.  As  with  the  low  freestream  turbulence  configu¬ 
ration,  an  increase  in  the  Reynolds  number  yielded  an  increase  in  heat  transfer  coeffi¬ 
cient,  illustrated  in  Figure  5.23;  however,  unlike  the  low  freestream  turbulence 
configuration,  manipulation  of  the  Nusselt  number  using  equation  2.17  did  not  force  the 
pressure  and  suction  surface  curves  to  coalesce,  illustrated  in  Figures  5.24a  and  5.24b. 
This  confirms  the  elimination  of  the  suction  surface  laminar  boundary  with  an  increase 
in  the  freestream  turbulence. 

Effect  of  Unlike  the  low  turbulence  configuration,  there  appears  to  be  no 
obvious  connection  between  the  heat  transfer  coefficient  and  the  AOI,  illustrated  in 
Figures  5.25a,  5.25b,  and  5.25c.  Whereas  for  the  low  turbulence  configuration  there 
was  a  distinct  decrease  in  the  heat  transfer  coefficient  with  an  increase  in  the  AOI,  there 
is  no  such  behavior  for  the  high  turbulence  configuration.  It  is  apparent  that  the 
freestream  turbulence  is  affecting  the  boundary  layer  in  such  a  way  as  to  eliminate  the 
distinct  variations  observed  for  the  low  turbulence  configurations.  "Normalization"  of  the 
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local  Nusselt  numbers  are  by  /?e°  ^  in  Figures  5.26a  5.26b  provides  no  useful  infor- 
mation.-The  curves  appear  to  coalesce  near  the  leading  edge  of  the  pressure  surface, 
but  this  is  a  known  region  of  flow  separation  and  reattachment. 

The  low  and  high  turbulence  Nusselt  numbers  are  illustrated  together  in  Figures 
5.27a  and  5.27b.  On  the  pressure  surface  of  Figure  5.27a,  the  increase  in  the 
freestream  turbulence  level  appears  to  translate  into  an  across  the  board  increase  in  the 
local  Nusselt  number.  On  the  suction  surface  of  Figure  5.27b,  the  increase  in  the 
freestream  turbulence  level  appears  to  eliminate  the  region  of  laminar  heat  transfer  and 
the  boundary  layer  separation,  with  the  level  of  heat  transfer  near  the  trailing  edge  equal 
to  the  low  turbulence  values. 

In  the  final  analysis  of  the  low  and  high  freestream  turbulence  heat  transfer,  the 
correlation  given  by  Simonich  and  Bradshaw  (Simonich  and  Bradshaw,  1978), 

Sf  =  5f„(l +57u)  (5.8) 

where  Sto  is  the  low  freestream  turbulence  Stanton  number,  was  applied  to  this 

experiment  as  illustrated  in  Figures  5.28a  -  5.28h.  For  most  high  Reynolds  number 
configurations,  the  predicted  curve  agrees  with  the  experimental  curve  remarkably  well 
except  in  the  known  regions  of  boundary  layer  separation  on  the  pressure  and  suction 
surfaces:  however,  for  the  low  Reynolds  number  configurations,  there  is  little 
agreement  between  the  measured  and  predicted.  This  would  indicate  that,  despite  the 
seemingly  irregular  behavior  of  the  high  freestream  turbulence  heat  transfer,  the 
Simonich  and  Bradshaw  correlation  would  allow  a  reasonable  prediction  of  the  high 
freestream  turbulence  heat  transfer  based  cn  the  low  freestream  turbulence  values  for 
the  high  Reynolds  number  configurations. 
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VI  Conclusions  and  Recommendations 

6.1  Conclusions 

There  are  many  analytical,  experimental,  and  computational  studies  dealing  with 
the  convective  heat  transfer  to  or  from  a  turbine  blade.  Convective  heat  transfer  is  a 
process  which  is  a  function  of  properties  of  the  heat  transfer  medium.  For  forced 
convection,  the  most  important  properties  have  to  do  with  the  behavior  of  the  boundary 
layer.  In  general,  the  behavior  the  boundary  layer  is  a  function  of  freestream,  or  edge, 
properties.  In  a  turbine  cascade,  the  edge,  or  passage,  flow  properties  are  a  function  of 
the  physical  aspects  of  the  linear  cascade  which  include  the  profile  of  the  blade,  the 
span,  and  the  spacing.  Background  research  for  this  experiment  included  reviewing 
many  publications  on  previous  analytical,  computational,  and  experimental  efforts.  The 
bottom  line  is  that  every  experiment  is  different  and  the  results  of  one  experiment  may 
not  provide  general  guidance  for  all  experiments. 

The  primary  objective  of  this  effort  was  to  determine  the  effect  of  the  angle  of 
incidence  to  a  linear  turbine  blade  cascade  on  the  turbine  blade  convective  heat 
transfer.  In  order  to  separate  the  effects  of  variations  in  the  model  Reynolds  number 
and  freestream  turbulence,  a  test  matrix  was  developed  which  included  these  factors  as 
configuration  parameters.  An  attempt  was  made  to  confirm  intermediate  conclusions  by 
comparison  with  previous  efforts.  Since  this  was  not  possible  in  all  cases,  the  following 
conclusions  apply  to  the  linear  turbine  cascade  as  installed  in  the  TCTF. 

1 .  For  the  low  turbulence  configuration,  the  pressure  surface  experimental  Nusselt 
number  values  exceed  the  turbulent  boundary  layer  predictions  yet  the  convective  heat 
transfer  coefficient  indicates  a  transitional  or  fully  turbulent  boundary  layer.  Relaminari- 
zation  of  the  boundary  layer  as  predicted  by  the  acceleration  parameter  does  not  seem 
to  occur.  The  suction  surface  heat  transfer  matches  the  laminar  predictions  to  the  point 


of  separation,  after  which  the  heat  transfer  lies  between  the  laminar  and  turbulent 
predictions.  • 

For  the  high  turbulence  configuration,  the  pressure  surface  heat  transfer  appears 
fully  turbulent  and  the  entire  suction  surface  heat  transfer  lies  between  the  laminar  and 
turbulent  predictions. 

2.  For  any  TCTF  configuration,  the  convective  heat  transfer  increases  as  the 
Reynolds  number  increases.  This  is  no  surprise  since  it  can  be  shown  that,  for  both 
laminar  and  turbulent  boundary  layers,  the  heat  transfer  is  proportional  to  the  local 
Reynolds  number.  Localized  boundary  layer  behavior,  such  as  separation,  may  vary 
with  the  Reynoids  number. 

3.  For  any  TCTF  configuration,  the  convective  heat  transfer  increases  as  the 
freestream  turbulence  level  increases.  The  primary  effect  of  increased  freestream  turbu¬ 
lence  was  the  apparent  elimination  of  a  region  of  boundary  layer  separation  on  the 
suction  surface,  dramatically  increasing  the  heat  transfer.  It  was  observed  that  the  high 
freestream  turbulence  heat  transfer  may  possibly  be  predicted  with  reasonable 
accuracy  from  the  low  freestream  turbulence  using  the  Simonich-Bradshaw  correlation. 

4.  For  a  given  low  turbulence  Reynolds  number,  the  convective  heat  transfer 
decreases  as  the  AOI  increases.  This  appears  to  be  partially  a  result  of  the  decrease  in 
the  streamwise  passage  and  local  velocity  for  increasing  AOI  and  partially  a  result  of  the 
local  and  regional  behavior  of  the  boundary  layer.  It  was  demonstrated  through  laminar 
boundary  layer  heat  transfer  analysis  that  the  variation  in  heat  transfer  was  not  solely  a 
result  of  the  variation  in  the  local  Reynolds  number. 
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5.  For  a  given  high  turbulence  Reynolds  number,  the  effects  of  increasing  the  AOl 
were  less  pronounced  and  difficult  to  isolate.  Unlike  for  the  low  freestream  turbulence 
configuration,  the  variation  in  the  heat  transfer  with  the  AOl  was  not  obvious  and 
appeared  more  or  less  random. 

6.  The  results  of  the  pressure  coefficient  scans  and  the  oil  drop  flow  visualization 
experiments  became  just  as  important  as  the  results  of  the  heat  transfer  experiments  in 
determining  and  understanding  the  behavior  of  the  heat  transfer  for  changes  in  the 
angle  of  incidence,  model  Reynolds  number,  and  freestream  turbulence  level.  Any 
attempt  to  make  conclusions  based  on  the  heat  transfer  data  alone  will  result  in 
unreliable  and  unsubstantiated  conclusions.  The  heat  transfer,  pressure  coefficient, 
and  oil  drop  data  combination  provides  an  almost  complete  package  of  undisputable 
information. 

6.2  Recommendations 
6.2.1  TCTF 

1 .  Fabricate  a  new  test  section  plexiglass  cover.  After  being  used  by  3  graduate 
students,  the  present  cover  looks  like  Swiss  cheese.  The  ideal  cover  would  provide 
access  for  the  necessary  instrumentation  to  the  necessary  locations.  The  concentric 
circle  configuration  developed  by  Capt  Acree  (Acree,  1 990)  provided  limited  access  to 
almost  any  location  in  the  cascade.  Since  for  this  series  of  experiments  data  was 
required  only  for  specific  locations,  this  cover  was  removed  and  replaced  with  a  cover 
providing  rapid  and  accurate  fixed  location  access.  The  Acree  cover  could  possible  be 
reutilized  and  reconfigured  to  provide  total  access  to  any  location  in  the  cascade, 
including  upstream  and  downstream. 
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2.  Fabricate  a  new  turbulence  generation  grid.  Since  adjusting  the  AOI  changed 
the  effective  spacing  between  the  tubes,  a  new  turbulence  grid  with  adjustable  tube 
spacing  should  be  installed  in  the  TCTF.  The  new  grid  should  be  designed  such  that 
parameters  of  the  turbulence,  such  as  amplitude  and  frequency,  are  variable.  A 
"squirrel  cage"  turbulence  generation  device  has  been  successfully  used  in  a  turbine 
cascade  (Bayley  and  Priddy,  1981). 

6.2.2  Turbine  Cascade 

1 .  Fabricate  a  new  heat  transfer  instrumented  turbine  blade  with  more  pressure 
taps  and  thermocouples,  particularly  around  the  leading  edge.  A  more  densely  instrum¬ 
ented  blade  would  reveal  more  particulars  about  the  heat  transfer  processes.  With  the 
current  thermocouple  and  pressure  tap  arrangement,  it  is  sometimes  difficult  to  follow 
the  movement  of  the  stagnation  point  and  regional  flow  anomalies.  The  heat  transfer 
surface  should  be  useable  throughout,  avoiding  areas  of  unknown  current  density. 

2.  Fabricate  a  heated  (not  necessarily  instrumented)  blade  with  a  temperature 
sensitive  liquid  crystal  coating  for  heat  transfer  visualization  purposes.  Experiments 
with  the  temperature  sensitive  liquid  crystal  material  (Hippensteele  et  al,  1985)  have 
demonstrated  the  feasibility  of  the  material  for  calculating  heat  transfer  coefficients.  On 
a  simpler  scale,  in  the  uncalibrated  mode  the  material  displays  cooler  and  warmer 
regions  as  different  colors,  revealing  flow  and  heat  transfer  patterns  not  visible  with  fixed 
location  thermocouples  and  pressure  taps.  The  addition  of  the  liquid  crystal  data  to  the 
pressure,  temperature,  and  oil  flow  visualization  data  would  provide  all  the  information 
needed  for  the  heat  transfer  analysis. 

6.2.3  Instrumentation 

1 .  Replace  the  Scanivalve  with  a  pressure  data  acquisition  system  similar  to  that 
used  in  the  5  ft  wind  tunnel.  The  Scanivalve,  although  relatively  inexpensive,  is  very 
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slow  and  somewhat  unreliable.  A  new  pressure  data  acquisition  system,  combined  with 
a  more  denseiy  instrumented  turbine  blade,  would  greatly  improve  the  experimental 
value  of  the  TCTF 

2.  Install  another  high  speed  voltmeter  in  the  HP3852A,  allowing  simultaneous 
wire1  and  wire2  readings.  This  would  allow  1 00,000  Hz  x-wire  readings,  necessary  for 
accurate  calculation  of  the  turbulence  scales  through  the  Eulerian  time  correlation  (see 
Appendix  A),  and  it  would  allow  a  greater  number  of  data  points  since  each  voltmeter 
has  its  own  memory.  There  is  a  high  speed  voltmeter  in  the  5  ft  wind  tunnel  HP3852A, 
and  if  is  not  being  used  to  acquire  hot-wire  data,  it  is  being  underutilized. 

3.  Fabricate  a  total  pressure  rake  for  both  the  cascade  entrance  and  exit.  A  total 
pressure  rake  for  each  AOI  of  interest  would  speed  up  the  total  pressure  data  acqui¬ 
sition  process  and  insure  repeatable  data.  The  current  total  pressure  data  acquisition 
procedure  requires  adjusting  the  total  pressure  tube  angle  and  moving  the  tube  to  a 
location  marked  on  the  test  section  plexiglass  cover.  This  method  results  in  unreliable 
and  unrepeatable  data.  Connected  to  the  new  pressure  data  acquisition  system,  total 
pressure  information  could  be  acquired  rapidly,  repeatably,  and  reliably. 

6.2.4  Procedures 

1 .  Utilize  the  most  advanced  x-wire  calibration  scheme  available  at  AFIT.  The.'’e 
was  no  mechanism  to  factor  in  the  ambient  temperature  to  the  hot-wire  calibrations. 
During  the  course  of  a  typical  summer  day,  the  temperature  in  Bldg  1 9  may  go  from  70 
degrees  to  95  degrees.  The  current  calibration  procedure  calls  for  performing  a 
calibration  for  every  5  degrees  F.  This  convention  complicates  the  data  acquisition  and 
reduction  procedure  since  a  different  calibration  file  must  be  accessed  for  a  change  in 
temperature.  A  single  calibration  file  would  eliminate  the  possibility  of  acquiring  or 
reducing  data  with  an  incorrect  calibration  file. 
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2.  The  software  provides  almost  fully  automated  data  acquisition  and  reduction; 
howeveFi  several  more  weeks  of  work  on  the  software  would  provide  a  complete,  idiot-  ^ 

proof,  fully  automated  data  acquisition  system  for  the  TCTF.  For  this  series  of  experi-  ^ 

ments,  the  data  reduction  process  was  removed  from  the  data  acquisition  process  due 
to  the  complexity  of  the  software  at  that  time.  Due  to  the  large  amount  of  data  acquired 
for  a  typical  test,  recombination  of  the  data  acquisition  and  reduction  software  in  a  user-  J 

friendly  package  may  reduce  the  time  required  for  data  acquisition  and  reduction.  | 

6.2.5  Follow-On  or  Supporting  Efforts 

1 .  Support  the  experimental  effort  with  a  computational  effort.  I’ve  always  ; 

preferred  “building  and  testing"  over  "modeling";  however,  a  good  model  of  the  TCTF 

would  help  in  formulating  conclusions.  Many  of  the  references  compared  their  exper-  | 

imental  results  with  predictions.  | 

2.  Perform  an  experiment  with  a  variable  turbulence  level.  By  being  able  to  vary  j 

the  turbulence  level  from  0.5  %  to  1 0%,  the  somewhat  erratic  heat  transfer  behavior  j 

observed  for  the  turbulent  configurations  in  this  experiment  may  be  eliminated  to  the  | 

point  where  the  effect  of  the  AOI  may  be  observed.  I 

i 

3.  Most  importantly,  perform  a  complete  study  of  the  behavior  of  the  flow  through  | 

the  cascade  passage,  including  the  turbine  blade  boundary  layer.  In  order  to  fully  j 

understand  the  convective  heat  transfer  mechanisms  for  a  given  configuration,  the  j 

behavior  of  the  flow  must  be  fully  understood.  Although  others  have  studied  the  flow  J 

I 

through  a  turbine  cascade,  their  results  do  not  help  in  completely  understanding  what  is  | 

going  on  in  the  TCTF.  s 
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Figure  2.1 .  Turbine  Blade  Heat  Transfer  Model 
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Figure  3.1  Turbine  Cascade  Test  Facility  (TCTF)  Layout 
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Figure  3.2  Adjustable  Inlet  Sideboards 
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Figure  3.6  Test  Section  Access  Configuration 
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Figure  3.7  Pressure  Measurement  Schematic 
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Figure  3.8  X-wire  Configuration 


bridge  compensotion 


Figure  3.9  Simplified  IFA-1 00  Schematic 


Figure  3.10  Instrumented  Blade  Power  Schematic 
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Figure  3.1 1  Rear  View  of  HP3852A 
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Figure  3.12  Data  Acquisition  System  Flow  Chart 
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Figure  3.13  Software  Flow  Chari 
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Figure  4.1  Thermocouple  Schematics 
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Figure  4.2  Data  Acquisition  System  Warm-up  Behavior 


Figure  4.3  Turbine  Blade  Warm-up  Behavior 
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Figure  5.1  a  Low  Turbulence  Inlet  Velocity  Profile 


Figure  5.1b  High  Turbulence  Inlet  Velocity  Profile 
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Figure  5.4  Typical  Acceleration  Parameter 
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Figure  5.5a  Variation  in  Cp  with  Re 
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Figure  5.5b  Variation  in  Local  Velocity  with  Re 
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Figure  5.6b  Low/High  i  urbulence  Local  Velocities 


Figure  5.9a  Cascade  Passage  Streamwise  Velocity,  Low  Tu 


Figure  5.9b  Cascade  Passage  Streamwise  Velocity,  High  Tu 
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Figure  5.10a  Cascade  Passage  Tu,  Low  Re 
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Figure  5.10b  Cascade  Passage  Tu,  High  Re 
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Figure  5.10c  Cascade  Passage  RMS  Velocity  Fluctuation,  Low  Re 
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Figure  5.1 1  Typical  Turbine  Blade  Temperature  Scan 
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Figure  5.12  Typical  Convective  Heat  Transfer  Coefficient.  Low  Tu 
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Figure  5.17a  Variation  in  Local  Nu  with  Re,  Low  Tu 
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Figure  5.17b  "Normalized"  Local  Nu,  LowTu 
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Figure  5.1 8a  Variation  in  h  with  AOI,  Low  Tu 
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Figure  5.18b  Experimental  and  Analytical  Nu,  LowTu 
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Figure  5.1 8c  Experimental  and  Analytical  Nu,  Low  Tu 
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Figure  5.19a  Variation  in  Local  Nu  with  AOI,  Low  Tu 
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Figure  5.19b  "Normalized"  Local  Nu,  LowTu 
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Figure  5.20  -2.5  deg  and  5.0  deg  Low  Turbulence  h 
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Figure  5.21  Typical  Low/High  Tu  Heat  Transfer  Coefficients 


Figure  5.22a  Experimental  and  Analytical  Nu,  High  Tu 


Figure  5.22b  Experimental  and  Analytical  Nu,  High  Tu 
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Figure  5.24a  Variation  in  Local  Nu  with  Re,  High  Tu 


Figure  5.24b  "Normalized"  Local  Nu,  High  Tu 
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Figure  5.25b  Experimental  and  Analytical  Nu,  High  Tu 
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Figure  5.25c  Experimental  and  Analytical  Nu,  High  Tu 
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Figure  5.27b  Low/High  Turbulence  Local  Nu,  Suction  Surface 
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Figure  5.28a  0  deg  AOI  Stanton  Numbers,  High  Re 
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Figure  5.28e  2.5  deg  AOI  Stanton  Numbers,  High  Re 
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Figure  5.28f  2.5  deg  AOi  Stanton  Numbers,  Low  Re 
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Figure  5.28g  5.0  deg  AOI  Stanton  Numbers,  High  Re 
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Figure  5.28h  5.0  deg  AOI  Stanton  Numbers,  Low  Re 
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Appendix  A:  Turbulence  Scales  Calculations 


1.  Turbulence  Space  and  Time  Correlations 

The  statistical  nature  of  turbulent  flow  lends  itself  to  the  existence  of  correlations, 
in  both  time  and  space,  which  may  be  used  to  describe  the  turbulent  flow.  These  corre¬ 
lations  can  be  used  to  determine  the  extent  to  which  the  “fluid  motions  have  a  certain 
coiinection  with  each  other",  often  expressed  in  terms  of  the  integral  scale  and  the 
microscale. 

Turbulent  fluid  motion  is  comprised  of  smaller  and  smaller  eddies  or  cells.  The 
size  of  the  largest  eddies,  the  integral  scale,  is  a  function  of  the  flow  passage  geometry: 
the  size  of  the  smallest  eddies,  the  microscale,  is  a  function  of  the  fluid  viscosity.  By 
applying  correlation  analysis  to  the  turbulent  flow  velocity,  relationships  between  the 
flow  velocity  at  given  points  and  times  can  be  manipulated  to  yield  the  dimensions  of 
integral  scale  and  the  microscale. 

Spatial  Correlations  The  spatial  correlations,  latitudinal  and  longitudinal,  are  illus¬ 
trated  in  Figure  A1 .  Latitudinal  correlation  is  determined  by  comparing  the  fluid  velocity 
at  one  point  to  the  fluid  velocity  at  a  transversely  located  point  at  the  same  instant; 
longitudinal  correlation  is  measured  by  comparing  the  fluid  velocity  at  one  point  to  the 
fluid  velocity  at  a  point  downstream  at  the  same  instant.  The  normalized  longitudinal 
correlation  coefficient  fix')  and  lateral  correlation  coefficient  gCy)  are  defined  as 


fix) 


a'iKWi^-x) 


(A.l) 


a'2 


(A.2) 


A-1 


As  mentioned,  manipulation  of  the  correlation  yields  the  dimensions  of  the  integral 
scale  and  the  micro  scale.  It  can  be  shown  (Hinze,  1 959)  that  about  the  origin,  where 
the  correlation  is  greatest,  /  (x )  [and  g(y )]  "approaches  a  parabolic  function"  with 
the  curvature  at  the  vertex  representative  of  the  rate  of  the  local  change  in  u ' . 
Assuming  that  the  smallest  eddies  in  the  flow  are  the  cause  of  this  change,  a  microscale 
K  f  can  be  defined  such  that 

/(x)-!-'^  (A.3) 


and 


dx‘ 


Jx-O 


(A.4) 


where  -4  is  obtained  from  the  "osculating  parabola"  in  Figure  A3, 

Jx'^ 


The  integral  scale  A  f  is  defined  as 

A/  =  J^  f{x)dx  (A. 5) 

Since  the  correlation  decreases  to  zero  as  the  distance  increases,  the  above  integral 
approaches  the  value  of  the  correlation  distance,  which  is  comparable  to  the  magnitude 
of  the  largest  eddies. 

Time  Correlations  The  time  correlation  is  measured  by  corri^aring  the  fluid 
velocity  at  one  point  with  the  fluid  velocity  at  the  same  point  h'.Jt  at  different  times,  as 
illustrated  in  Figure  A2.  The  nondimensioii.al  time  correlation  coefficient  R e{x),  also 
known  as  the  Eulerian  time-correlation  coefficient,  is  defined  as 


A-2 


(A.6) 


a'(x)u'(x  -t) 


u 


Unlike  the  spatial  correlations,  manipulation  of  /?£(t)  /ields  time  scales.  Since  the 


same  logic  is  applied  to  the  definitions  of  both  the  spr-’i-i.  .nc  time  correlations,  the 
scales  are  similarly  defined.  Whereas  Xf  was  a  me.>-^e  .v‘  'hu  smallest  eddies,  the 
Eulerian  micro  time  scale 


1 

2 


(A.7) 


is  a  measure  of  smallest  correlation  times.  Similar'  >,  the  Eulerian  integral  time  scale. 


T 


f 


(A.8) 


is  a  measure  of  the  longest  correlation  times. 

The  spatial  and  time  correlations  are  related  by  Taylor’s  hypothesis 


£ 

dt 


(A.9) 


which  is  based  on  the  assumption  that  "the  fluctuations  are  too  weak  to  induce  any 
significant  motion  of  their  own,  so  that  disturbances  are  convected  along  at  the  mean 
stream  velocity"  (Cebeci, ).  In  cases  where  Taylor’s  hypothesis  holds,  it  can  be  shown 
(Hinze,  1 959)  that 

/(x)-/?,(0  (A. 10) 

2.  energy  Relationships 

The  turbulence  may  also  be  analyzed  through  the  relationship  between  the 
frequencies  and  associated  kinetic  energies,  or  strengths,  of  the  turbulent  motions.  The 
frequency/energy  relationship  has  been  defined  as 


A-3 


u'2=  E^.ia)dii  (A. 11) 

•Jo 

where  £„-(«)  is  the  "energy"  and  E n-{n)dn  is  the  nth  contribution  to  the  average 
kinetic  energy' ,  expressed  as  . 

There  exists  a  logical  relationship  between  the  correlations  and  the  energy  distri¬ 
bution.  If  the  micro  correlations  are  a  measure  of  'ie  smallest  disturbances,  both  in 
space  ar,  j  time,  ‘hey  must  be-  associated '  ;Ilh  the  high  frequency  energy  cont  ibutions; 
by  the  .same  reasoning,  the  integral  scales  are  associated  with  the  low  freque.  icy 
energy  contributions.  Through  Fourier  analysis  (Hinze,  1959),  these  relationships  are 

/Cx)  =  ==  £^.(.n)cos  -  _  dn  (A. 12) 

r 

and 

1  r“ 

Eg(t)  =  ==  dnEu.(n)cos(2nnt)dn  (A. 13) 

U'*'  VO 

With  the  energy  distribution  £^.(n)  e  common  lactor  in  the  equations  above,  the 

space  and  time  scales  can  be  expresses  in  terms  of  either  the  space  correlation  /(x), 
the  time  correlation  /?  £  ( f ) ,  or  the  energy  distribution  Eu-(n).  Manipulation  of  the 
definitions  yields 

VO*  =  4T,  (A.  14) 

=  =  (A. IS) 

and 
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dt^ 


Pjt2  /'<» 

=  =  n^EMdn  (A. 16) 
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2n2  r 
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(A. 17) 


Further  manipulation  yields 


1_  1 

2 
■/ 


2U^ 


d^Ri 

TF 


(A. 18) 


Jt-O 


and 


A/ 


(A. 19) 


Thus,  both  the  micro  and  integral  longitudinal  ••  caies  can  be  determined  from  the 
Eulerian  time  correlation  coefficient.  This  is  very  useful  since  it  is  relatively  simple  to 
obtain  and  manipulate  the  time  correlation  data. 

3.  Turbulence  Scales  using  the  Eulerian  Correlation  Coefficients 

As  mentioned,  the  time  correlation  is  measured  by  comparing  tho  fluid  velocity  at 
one  point  with  the  fluid  velocity  at  the  same  point  but  at  different  times,  such  as  with  a 
single  hot-wire  anemometer  sampling  at  a  very  high  rate.  The  velociry  O' (t)  data  is 
then  statistically  reduced  to  yield  the  mean  velocity  U  and  the  instantaneous 
fluctuation  u'(0  such  that 


U{l)  =  U  +  a'{t) 

The  Eulerian  time  correlation  coefficient  is  defined  as 


a’{x)u\x-t) 


(A. 20) 


(A.21) 


where  the  numerator  may  be  calculated  from 


A-5 


The  integral  term  is  almost  the  correlation  integral.  By  performing  the  following  change 
of  variables 


a  =  x-i-t 


x  =  t  +  a 


u'(x)  =  u'(t  +  o) 


u'(x-t)  =  u'(a) 


dx  =  da 


the  integral  term  in  the  numerator  reduces  to 


(A. 23) 
(A. 24) 
(A. 25) 
(A. 26) 
(A.27) 


(A. 28) 


which  is  the  autocorrelation  function  for  u'(t).  The  autocorrelation  function  has  the 
Fourier 

Hu<n)-HAny  (A.29) 

where  //u.(n)  =  ^e(n)  +  t/m(n)  is  obtained  from  the  Fourier  transform  of  u'(0 

such  that  the  inverse  Fourier  transform  of  H •  H will  yield  the  integral 
value. 

In  summary,  the  Eulerian  time  correlation  coefficient  7?  f  ( O  may  be  written 


/?f(0 


^ u'(a)u'(<  +  a)da 
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'2 


(A.30) 


For  discrete  data,  such  as  that  obtained  with  a  hot-wire. 


(A.31) 


where 


and 
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vY.‘^'U)u-U*0 
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n  =  iAn 
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N 
t  =  iAt 
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At  = 
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sampta 


To  obtain  the  summation  quantity  in  the  numerator  of  equation  A.31 : 

1  perform  the  discrete  Fourier  transform  of  u'(t)  to  obtain 
H  u-in)  =  RQ(n)  +  il  m{n) 


(A. 32) 
(A. 33) 
(A.34) 
(A.35) 


2  form  the  product  H^  iny  ==  Re{n)^  + 

3  perform  the  inverse  discrete  Fourier  transform  of  H H to 
obtain  the  summation  quantity 

With  the  discrete  values  of  ^^(0  determined,  the  microscale  and  integral  scale  can 

be  calculated.  A  program  for  calculating  the  scales  based  on  this  method, 
TAYLOR.bas,  is  attached. 
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4.  Turbulence  Scales  using  the  Energy  Distribution 

The  magnitude  of  the  scales  may  also  be  obtained  through  manipulations  of  the 
energy  distribution  relations  above  yielding 

1  2ii^  r” 

-^  =  =p==  n^E^.(n)cln  (A. 36) 

and 

A/=  lim  -^EAn)  (A.37) 

n-*0  4(i'2 

where  E u-(.ri)  is  derived  from  the  Fourier  transform  of  u" .  For  discrete  data,  the 
definition  of  E^-in.)  can  be  written 

iXa'2(i)=  (A.38) 

^  i-  I  n-  1 

Combining  this  with  the  discre^  form  of  Parseval’s  theorem 

=  (A.39) 

i-l  'V/i-l 

the  energy  distribution  is  obtained  from  the  Fourier  transform  of  u '  as 


where 


_HAn)l 

iV 


_1_ 

An 


(A. 40) 
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(A.41) 
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5.  Experimental  Calculations 


General.  The  concepts  presented  above  were  applied  to  turbulence  data 
acquired  in  the  TCTF  during  the  heat  transfer  experiments.  Although  the  turbulence 
scales  were  not  required  for  any  heat  transfer  analysis,  the  turbulence  data  was 
acquired  for  determining  the  turbulence  intensities,  discussed  in  Section  VI. 

Correlation  Method  Experimental  Calculations.  The  turbulence  scales  were  calcu¬ 
lated  from  equations  A.18  and  A.19  where  the  discrete  values  of  the  correlation  coeffi¬ 
cient  were  computed  using  equations  3.31  -  3.35,  incorporated  in  the  QuickBASIC 
program  TAYLOR.bas  (attached). 

Micro  Scale.  Using  the  discrete  correlation  coefficients  computed  by  the 
TAYLOR.bas  program,  the  equation  of  the  "osculating  parabola",  required  for  the  calcu¬ 
lation  of  the  micro  scale,  must  be  estimated.  Suppose  the  osculating  parabola  has  the 
form 

fCt)  =  a  +  bt*ct^  (A.42) 

Since 


and 


/(0)=1 


(A. 43) 


ilLL) 

dl 


=  b  =  0 


(A. 44) 


then 


and 


/(0=  1  +2c£^ 


(A. 45) 
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(A. 46) 


^V(0 

dt^ 


t-o 


=  2c 


Substituting  equation  A.46  into  equation  A.18 


(A.47) 


For  the  discrete  f  (  O ,  the  vaiue  of  c  corresponds  to  the  vaiue  of  the  correlation 


coefficient  at  t  =  At  such  that 


_^£(A0-I 

(AO" 


(A.48) 


Figure  A4  illustrates  the  discrete  Eulerian  time  correlation  coefficient  about 


f  =  0  for  data  acquired  at  20000  Hz.  The  curve  appears  convex  at  f  =  0  where  it 
should  be  concave;  however,  the  equation  for  the  osculating  parabola  can  be 
estimated  using  the  method  described  above  to  obtain  the  value  of  c ,  As  the  data 
acquisition  rate  is  increased,  the  correlation  coefficient  curve  more  and  more  approxi¬ 
mates  a  parabola  about  t  =  0  yielding  more  and  more  accurate  results.  When  the  data 
acquisition  rate  is  increased  to  50000  Hz,  the  maximum  attainable  for  an  X-wire  in  the 
TCTF,  a  glimpse  of  concavity  appears,  as  illustrated  in  Figure  A5.  By  going  to  a  1  -wire 
anemometer  system  and  maximizing  the  data  acquisition  rate  at  1 00000  Hz,  a  clearly 
vi/ell-behaved  correlation  coefficient  curve  appears,  as  illustrated  in  Figure  A6.  Since 
these  curves  do  not  correspond  to  identical  flow  conditions,  the  microscale  values 
cannot  be  compared  to  each  other;  however,  the  microscale  values  obtained  with  the 
correlation  method  can  be  compared  to  the  values  obtained  with  the  energy  method  of 
equations  A.36  -  A.41 . 
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Integral  Scale.  Using  the  discrete  correlation  coefficients  computed  by 
TAYLOftbas  and  the  discrete  form  of  equation  A.1 9, 

_  ^ 

(A. 48) 
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the  value  of  the  integral  scale  can  theoretically  be  determined;  however,  a  typical 
experimental  correlation  coefficient  function  is  illustrated  In  Figure  A7.  The  oscillating 
motion  of  the  function  is  probably  an  artifact  of  "noise",  either  from  the  data  acquisition 
equipment  or  the  geometry  of  the  TCTF.  Summing  over  the  entire  function  has  yielded 
Incorrect  results;  therefore,  it  was  determined  to  sum  the  coefficients  until  the  value  of 
the  coefficients  became  negative,  illustrated  in  Figure  A8.  Although  this  method  is  not 
completely  accurate,  it  approximates  a  "clean"  correlation  coefficient  function.  A  clean 
correlation  function  will  have  negative  values  but  will  not  oscillate. 

Energy  Method  Experimental  Calculations.  The  methods  developed  in  equations 
A.38  through  A.39  were  incorporated  into  the  QuickBASIC  VELREDUC.bas  program, 
listed  in  Appendix  F.  This  program  computes  the  discrete  Energy  values  for  discrete 
frequency  values.  A  typical  discrete  Energy  distribution  is  illustrated  in  Figure  A9. 

Micro  Scale.  The  microscale  calculation  was  straightforward,  based  on  the 
discrete  form  of  equation  A.36 


2n2 


I- 1 


(A. 49) 


Integral  Scale.  The  integral  scale  calculation  appears  straightforward,  based 
on  the  discrete  Energy  value  obtained  as  n  approaches  0  but  not  necessarily  equal  to 
0 .  There  are  many  methods,  which  all  yield  different  answers,  to  obtain  this  value.  The 
method  that  was  decided  upon  was  to  take  the  average  of  values  near  the  n  =  0  point: 
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(A. 50) 


f(0)  =  ^|f(0 

Correlation  Method  ys  Energy  Method.  The  following  comparisons  are  based  on 
data  acquired  for  different  tunnel  configurations;  therefore,  there  is  no  reiationship 
between  each  set  of  data,  only  between  scales  calculated  by  each  method  for  each  set 
of  data. 

From  a  typical  data  set  of  1 6384  points  acquired  at  20000  Hz  with  an  X-wire,  the 
following  numerical  experimental  results  were  obtained: 


micro  scale 
integral  scale 


Although  comparable,  the  results  of  the  two  methods  do  not  verify  the  assumption  of 
Taylor’s  hypothesis.  As  the  data  acquisition  rate  is  increased  to  50,000  Hz,  the  micro¬ 
scales,  listed  below,  remain  comparable  but  the  integral  scales  agree  less. 


micro  scale 
integral  scale 


At  the  maximum  data  acquisition  rate  of  100000  Hz  with  the  single  hot-wire,  the  micro¬ 
scales  computed  by  both  methods  are  nearly  equal;  however,  the  integral  scale  still 
shows  considerable  disagreement. 


correlation  energy 


2.509  mm 

2.186  mm 

4.025  cm 

3.254  cm 

correlation  energy 


3.772  mm 

3.082  mm 

1.481  cm 

1.010  cm 
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correlation  energy 


micro  scale 
integral  scale 


3.832  mm 

3.788  mm 

4.794  cm 

2.605  cm 

The  examples  presented  above  are  random  samples  of  the  data  acquired  for  this 
experiment;  some  data  sets  showed  better  agreement  while  others  showed  worse.  In 
general,  the  micro  scales  are  in  closer  agreement  than  the  integral  scales,  primarily  due 
to  the  difficulties  in  calculating  the  integral  scale,  described  above.  As  mentioned,  there 
are  many  methods  of  determining  Energy  value  as  n  0,  some  of  which  yield  more 
accurate  results. 

It  is  apparent  that,  as  the  data  acquisition  frequency  is  increased,  the  value  of  the 
microscale  obtained  with  the  correlation  method  approaches  that  obtained  with  Energy 
method;  in  other  words,  the  Energy  method  yields  a  more  accurate  value  regardless  of 
the  data  acquisition  frequency.  Since  it  is  desirable  to  acquire  data  at  20,000  Hz,  twice 
the  accepted  highest  fluctuation  frequency  of  1 0,000  Hz,  the  Energy  method  is 
preferred  over  the  correlation  method.  In  addition,  for  a  given  number  of  data  points,  a 
lower  data  acquisition  frequency  yields  better  resolution.  The  Energy  method  is  also 
faster  since  it  requires  performing  a  Fourier  transform  only  in  the  forward  direction;  no 
inverse  Fourier  transform  is  required  as  in  the  correlation  method.  For  these  reasons, 
the  Energy  method  was  used  to  analyze  the  turbulence  scales  in  the  TCTF. 

Turbulence  Scales  in  the  TCTF.  Figures  A1 0  and  A1 1  illustrate  typical  behaviors 
of  the  freestream  turbulence  scales  for  both  low  and  high  freestream  turbulence  at  both 
low  and  high  model  Rq.  In  general,  there  seems  to  be  no  discernable  relation  between 
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the  AOI  and  the  scale  for  both  the  microscale  and  the  integral  scale.  One  would  expect 
some  sort  of  relation  for  the  high  turbulence  configuration  since  the  effective  spacing 
between  the  tubes  in  the  fixed  turbulence  grid  decreases  as  the  AOI  increases. 

One  obvious  relation  exists  between  the  turbulence  level  and  the  scale, 
independent  of  AOI.  The  scales  are  larger  for  lower  turbulence.  This  makes  sense 
since,  for  a  fixed  velocity,  the  scales  generally  become  larger  as  the  distance  from  the 
turbulence  source  increases.  For  the  low  turbulence  configuration,  the  turbulence 
source  consists  of  the  tunnel  inlet  and  walls.  When  the  grid  is  inserted,  the  scales 
decrease  in  size. 

Another  obvious  relation  exists  between  the  freestream  velocity  and  the  scale.  For 
the  high  turbulence  configuration,  both  the  micro  and  integral  scale  increase  nearly 
linearly  with  the  freestream  velocity;  however,  for  the  low. turbulence  configuration,  the 
increase  appears  more  exponential  than  linear,  probably  due  to  the  variation  in  the 
tunnel  flow  behavior  as  the  tunnel  velocity  increases. 

The  streamwise  scale  behavior  through  a  cascade  passage  is  illustrated  in  Figures 
A1 2  through  A1 6.  The  variation  in  the  scales,  particularly  the  microscale,  mirrors  the 
variation  in  streamwise  velocity  through  the  cascade  passage,  reillustrated  in  Figure 
A17.  For  the  microscale,  the  overall  trend  is  increasing  scale  corresponding  to  the 
increasing  velocity  through  the  cascade  passage.  In  the  high  turbulence  configuration, 
the  variation  is  minimized  and  the  relationship  appears  linear  as  observed  for  the 
freestream  micro  scale.  The  integral  scale  behavior  deviates  from  the  velocity  trend  as  it 
decreases  near  the  exit  of  the  cascade  at  the  point  of  maximum  velocity:  as  with  the 
microscale,  the  variation  in  scaie  is  minimized  for  the  high  turbulence  configuration. 
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6.  Conclusion 


Since  it  was  not  the  intention  of  this  experiment  to  fuiiy  analyze  the  scale  behavior, 
no  attempt  was  made  to  relate  the  scale  behavior  to  the  turbine  blade  heat  transfer. 

The  turbulence  scale  behavior  described  above  applies  only  to  the  TCTF.  Further 
studies  involving  controlled  flow  conditions  are  recommended.  The  variation  in  velocity, 
geometry,  and  turbulence  level  did  not  allow  any  the  formation  of  absolute  conclusions 
concerning  relationships  between  the  velocity,  turbulence  level,  and  turbulence  scales. 

The  Eulerian  time  correlation  approach  to  calculating  the  turbulence  scales 
appears  to  be  in  agreement  with  the  energy  method,  thereby  "confirming"  Taylor’s 
hypothesis.  The  method  appears  to  be  more  sensitive  to  the  data  acquisition  rate  and 
"noise"  in  the  system;  in  addition,  the  Eulerian  method  is  more  computational  intensive, 
requiring  the  calculation  of  both  a  forward  and  inverse  Fourier  transform.  For  these 
reasons,  the  energy  method  is  preferred. 
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Figure  A4  Eulerian  Correlation  for  20, OCX)  ’  !z 


t  (seconds) 


Figure  A5  Eulerian  Correlation  for  50,000  Hz 


Figure  A6  Eulerian  Coneiation  for  100,000  Hz 
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Figure  A7  Typical  Euierian  Correlation  Function 
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Figure  A8  Point  Where  the  Euierian  Correlation  Crosses  the  X-Axix 
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Figure  A1 1  TCTF  Freestream  Integral  Scale  Behavior 
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Figure  A15  Low  Turbulence  Passage  Integral  Scale  Behavior 


Figure  A16  High  Turbulence  Passage  Integral  Scale  Behavior 
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TAYLOR.BAS 


This  QUICKBASIC  program  calculates  the  micro  and  integral  scale  based  on  the  corre¬ 
lation  method  described  in  Appendix  A.  Refer  to  Appendix  F  (Data  Reduction  Software) 
for  complete  subroutine  listings. 


DECLARE  SUB  STATS  (Nl,  El  IQ.  E2I0.  MEANIO.  STDDEVIO) 
DECURE  SUB  REDUCE  (EIQ,  E20,  N,  infile1$,  infile2$,  calcoefS) 
DECURE  SUB  RealFT  (datQ,  N,  ISIGN%) 

DECURE  SUB  Fourl  (datO,  NN,  ISIGN%) 

DEFINTW 

'  sea  subprogram  listings  under  VELREDUC.bas  in  Appendix  F 

'  Main  Program  Block 
’SDYNAMIC 

DIM  E1  (16384),  E2(16384),  frequ(16384) 

DIM  MEAN(2),  STDDEV(2),  MICROI(2).  INTEGRAL!(2) 

N  =  16384  'samples 

fsample  =  500001  'hertz 

dt  -  1  /  fsample  'seconds 

df  =  fsample  /  N  'hertz 

runs  =  •hxIA* 

calcoefS  =  *e:\steve\cal\044_84b.cof 
infilelS  =  *a;\turb\*  +  runS  + 'el.dat* 
infile2S  =  *a;\turb\*  +  run$  +  •e2.dat' 

RefileS  =  'e:\steva\scalas\‘  +  runS  +  'nRE.daf 
'  read  In  the  x-wire  voltages  and  converted  to  velocities 

CAU.REDUCE(E10.  E20,  N,  infilelS,  infile2S,  calcoefS) 
CALLSTATS(N,  EIQ,  E20,  MEANQ.  STDDEVQ) 

’  obtain  the  FFT  of  u'(Q  :  H(f)  stored  in  El  (i) 

PRINT :  PRINT  •  performing  FFT  on  u1  '• 
hatfn  =  N/2 

CALL  RealFT(E10,hatfn,1) 

'  calculate  H(f)  x  H(f)* 

PRINT ;  PRINT '  calculating  H(f)  X  H(0** 

E1(1)  =  E1(1)  2 

El  (2)  =  El  (2)  ^  2 
FORi  =  3TON-1  STEP2 

E1(i)  =  E1(i)  ''2  +  E1(i  +  1)  '‘2 
E1(i+  1)  =  0! 

NEXTi 

'  obtain  the  inverse  FFT  of  H(f)  x  H{f)*  to  get  the  convolution  u'(t)*u'(t) 

PRINT :  PRINT  *  performing  the  inverse  FFT  on  H(f)  X  H(0*' 

CALL  RealFT(E10,  halfn,-1) 

'  divide  the  result  of  inverse  FFT  by  n/2  factor  (algorithm  requirement) 
FOR  i  =  1  TO  N 

El  (I)  =  E1(i)/hatfn 

NEXTi 

'  calculate  the  correlation  function  by  Taylor's  Hypothesis 
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FOR i  =  1  TO  N 

=  El  (i)  /  (N  *  STDDEV(I)  2) 

'  output  the  first  25  to  a  file  for  parabola  fitting 

OPEN  RefileS  FOR  OUTPUT  AS  #1 
time  =  01 
FOR  i  =  1  TO  25 
PRINT  #1,  time.  E1(i) 
time  =  time  +  dt 
NEXT! 

CLOSE  #1 

'  calculate  the  turbulence  scales 

PRINT :  PRINT '  calculating  the  integral  scale* 

I  =  1 
sum  =  0 
WHILE  E1(0  >0 
sum  =  sum  +  El  (I) 

1  =  1  +  1 
WEND 

INTEGRALI(I)  =  MEAN(1)  *  sum  *  dt 

MICRO(I)  =  SQR((MEAN(1)  ^  2)  /  (-(El  (2)  - 1)  /  (dt  2))) 

PRINT  USING  'micro  =  #.######  integral  =  #.######*:  MICRO(I):  INTEGRAL!(1) 
END 
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Appendix  B:  Turbine  Blade  Profile  Data 


X  Y 

X  Y  (cont) 

-0.0023  0.2012 

-3.7395  -2.0635 

-0.0126  0.2516 

-3.7282  -2.1139 

-0.0352  0.3020 

-3.6867  -2.1566 

-0.0629  0.3523 

-3.6388  -2.1697 

-0.0956  0.4026 

-3.5684  -2.1566 

-0.1258  0.4530 

-3.4954  -2.0660 

-0.1598  0.5033 

-3.4400  -1.9693 

-0.2340  0.5953 

-3.3771  -1.8647 

-0.3309  0.7046 

-3.3192  -1.7615 

-0.4479  0.8053 

-3.2664  -1.6715 

-0.5914  0.9059 

-3.1997  -1.5602 

-0.7839  1.0066 

-3.1431  -1.4633 

-1 .0481  1 .0846 

-3.0827  -1.3656 

-1.2180  1.1085 

-3.0173  -1.2608 

-1.3891  1.0972 

-2.9481  -1.1601 

-1.5200  1.0659 

-2.8814  -1.0594 

-1.7011  1.0066 

-2.8059  -0.9525 

-1.8899  0.9059 

-2.7405  -0.8622 

-2.0384  0.8053 

-2.6650  -0.7626 

-2.1591  0.7046 

-2.5769  -0.6543 

-2.2598  0.6040 

-2.4913  -0.5536 

-2.3504  0.5033 

-2.4007  -0.4505 

-2.4309  0.4052 

-2.3101  -0.3598 

-2.5165  0.2968 

-2.1994  -0.2553 

-2.5870  0.2012 

■2.0748  -0.1497 

-2.6549  0.1007 

-1.9415  -0.0478 

-2.7254  -0.0025 

-1.7804  0.0503 

-2.7883  -0.0981 

-1.5653  0.1510 

-2.8487  -0.2051 

-1.4243  0.2013 

-2.9040  -0.3020 

-1.2532  0.2391 

-2.9644  -0.4026 

-1.0620  0.2504 

-3.0185  -0.5033 

-0.8845  0.2340 

-3.0701  -0.6040 

-0.7197  0.2051 

-3.1204  -0.6973 

-0.5285  0.1460 

-3.1745  -0.8053 

-0.3322  0.0453 

-3.2249  -0.9059 

-0.2655  0.0101 

-3.2727  -1.0066 

-0.1988  0.0000 

-3.3218  -1.1073 

-0.1258  0.0214 

-3.3728  -1.2101 

-0.0755  0.0541 

-3.4199  -1.3086 

-0.0289  0.1057 

-3.4627  -1.4025 

-0.0025  0.1762 

-3.5130  -1.5099 

-0.0023  0.2013 

-3.5558  -1.6079 

-3.6024  -1.7112 

-3.6439  -1.8009 

-3.6892  -1.9125 

-3.7307  -2.0132 
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Appendix  C:  Thermocouple  and  Pressure  Tap  Locations 


x/c 

s(m) 

s/c 

surface 

0.0174 

0.0100 

0.0875 

suction 

0.0652 

0.0010 

0.0087 

leading  edge 

0.1243 

0.0137 

0.1199 

pressure 

0.1913 

0.0280 

0.2450 

pressure 

0.3157 

0.0420 

0.3675 

pressure 

0.4365 

0.0615 

0.5381 

pressure 

0.6000 

0.0825 

0.7218 

pressure 

0.7687 

0.0930 

0.8136 

pressure 

0.8522 

0.1023 

0.8950 

pressure 

0.9235 

0.1085 

0.9493 

pressure 

0.9913 

0.1257 

1.0997 

trailing  edge 

0.9026  . 

0.1363 

1.1925 

suction 

0.7643 

0.1240 

1.0849 

suction 

0.6670 

0.1120 

0.9799 

suction 

0.5817 

0.1005 

0.8793 

suction 

0.4870 

0.0885 

0.7743 

suction 

0.3852 

0.0820 

0.7174 

suction 

0.3322 

0.0683 

0.5976 

suction 

0.2191 

0.0545 

0.4768 

suction 

0.1078 

0.0402 

0.3517 

suction 

0.0357 

0.0238 

0.2082 

suction 

0.0000 

0.0175 

0.1531 

suction 

0.0087 

0.0135 

0.1181 

suction 
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Appendix  D:  Scanivalve  Port  Assignments 


PORT 


ASSIGNMENT 


tap  1 


tap  2 


tap  3 


tap  4 


tap  5 


tap  6 


tap  7 


tap  8 


tap  9 


tap  10 


atmosphere 


tap  11 


tap  12 


tap  13 


tap  14 


tap  15 


tap  16 


tap  17 


tap  18 


tap  19 


tap  20 


atmosphere 


tap  21 


tap  22 


tap  23 


pitot  tube  stagnation 


atmosphere 


pitot  tube  static 


Appendix  E:  HP3852A  Channel  Assignments 


CHANNEL 

ASSIGNMENT 

000  -  006 

not  used 

007 

blade  thermocouple  1 7 

008 

blade  thermocouple  1 8 

009 

blade  thermocouple  1 9 

010 

blade  thermocouple  20 

Oil 

blade  thermocouple  21 

012 

blade  thermocouple  22 

013 

blade  thermocouple  23 

014 

blade  thermocouple  24 

015 

blade  thermocouple  25 

016 

blade  thermocouple  26 

017 

blade  thermocouple  27 

018 

blade  thermocouple  28 

019 

room  temperature  thermocouple 

020 

tunnel  temperature  thermocouple 

021  -  023 

not  used 

100 

Scanivalve  step  control 

101 

Scanivalve  home  control 

102-103 

not  used 

104 

blade  surface  power  interrupt 

105-107 

not  used 

200 

not  used 

201 

blade  thermocouple  1 

202 

blade  thermocouple  2 

203 

blade  thermocouple  3 

204 

blade  thermocouple  4 

205 

blade  thermocouple  5 

206 

blade  thermocouple  6 

207 

blade  thermocouple  7 

208 

blade  thermocouple  8 

209 

blade  thermocouple  9 

J 


CHANNEL 

ASSIGNMENT 

210 

blade  thermocouple  1 0 

211 

blade  thermocouple  1 1 

212 

blade  thermocouple  12 

213 

blade  thermocouple  1 3 

214 

blade  thermocouple  14 

215 

blade  thermocouple  1 5 

216 

blade  thermocouple  16 

217-220 

not  used 

221 

Scanivalve  pressure  transducer 

222  -  223 

not  used 

300  -  320 

not  used 

321 

X-wire  channel  1 

322 

X-wire  channel  2 

323 

not  used 

Appendix  F:  Software 


1.  Calibration  Software 


PRESSCALBAS 


■PRESSCALbas 

'  This  program  performs  the  calibration  of  the  Scanivatve  pressure 
'  transducer.  Before  executing,  prepare  for  the  calibration  by; 

'  1.  set  the  Scanivalve  on  port  0; 

'  2.  hook  up  a  pressure/vacuum  source  and  gage  to  port  0; 

'  3.  do  not  overpressurize  the  gage  (+/•  2  psi  max); 

'  4.  incrementally  apply  pressure^acuum  while  acquiring  data; 

'  5.  when  the  calibration  is  complete,  curve  fit  the  data 


DIM  CAL(30)  AS  SINGLE 

'  initialize  the  screen  and  assign  datafile  and  directory 
COLOR  7, 1 
LOCATE  3,1 

PRINT  (■  SCANIVALVE  CALIBRATION  PROGRAM*) 

VIEW  PRINT  5  TO  24 

INPLIT  'Directory  used  to  store  the  data  file?  (c:\stevo\carO  *,  ditS 
IF  dir$  =  ••  THEN  dir$  =  •c;\steve\cal\' 

INPUT  'Name  of  file  used  to  store  data?  (calib)  *,  fInameS 
IF  fInameS  =  "  THEN  fInameS  =  'calib* 

COLOR  1.7 
CLS2 

OPEN  dir$  +  fInameS  +  '.dat*  FOR  OUTPUT  AS  #1 
'  configure  the  HP38^  for  reading  the  Scanivalve  voltages 

CAaiBFIND('HP3852',  dvm%) 

CALLIBWRT(dvm%,  'RSP) 

CALL  IBWRT(dvm%,  *RST  600*) 

CAU  IBWRT(dvm%,  'REAL  C(49).  L,H,M.S') 

CALL  IBWRT(;^m%.  'USE  600;  NPLC  1*) 

‘  create  the  HP3S52  measurement  subroutine 

CALL  IBWRT(dvm%.  'SUB  VOLTAGE*) 

CALL  IBWRT  dvm%.  'AZERO  ONCE.USE  600') 

CALL  IBWRT|dvm%,  'CONF  OCV;  NROGS  50') 

CALL  IBWRT(dvm%,  'MEAS  DCV,  221 .  USE  600,  INTO  C) 

CALL  IBWRT  dvm%,  'STAT  UH.M.S.C*) 

CALL  IBWRT(dvm%.  'OISP  M*) 

CALL  IBWRT(dvm%.  'VREAD  M*) 

CAU  IBWRT(dvm%,  'SUBEND*) 

’  begin  the  main  program 

1  =  1 

PRIIVT  '<Pros8>  any  key  to  measure  the  pressure' 

DO:  LOOP  WHILE  INKEYi  =  " 

DO 

PRINT*  MEASURING  PRESSURE* 

'  measure  the  voltage,  convert  to  kPa,  and  display  to  screen 

CALL  IBWRT{dvm%.  'CALL  VOLTAGE') 
rdS  =  SPACE$(16) 

CALL  ibrd(^m%,  rd$) 
vt!  =  VALfrdS) 

PRINT  USING  *Voltage  output  is  ##.####  volts';  vt! 

INPUT  'Manometer  height  difference? dh! 
p!  =  dhl  *  .2483 

PRINT  #1,  USING  '##.####  ,  ###.#####  .  ###.##)!<#######';  dhl;  pi;  vt! 
INPUT  'Continue?  (y  or  n) ',  ans$ 
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IF  an8$  =  'n*  OR  ans$  =  'N'  THEN 
EXIT  DO 
END  IF 

PRINT :  PRINT  '<Pre88>  any  key  to  measure  the  pressure* 
DO:  LOOP  WHILE  INKEYS  =  " 

1  =  1  +  1 
LOOP 

’  reset  the  HP3852 

CALL  IBWRT(dvm%.  'RST:  CLR') 

END 
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XWIRECALBAS 


DECLARE  SUB  E1E2  (VQ.  Eft.  NPOINTS.  SWrTCH,  filenames) 

DECLARE  SUB  CONSTAI^S  (VQ,  EQ.  NPOINTS.  SWITCH,  filenames) 
•XWIRECALbae 

'Thia  calibration  program  performs  the  hot-wire  calibration  by  reading  the 
'  voltage  output  by  each  wire  and  calculating  the  corresponding  velocity 
'  based  on  the  calibration  tank  pressure  (by  manual  manometer  input)  and 
'  temperature  (by  automatic  thermocouple  input).  Before  executing,  prepare 
'  for  the  calibration  by; 


’  1.  attaching  the  tank  Tcouple  to  Tcouple  wire  #29  (Channel  019  on  the  HP3852); 
’  2.  configuring  the  IPA100  for  the  xwire  in  use; 

'  3.  connecting  the  micromanometer  to  the  calibration  tank 
DIM  V(50),  E(2.  50) 

’  initialize  the  screen 

CLS 

LOCATE  5. 7:  PRINT  •  X-wire  calibration  program  for  the  HP  3852A' 

PRINT 

DO 

CLS  :  PRINT ;  PRINT  *  MAIN  MENU’;  PRINT 
PRINT STRINGS(80.  •_■);  PRINT 
PRINT  *1)  obtain  caiitfration  voltages' 

PRINT  *2)  obtain  calibration  constants* 

PRINT '3]  EXIT 

PRINT :  PRINT  Type  your  selection  (1 , 2,  or  3);* 
chsS  =  INPUTS(I) 

'  Use  SELECT  CASE  to  process  response. 

SELECT  CASE  chsS 

CASE  ‘I*;  CALL  E1E2(V0.  EO.  NPOINTS,  SWITCH,  filenames) 

CASE  *2*;  CALL  CONSTANTS(V0.  EQ.  NPOINTS,  SWITCH,  filenames) 

CASE  *3*;  EXIT  DO 
END  SELECT 
LOOP 

END 


SUB  CONSTANTS  (VQ.  EO.  NPOINTS.  SWITCH,  filenames) 

’  this  subroutine  calculates  the  X-wire  calibration  constants 

CLS 

DIM  X(2, 50),  Y(2. 50),  a(2),  B0(2),  B1  (2) 

DIM  ueff(2,  50) 

DIM  Eout(2,  50) 

'  if  calibration  data  was  just  collected,  go  to  line  20  and  continue 
IF  SWITCH  =  1  GOTO  20 

LOCATE  1,7;  INPUT  'enter  the  calibration  data  filename';  filenames 

'  read  in  previously  collected  data  from  a  data  file 

OPEN  'e;\8teve\cal\*  +  filenames  +  '.dat*  FOR  INPUT  AS  #1 
i  =  0 

WHILE  NOTEOF{1) 
i  =  i  +  1 

INPUT  #1,  dummy,  V(i),  E(1,  i),  E(2,  Q 
WEND 
CLOSE  #1 
NPOINTS  =  i 

20  OPEN  'e;\steve\cal\'  +  filenames  +  '.cal'  FOR  OUTPUT  AS  #2 

CLS  ;  LOCATE  1,7;  PRINT  USING  there  are  ##  data  points';  NPOINTS 
PRINT 

'  remove  the  IFA  gain  and  offset 

gain  =  101;  offset  =  1! 

FOR  i  =  1  TO  NPOINTS 

tempi  =  E(1,  i):  temp2  =  E(2,  i) 

E{1 , 1)  =  E(1 , 1)  /  gain  +  offset 
E(2,  i)  =  E(2,  i)  /  gain  +  offset 
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LOCATEi  +  1.7 

PRINT  USING  •##.##  ##.###  ##.###  ##.###  ##.###•;  V(i):  tempi:  E(1.  i):  temp2:  E(2.  i) 
PRINT  #2.  USING  •##.##  ##.###  ##.###  ##.###  ##.###■;  V(i):  tempi;  E(1.  i):  temp2:  E{2. 0 
NEXTr 

INPUT  'press  ENTER  to  continue';  ans$ 


PI  =  3.1416 
mu  =  (PI  /  2!)  'radians 
a(1)  =  mu/2! 
a(2)  =  a(1) 


'  calculate  the  XWIRE  effective  velocities 


2)  +  ((.15  • 
2) +  ((.15* 


SIN 


^  2) 


SIN(a(2)))  -  2) 


cl  =  (COS(a(1 
c2  =  (COS(a  2] 

PRINT  cl,  c2 
FOR  i  =  1  TO  NPOINTS 
ueff(1,i)  =  cl  *V(i)  ''  2 
ueff(2,i)  =  c2*V(i)  ^2 
ueff(1,0  =  SQR(ueff(1.i)) 
ueff(2, 0  =  SQR(ueff|2,  i)) 

PRINT ueff(l.i),  ueff(2.  0 
PRINT  #2.  ueff(1,0.  ueff(2.  i): 

NEXT! 

INPUT :  'hit  enter  key  to  continue',  ans$ 

'  ftt  the  data  to  a  straight  line  using  King's  Lavir 
F0Rj  =  1T02 
FORi  =  lTONPOINTS 
Yq,i)  =  E0.0  2 

XU,  0  =  ueffO,  0  ''  -45 
NEXT! 

SUMX  =  0! 


SUMY  =  0! 

SUMXX  =  0! 

SUMXY  =  0! 

FOR  i  =  1  TO  NPOINTS 
SUMX  =  SUMX  +  Xq,  i) 

SUMY  =  SUMY  +  Yq,  i) 

SUMXX  =  SUMXX  +  Xq.  0  •'  2 
SUMXY  =  SUMXY  +  Xq,  i)  *  Yq,  0 
NEXT! 

B1  q)  =  NPOINTS  ‘  SUMXY  -  SUMX  *  SUMY 
B1 0  =  B1  q)  /  (NPOINTS  *  SUMXX  -  SUMX  ''  2) 
800  =  (SUMY  -  SUMX  *  81 01)  /  NPOINTS 


'  print  the  calibration  data  to  files 


PRINTUSING'80  =  +##.####,81  =  +##.####';  800;  810 
PRINT  #2.  USING  '80  =  +##.####.  81  =  +##.####';  800;  81 0 
NEXT] 

INPUT :  'hit  ENTER  to  continue',  ans$ 

PRINT 

PRINT*  vel  El  exp  El  fit  E20xp  E2rif 

PRINT 

F0Rj  =  1T02 
FORi  =  1TONPOINTS 

Eoutq,  n  =  800  +  81 0)  *  ueffq,  i)  .45 
Eout(j,  i]  =  SQR(EoutO.  i]) 

NEXT! 

NEXT) 

FOR  i  =  1  TO  NPOINTS 

PRINT  V0,  E(1,  0.  Eoutd,  i),  E(2.  0.  Eout{2,  i) 

PRINT  #2.  V(i),  E(1.  i),  Eout(1, 0.  E(2.  i),  Eout(2.  0 
NEXT! 

'  display  results  to  screen  using  the  BASIC  graphics 
INPUT :  'continue  (Y/N) ans$ 

CLS 

SCREENS 

Xmin  =  0:  Xmax  =  639 
Ymin  =  0:  Ymax  =  319 
UNE  (Xmin.  Ymin)-(Xmax,  Ymax),  5, 8F 
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Emax  =  1.2  »  E(1,  NPOINTS) 

Umax  =  1.2  *  ueff(1,  NPOINTS) 

F0Rj  =  1T02  . 

x1  =  0^y1  =  0 

x1  =  639  *  uefffi,  1)  /  Umax 

y1  =  (1!  -  EoutO.  1)  /  Emax)  •  319 

FOR  i  =  2TO  NPOINTS 

x2  =  639  *  uefffi.  0  /  Umax 
y2  =  (1!  -  Eout(j,  i)  /  Emax)  *  319 
UNE{x1,y1)-(x2,y2) 

v1  •=  yO'  \i1  =  x/O 


nbox  =  2 

FOR  i  =  1  TO  NPOINTS 

x2  =  639  •  ueffO,  n  /  Umax 
y2=  {1!-E(j.i) /Emax)  *319 
x1  =  x2  -  nbox:  x2  =  x2  +  nbox 
y1  =  y2  -  nbox:  y2  =  y2  +  nbox 
lJNE(x1.y1)-(x2,y2),2.BF 
NEXT! 

NEXn 

INPUT ‘press  ENTER  to  continue*;  ans$ 

SCREEN  0:  CLS 

'  output  the  calibration  constants  to  the  data  file 

OPEN  •e:\steve\ca^•  +  filenames  +  ‘.cof  FOR  OUTPUT  AS  #1 
FOR  i  =  1  TO  2 

PRINT  #1.  BO (i),B1(0 
NEXTi 
CLOSE  #1 
CLOSE  #2 
END  SUB 


SUB  E1E2  (VQ.  EO.  NPOINTS,  SWITCH,  filenames) 

'  This  subroutane  performs  the  voltage  scan  on  the  hot  wire  for 
'  conversion  to  velocities.  Two  channels  are  scanned  consecutively 
'  with  the  HP44702B  and  the  FET  MuHiplexer.  The  FET  is  connected 
‘  internally  to  the  HP44702B  via  a  ribbon  cable. 

CLS 

DEFINTI-M 

'  initialize  the  HP3852 

CALL  IBFINDrhp3852‘.  dvm%) 

CALL  IBWRT(dvm%,  ‘RST;  CLR") 
rdS  =  3PACES(16) 

'  measure  the  calibration  tank  temperature 

LOCATE  7, 7:  PRINT  'measuring  calibration  tank  air  temperature* 

CALL  IBWRTfdvmX,  'USE  60C:RST  600;REAL  TT;  A2ERO  ONCE*) 

CALL  IBWRT(dvm%,  'CONFMEAS  TEMPJ,  019,  INTO  TT) 

CALL  IBWRT(dvm%.  'VREAD  TT) 

CALL  IBWRT(dvm%,  ’disp  on*) 

CALLIBRD(dvm%,rdS) 

TTC  =  VAL(rdS) 

1TF  =  TTC*{9!/5!)  +32! 

LOCATE  9, 7 

PRINT  USING  the  tank  temp,  is  ##.##  degC  or  degF*:  TTC:  TTF 

'  input  the  probe  and  file  information 

LOCATE  11,7:  INPUT  'enter  the  probe  serial  #*;  snS 

LOCATE  1 3, 7:  INPUT  'enter  the  filename  for  data  output*:  filenames 

caKileS  =  'c:\steve\caA*  +  filenames  +  *.dat" 

LOCATE  15, 7:  PRINT  USING  the  calibration  data  for  probe  &  is  in  &*:  snS:  caifileS 

OPEN  calfileS  FOR  OUTPUT  AS  #1 

LOCATE  19, 7:  INPUT  'press  ENTER  to  continue':  ansS 

CLS 

'  compute  current  air  density,  rho! 
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rholowl  =  1.3947:  drhol  =  .4666:  TIowl  =  23.15:  dTomp!  =  100! 
rhol  =  rholowl  •  drhol  *  (TT  +  TIowl)  /  dTemp! 

'  main  program 

NPOINTS  ==  0 

10  PRINT 'Press  any  k«y  to  begin  scan  ■ 

DO:  LOOP  WHILE  INKEY$  =  " 

NPOINTS  =  NPOINTS  +  1 

PRINT  •  Scanning  hot  wire  voltages.' 

'  proceed  with  channel  1  and  channel  2  voltage  scan. 
nminua1$  =  'REAL  WAVE(1 999),  CHAN1  (999),  CHAN2(999)' 
nscansS  =  'PRESCAN  1000;  POSTSCAN  O' 
nsamplesS  =  'XRDQS  400.  2000  INTO  WAVE' 

CLS  :  LOCATE  7, 7:  PRINT  'voltage  scan  in  progress' 

'  read  hotwires 

CALL  IBWRT(dvm%,  'RST;  CLR;  RST  400') 

CALL  IBWRT(dvm%,  'DISP  OFF') 

CALL  IBWRT(dvm%,  'USE  400') 

CALL  iBWRT(dvm%,  nminusIS) 

CALL  IBWRT(dvm%,  'REAL  LI  .L2,H1  ,H2,AV1  ,AV2,S1  ,S2') 

'  configure  the  HP44702B  for  high-speed  scanning, 

CALL  iBWRT{dvm%,  'SCANMODE  ONiCONF  DCVjARMODE  BEFORE') 

CALL  IBWRT(dvm%,  TERM  RIBBON;range  9:RDGSMODE  COMPLETE') 

CALL  iBWRT(dvm%,  nscansS) 

'  set  the  time  (in  seconds)  between  samples 

sperS  =  'SPER  0.001' 

CALL  IBWRT(dvm%,  spetS) 

CALL  IBWRT(dvm%,  'STTRIG  INF) 

'  scan  the  X-wires 

CALL  IBWRT(dvm%,  'CLWRfTE  SENSE, 321 -322:ASCAN  ON:AZERO  ONCE:SCTRIG  INF) 
CALL  IBWRT(dvm%,  nsamplesS) 

'  dump  the  scanned  voltages  from  the  voltmeter  buffer  to  the  HP3852  mainframe 

CALL  IBWRT(dvm%,  'SUB  SEPARAF)  'Separate  the  readings  into  2  arrays 
CALL  IBWRT  dvm%,  'INTEGER  l,J') 

CALLIBWRT  dvm%,  'J  =  O') 

CALL  IBWRT(dvm%,  'FOR  I  =  0  TO  1999  STEP  2') 

CALL  IBWRT  dvm%,  'CHANI  (J)  =  WAVE(I)') 

CALLIBWRT(dvm%.  'CHAN2(J)  =  WAVE(l+1)') 

CALL  IBWRT(dvm%.  'J  =  J+1') 

CALL  IBWRT(dvm%,  'NEXT  H 
CALL  IBWRT(dvm%.  'SUBEND') 

CALL  IBWRT(dvm%,  'CALL  SEPARATE 
CALL  IBWRT  dvm%,  'stat  LI , HI , AVI, SI. CHANI ') 

CALL  IBWRT(dvm%.  'stat  L2.H2,AV2,S2.CHAN2') 
rd$  =  SPACE$(16) 

CALL  IBWRT(dvm%,  \read  AVI ') 

CALLIBRD(dvm%.  rd$) 

E(1,  NPOINTS)  =  VAL(rd$) 

CALL  IBWRT(dvm%,  Vread  AV2') 

CALL  IBRD(dvm%.  rd$) 

E(2.  NPOINTS)  =  VALfrdS) 

CALL  IBV;RT(dvm%.  'SCTRIG  HOLD') 

CALL  IBWRT(d-vm%.  'DISP  ON;  RST;  CLR') 

'  Comoute  approximate  velocity  at  orifice 

PRINT  .  INPUT 'input  current  tank  pressure  (in  H20)':  Ptank! 

V(NPOINTS)  =  SQR(Ptank!  •  246.7  •  2!  /  rho!) 

PRINT :  PRINT  'current  air  speed  is: ',  V(NPOINrS), '  m/sec' 

PRINT:  PRINT  Ptank,  V(NPOINTS).  E(1,  NPOINTS),  E(2.  NPOINTS) 

PRINT  #1 .  Ptank!,  V(NPOINTS),  E(1 .  NPOINTS).  E(2.  NPOINTS) 

PRINT :  INPUT  'quit?  (y  or  n) '.  q$ 

IF  qSoyTHEN  GOTO  10 

CLOSE  #1 
SWITCH  =  1 

END  SUB 
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2.  Data  Acquisition  Software 
ACQUIRE.BAS  -  main  program 


■ACQUIRE.BAS 

’  This  QuickBASIC  program  performs  automated  data  acquisition  for  the 
'  Turbine  Cascade  Test  Facility  in  Bldg  19  through  the  HP3852A. 

'  To  execute,  at  the  c:\stevo  prompt  enter  QB/L  QBIB4  /AH. 

'  This  will  access  the  GPIB  library  and  allow  large  data  matrix 
’  manipulation  for  turbulence  data  acquisition. 

'  user  written  subroutines 
DECURE  SUB  TEMPACQ  (dvm%,  file$) 

DECURE  SUB  PRESSACQ  (dvm%,  Ttunnel,  file$) 

DECURE  SUB  VEUCQ  (dvm%,  fileS) 

DECURE  SUB  PLOTT  (xdatQ  AS  SINGLE,  ydatQ  AS  SINGLE,  xmini,  xmax!,  ymini,  ymaxi,  xlab$,  ylabS,  pItlabS, 

ntot%) 

DECURE  SUB  TEMPS  (dvm%,  Troom,  Ttunnel) 

DECURE  SUB  VSET  (dvm%,  Ttunnel) 

DECURE  SUB  TSET  (dvm%) 

'  read  in  the  thermocouple  and  pressure  port  chord  locations  for  plotting  purposes 

DIM  chord(24)  AS  SINGLE  'chord  locations  are  common 
OPEN  ‘o:\steve\dat\chord.dat*  FOR  INPUT  AS  #1 
FOR  I  =  1  TO  23 
INPUT  #1,chord(i) 

NEXT  I 
CLOSE  #1 

'  Initialize  the  HP3852  data  acquisition  unit 

CALL  IBFINDfhpSSSS',  dvm%) 

CALL  IBWRT(dvm%,  ‘RST;  CLR') 

'  initialize  the  screen 

CLS  :  COLOR  3,  0:  LOCATE  2, 12 

PRINT  •  HP3852  Controller  Program* 

LOCATE  4, 12 

PRINT  •  Lt  Galassi  6/89,  Capt  Acree  7/90,  Capt  Meschwitz  7/91  * 

LOCATE  6, 12 

PRINT  This  program  controls  the  HP3852  Data  Acquisition  System* 

LOCATE  12, 1 

CALLTEMPS(dvm%,  Troom,  Ttunnel) 

PRINT 

'  input  the  data  run  filename 

INPUT  'Run  Designation  for  all  datafiles  ?  (enter  run  #):*,  file$ 

'  routine  Branching  for  the  MAIN  MENU 
DO 

CLS:  PRINT:  PRINT*  MAIN  MENU*:  PRINT 
PRINT  STRING$(80,  *_*):  PRINT 
°RINT  *1)  Temperature  Run* 

PRINT  *2)  Pressure  Run* 

’■RINT  *3)  Turbulence  Run* 

PRINT  *4)  Position  X-wire* 

PRINT  *5)  Monitor  Blade  Temperature* 

PRINT  *6)  Change  File  Designator* 

PRINT  *^  EXIT* 

PRINT :  PRINT  USING  *the  current  file  designator  is  &*;  file$ 

PRINT :  CALLTEMPS(dvm%,  Troom,  Ttunnel) 

PRINT :  PRINT  Typo  your  selection  (1-7):* 
chs$  =  INPUT$(1) 

'  Use  SELECT  CASE  to  process  response. 

SELECT  CASE  chs$ 

CASE  *1*:  CALL  TEMPACQ(dvm%,  fileS) 

CASE  *2*:  CALL  PRESSACQ(dvm%,  Ttunnel,  fileS) 

CASE  *3*:  CALL  VEUCQ(dvm%,  fileS) 

CASE  *4*:  CALL  VSET(dvm%,  Ttunnel) 

CASE  *5*:  CALLTSET|dvm%) 

CASE  *6*:  CLS  :  INPUT  *New  file  designation:*,  fileS: 

CASE  *7*:  CLS  :  EXIT  DO 


CASE  ELSE 
END  SELECT 
LOOP 


END 


Subroutine  PLOTT 


SUffPLOTT  (xdatO  AS  SINGLE,  ydatQ  AS  SINGLE,  xmlnl,  xmaxl,  ymini,  ymaxl,  xlab$,  ylab$,  pKlab$,  ntot%) 

'  this  section  of  the  code  creates  a  plot  of  ydat  vs  xdat 

CLS 

SCREEN  9  'Hl-res  graphics  mode 
VIEW  (120. 10)-(570. 290). .  1 
WINDOW  (xminl,  yminl)-(xmaxl,  ymaxl) 
style%  =  &HFF00 
LOCATE  23,  42:  PRINT  pttlabS 

LOCATE  1.9:  PRINT  ymaxl:  LOCATE  10, 7:  PRINT  ylab$:  LOCATE  21. 9:  PRINT  yrnini 
LOCATE  22, 14:  PRINT  xminl:  LOCATE  22, 43:  PRINT  xlab$:  LOCATE  22.  71 :  PRINT  xmax! 

VIEW  PRINT  24  TO  25: 

PRINT  SPC(14); ,  'press  any  key  to  return  to  Main  Menu' 

CLS 

LINE  (1.  0):(0.  0), . .  8tyl0% 

FOR  i  =  1  TO  ntot% 

X  =  xdat(i) 

Y  =  ydatli) 

LINE-{X,Y) 

NEXT! 

DO:  LOOP  WHILE  INKEYS  =  " 

CLS 

SCREEN  0:  COLOR  3, 0 
END  SUB 


Subroutine  PRESSACQ 


SUB  PRESSACQ  (dvm%,  Ttunnel,  fileS) 

'  This  subroutine  performs  a  pressure  scan  on  the  turbine  biade 
'  using  the  HP3852A  data  acquisition  system  and  low  speed  voltmeter. 
'  The  voltmeter  takes  and  averages  50  readings  per  port 

DiM  p8can(28),  ydat!(28).  xdat!(28) 

SHARED  chordO  AS  SiNGLE 
DEFihfTI-M 


’  initialize  the  screen 
CLS 

LOCATE  8, 7:  PRINT  'BLADE  PRESSURE  SURVEY  (approx  115  sec.  duration)' 
’  set  up  the  HP44701  for  voHage  measurements 
CALL  iBWRT(dvm%.  'RST  600') 

CALL  iBWRT(dvm%,  'REAL  PSCAN(28),C(49),L,H,M,S') 

’  perform  the  actual  pressure  scan 

CALL  IBWRT(dvm%,  'USE  600:AZERO  ONCE;  NPLC  16') 

CALL  IBWRT(dvm%,  'CLOSE  101')  ’call  homes  the  scanivalve 

CALLIBWRT(dvm%,  'OPEN  101;CONF  DCV) 

CALL  IBWRT  dvm%.  'NRDGS  50') 

CALL  IBWRT(dvm%,  'MEAS  DCV,  221  .INTO  C)  ’read  atmospheric  pressure 
CALL  iBWRT(dvm%,  'STAT  L.H,M,S,C') 

CALLiBWRT(dvm%.  'PSCAN(0)  =  M') 

CALLiBWRT  dvm%,  'DISP  PSCAN(0)') 

CALL  iBWRT(dvm%,  'SUB  LOOP')  ’HP  3852A  subroutine 

CALL  iBWRT(dvm%.  'iNTEGER  n 
CALL  iBWRT  dvm%,  'FOR  I  =  1  TO  28') 

CALL  iBWRT(dvm%,  'CLOSE  100') 

CALL  iBWRT(dvm%,  'OPEN  100  :CONF  DCV:  NRDGS  50') 

CALL  IBWRT(dvm%.  'MEAS  DCV.  221,  USE  600,  INTO  C) 

CALL  IBWRT(dvm%,  'STAT  L.H.M.S.C) 

CALLIBWRT(dvm%,  'PSCAN(I)  =  M') 

CALL  IBWRT  dvm%,  'DISP  PSCAN(I)') 

CALL  IBWRT(dvm%.  'NEXT  ^ 

CALL  IBWRT(dvm%,  'SUBEND') 

CALL  IBWRT(clvm%,  'CALL  LOOP') 

SLEEP  1 20  ’allow  1 20  seconds  for  scan 

CALL  IBWRT(dvm%.  'VREAD  PSCAN') 

’  write  pressure  voltage  readings  onto  the  data  file  on  the  hard  disk 
outfileS  =  •c:\8teve\pre8sure\'  +  file$  +  'pv.daf 
CALL  IBRDF(dvm%,  outfileS) 

CALL  IBWRT(dvm%,  'RST;  CLR') 

’  write  the  tunnel  temperature  to  the  pressure  data  file 
OPEN  outfileS  FOR  APPEND  AS  #1 
PRINT  #1,  Ttunnel 
CLOSE  #1 

’  read  in  the  file  for  display 

OPEN  outfileS  FOR  INPUT  AS  #1 
FORi  =  0TO28 
INPUT  #1,  pscan(i) 

NEXT! 

CLOSE  #1 

’  calculate  the  tunnel  q 
q  =  p8can(26)  -  p8can(28) 

’  perform  shift  on  PSCANQ  to  collapse  the  vector. 

’  this  technique  is  recommended  by  the  Scanivalve  manufacturer  to  help 
’  identify  step  motor  miscues. 

’  array  locations  are  assigned  as  follows: 

’  0,1 1 ,22,29-36  -  Atmospheric  pressure 
’  1-10,12-21,23-25  -  Blade  static  pressures 
’  26  -PHot  head  press. 

’  28-  Pitot  static  press. 


FOR  i  =  11  TO  20 

=  p8can(i  +  1) 

FOR  i  =  21  TO  23 
pscan(i)  =  pscan(i  +  2) 

NEXT  i 

'  this  section  drives  the  piotter. 

ntot%  =  23  This  no.  co'.'esponds  to  the  total  no.  of  pressure  ports 
yiabS  =  'Cp* 
xiabS  =  •x/c’ 

pitiab$  =  'Cp  on  the  turbine  biade' 
xmin!  =  Oi 
xmaxi  =  1 1 
ymin!  =  -5 
ymax!  =  5 
FOR  i  =  1  TO  ntot% 
ydat!(0  =  (pscan(i)  -  pscan(28))  /  q 
xdatl(i)  =  chord(i) 

NEXTi 

CALL  PLOTT(xdat!0,  ydatlQ,  xmin!,  xmaxi,  ymin!,  ymax!,  x!ab$,  y!ab$,  pH!ab$,  ntot%) 
END  SUB 
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Subroutine  TEMPACQ 


DEFSNG  l-M 

SUB  TEMPACQ  (dvm%,  fileS) 

'  This  subroutine  performs  the  temperature  scan  on  the  turbine  biade 
'  and  stores  the  data  in  c:\stove\temp\file$.dat  The  HP3852  is  used  to 
'  acquire  data  through  the  HP  4471 3A  FET  Muitiplexers. 

'  The  program  momentarily  interrupts  the  power  to  the  turbine  blade 
'  with  relay  channel  104  so  as  not  to  overload  the  thermocouples 
SHARED  chordO  AS  SINGLE 
DIM  temp(30),  xdat(30),  ydat(30) 

OEFINT  l-M 

'  read  the  blade  temperatures,  and  Ttunnel  and  Troom. 

'the  first  ISTcouples  (wire  #8  1-16)  are  on  slot  2  (Channels  201-216), 

'  and  the  next  14  (wire  #8  17-30)  are  on  slot  0  (Channels  007-020). 

CLS 

'  configure  the  HP3852  and  measure  the  temperatures 

LOCATE  7, 7:  PRINT  'measuring  temperatures' 

CALL  IBWRT(dvm%,  'RST600:  USE  600') 

CALL  IBWRT(dvm%,  'AZERO  ONCE') 

CALL  IBWRT(dvm%,  'REAL  A(29)') 

CALL  IBWRT(d-vm%,  'DISP  OFF') 

CALL  IBWRT(dvm%,  'CONF  TEMPJ') 

CALL  IBWRT(dvm%.  'CLOSE  104')  'blade  current  off 

CALL  IBWRT(dvm%,  'MEAS  TEMPJ,  201-216,  007  -  020,  INTO  A') 

CALL  IBWRT(dvm%,  'OPEN  104')  'blade  current  on 

CALL  IBWRT(dvm%.  'VREAD  A') 

'  transfer  the  readings  to  the  PC  for  temporary  storage 

LOCATE  9, 7:  PRINT  transferring  data  to  PC  hard  disk  storage' 

CALL  iBRDF(dvm%,  'c:\8teve\temp\bladotem.daf) 

CALL  IBWRT(dvm%,  'DISP  ON:  RST;  CLR') 

'  create  permanent  storage  and  display  readings  to  the  screen 

LOCATE  11,7 

PRINT  'creating  permanent  storage' 

OPEN  •c:\steve\temp\bladetem.daf  FOR  INPUT  AS  #1 
OPEN  'o:\steve\temp\'  +  fileS  +  T.daf  FOR  OLfTPUT  AS  #2 
LOCATE  13, 7:  PRINT '  TC#  Temp' 

LOCATE  1 4,  7:  PRINT '  (degC)' 

FOR  I  =  1  TO  30 
INPUT  #1,temp(i) 

PRINT  #2,  temp(i) 

NEXT! 

CLOSE  #1 
CLOSE  #2 

FORi=1TO30STEP2 

PRINT  USING'##  ###.###  ##  ##.###':  I;  temp(i):  I  +  1:  tempo  +  1) 

NEXT! 

INPUT '  press  ENTER  to  continue';  rtn$ 

'  create  a  file  for  'GRAPHER' 

OPEN  'c:\steve\plt\'  +  fileS  +  TGR.daf  FOR  OUTPLR'  AS  #1 
FOR  I  =  1  TO  23 

PRINT  #1,  USING  '#.#######  ##.#####':  chord(i):  temp(i) 

NEXT! 

CLOSE  #1 

'  plot  the  temperature  readings 

ntot1%  =  8  This  no.  is  the  total  no.  of  pressure  side  thermocouples 
ntot2%  =  12  This  no.  is  to  the  total  no.  of  suction  side  thermocouples 
ylabS  =  Temp.  (C)':  xIabS  =  'x/C 
phlabS  =  Temp,  on  turbine  blade  ‘ 
xmin!  =  01:  xmaxi  =  1 !:  ymin!  =  20:  ymax!  =  65 
FOR  I  =  1  TO  ntot1% 
ydatin  =  temp(l) 
xdat(i)  =  chordfi) 

NEXT! 
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FOR  i  =  ntot1%  +  1  TO  ntot2% 
xdat(i)  =  chord(n 
ydat(i)  =  0 

N0ai 

FOR  i  =  ntot2%  +  1  TO  ntot1%  +  ntot2%  +  3 
ydat(i)  =  temp(i  -  4) 
xdat(i)  =  chord(i) 

NEXT! 

ntot%  =  ntot1%  +  ntot2%  +  3 

CALL  PLOTT(xdat!0,  ydatlQ,  xmin!,  xmax!,  yminl,  ymax!,  xlab$,  ylab$,  pltlab$,  ntot%) 
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Subroutine  TEMPS 


DEFSNG  I'M 

SUB  TEMPS  (dvm%,  Troom,  Ttunnel) 

'  This  subprogram  configures  the  HP3852  to  read  the  room  temperature 
'  on  Channel  012  (wire  29)  and  the  tunnel  temperature  on  Channel  020  (wire  30) 

CALL  IBWRT(dvm%.  'USEeoOiRSTeoOiREAL  Wm.TtunneLAZERO  ONCE") 

'  read  the  room  temperature 

CALL  IBWRT(dvm%,  ‘CONFMEAS  TEMPJ,  019,  INTO  Troom’) 

CALL  IBWRT(dvm%,  ’DISP  Troom’) 

CALL  IBWRTtdvm*,  ’VREAO  Troom’) 
rd$  =  SPACE$(16) 

CALLIBRD(dvm%,rd$) 

Troom  =  VAL(rd$)  'convert  the  string  value  into  a  number 
'  read  the  tunnel  temperature 

CALL  IBWRT(dvm94.  ’CONFMEAS  TEMPJ,  020,  INTO  Ttunnel’) 

CALL  IBWRT  dvm%,  ’DISP  Ttunnel’) 

CALL  IBWRT(dvm96,  ’VREAD  TtunneH 
rd$  =  SPACE$(16) 

CALL  IBRD(dvm%,  rd$) 

Ttunnel  =  vAL(rd$)  'convert  the  string  value  into  a  number 

PRINT  USING  the  room  temp  is:  ##.##  degF  (##.##  degC)’;  (Troom  •  9!  /  5!)  +  32!;  Troom 
PRINT  USING  the  tunnel  temp  is:  ##.##  degF  (##.##  degC)’;  (Ttunnel  •  9!  /  5!)  +  321;  Ttunnel 
CALL  IBWRT(dvm%,  ’RST;  CLR') 

END  SUB 
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Subroutine  VSET 


SUBTSET(dvm%) 

'  This  subprogram  screon  displays,  in  bar  graph  form,  the  room  temperature 
'  and  the  5  internal  blade  temperatures 

DIM  temp(6) 

CLS 

tempO$  =  SPACE$(16) 
tempts  =  SPACE$(1 6) 
temp2$  =  SPACE$(16) 
temp3$  =  SPACE$(16) 
temp4S  =  SPACE$(16) 
temp5$  =  SPACE$(16) 

’  set  the  screen  graphics 

SCREEN  9 

ymin  =  Ohymax  =  1001 
UNE  (0, 0H639.  319),  3,  BF 
FORi  =  0TO10 
Y=  (1  .i/10)*319 
LINE  (0,  Y)-(639,  Y) 

NEXT! 


’  outline  the  T  =  60  degC  line 

yset  =  601 

ytol  =  .5 

yO  =  yset  -  ytol 

y1  =  yset  +  ytol 

yO  =  (1 1 .  (yO  -  y min)  /  (ymax  -  y min))  *  31 9 
y1  =  (1 1  •  &1  •  ymin)  /  Nmax  -  ymin))  *  319 
LlNE(0,y0M639,y1),6,  BF 
yset  =  (1  •  (yset  -  ymin)  /  (ymax  -  ymin))  *  319 
LINE  (0,  ys^-(639,  yset) 

'  initialize  the  HP 


CALL  IBWRT(dvm%,  'RSTOOO:  USE  OGOT 
CALL  IBWRT  dvm%,  'AZERO  ONCET 
CALL  IBWRT  dvm%,  'DISP  OFF’) 

DO 

'  take  the  temperature  measurements 
CALL  IBWRT(dvm%,  ’REAL  A(5)*) 

CALL  IBWRT(dvm%,  ’CONFMEAS  TEMPJ,  020, 014  •  018,  USE  600  INTO  A*) 
'  transfer  the  measurements  to  the  PC 


CALL  IBWRT(dvm%,  -VREAD  A(0)‘) 
CALL  IBRD(dvm%,  tempOS) 
tomp(1)  =  VAL(tempO$) 

CALL  IBWRT(dvm%,  “VREAD  A(l)’) 
CALL  IBRD(dvm%,  tempi  $) 
temp(2)  =  VAL(temp1$) 

CALL  IBWRT(dvm%,  “VREAD  A(2)') 
CALL  IBRD(dvm%,  temp2$) 
temp(3)  =  VAL(tomp2$) 

CALL  IBWRT(dvm%,  'VREAD  A(3)') 
CALL  IBRD(dvm%,  temp3$) 
temp(4)  =  VAL(tomp3$) 

CALL  IBWRT(dvm%,  'VREAD  A(4)') 
CALL  IBRD(dvm%,  temp4$) 
temp(5)  =  WL(tomp4S) 

CALL  IBWRT(dvm%,  'VREAD  A(5)') 
CALL  IBRD(dvm%,  temp5$) 
temp(6)  =  VAL(temp5$) 


'  plot  the  measurements  on  the  PC  graphics  screen 


nplot  =  6 

ymin  =  0!:  ymax  =  100! 
xmax  =  nplot 
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nbox  =20 
yO  =  0 

FOR  i  =  1  TO  nplot 
71  =  (1  -temp{i)  /ymax)  *319 
Xl  =  {(D  /  (xmax  +  1))  *  639 
xO  =  x1  -  nbox:  x1  =  x1  +  nbox 
UNE(xO.  y0)-(x1,y1),3,  BF 
NEXT! 
yO  =  319 

FOR  i  •■=  1  TO  nplot 
y1  =  (1  -  temp(l)  /  ymax)  *  31 9 
x1  =  ((0  /  (xmax  +  1))  *  639 
xO  =  x1  -  nbox:  x1  =  x1  +  nbox 
UNE{xO.y0).(x1,y1).5,  BF 
NEXT! 

LOOP  UNTIL  INKEYS  <>•• 

SCREEN  0:  COLOR  3,0 
CALL  IBWRT(dvm%,  'CLR;  RSP) 
END  SUB 
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Subroutine  VELACQ 


SUB  VELACQ  (dvm%.  (i!e$) 

DIM  El  (16384),  E2(1 6384) 

'  This  subroutine  performs  the  voKage  scan  on  the  hot  wire  for 
'  conversion  to  velocities.  Two  channels  are  scanned  consecutively 
'  with  the  HP44702B  and  the  FET  Multiplexer.  The  FET  is  connected 
’  internally  to  the  HP44702B  (high  speed  voltmeter)  via  a  ribbon  cable. 

DEFIMTI-M 

'  set  the  number  of  scans  per  channel 
nscans  =  16384 
N  =  nscans 
nsamples  =  2  *  nscans 
nminusi  =  nsamples  - 1 
nscans$  =  STR$(nscans) 
nsamplesS  =  STR$(nsamples) 
nminus1$  =  STR$(nminus1) 

nminusIS  =  'REAL  WAVEC  +  nminusIS  +  •).  OUr(8191)* 
nscansS  =  'PRESCAN  •  +  nscansS  +  *;  POSTSCAN  O’ 
nsamplesS  =  ‘XRDGS  400,  ’  +  nsamplesS  +  *  INTO  WAVE’ 

CIS  :  LOCATE  2, 7;  PRINT  ’  Turbulence  scan  in  progress’ 

'  configure  the  HP3852 

CALL  IBWRT(dvm%,  ’CLR;  RSr) 

CALL  IBWRT  dvm%,  ’RST  400;  USE  400’) 

CALL  IBWRT(dvm%,  nminusIS) 

CALL  IBWRT(dvm%,  ’DISP  OFF’) 

CALL  IBWRT(dvm%,  ’SCANMODE  ON;  CONF  DCV;  ARMODE  BEFORE’) 

CALL  IBWRT(dvm%,  TERM  RIBBON;  range  9;  RDGSMODE  COMPLETE’) 

CALL  IBWRT(dvm%,  nscansS) 

CALL  IBWRT(dvmX,  ’STTRIG INT) 

‘  set  the  sample  period  s  1  /  fsample 
period  =  .00005 
periods  =  STR$(period) 
sperS  =  ’SPER  ‘  +  periods 
CALL  IBWRT(dvm%,  spetS) 

'trigger the  scans 

LOCATE  4, 7:  PRINT  ’Acquiring  Data* 

CALL  IBWRT(dvm%,  ’CLWRiTE  SENSE,321-322:  ASCAN  ON;  AZERO  ONCE;  SCTRIG  INT) 

'  transfer  the  data  from  the  44702B  VM  to  the  hard  disk  via  HP  mainframe 
FORi  =  0TO3 

LOCATE  6, 7:  PRINT  USING  ‘transfering  ####  data  points  to  PC  file  FREQ#.dat';  N/2;  i 
CALL  IBWRT(dvm%,  ’XRDGS  400.  8192  INTO  OUT) 

CAU  IBWRT(dvm%,  -VREAD  OUT) 

filenames  =  •c:\8teve\dat\FREQ’  +  LTRIM$(STRS(0)  +  ’.DAT 

CALL  !BRDF(dvm%,  filenames) 

NEXT! 

CALL  IBWRT(dvm%.  ’DISP  ON;RST:  CLR’) 

'  read  in  hot  wire  data  from  temporary  storage 

k  =  0 

FORi  =  0TO3 

infileS  =  ’ciWeveWaftfreq’  +  LTRIM$(STRSO))  +  ’.DAr 
LOCATE  8, 7:  PRINT  USING  ’reading  in  rawX-wire  data  from  infileS 

OPEN  infileS  FOR  INPUT  AS  #1 
FORi  =  1  TON/4 
k  =  k  +  1 

INPUT  #1.  El  (k).E2(k) 

NEXT! 

CLOSE  #1 
NEXT) 

'  write  anemometer  voltages  to  hard  disk  in  binary  files 


F-17 


e1file$  =  'c:\8teve\daA'  +  file$  +  'el.dat* 
e^le$  =  'c:\8teve\dal\*  +  file$  +  •e2.dat‘ 

LOCATE  10,  7:  PRINT  'writing  data  in  binary  filea' 

LOOVTE  12,  7:  PRINT  USING  "writing  files  &  and  elfileS;  e2file$ 

OPEN  elfileS  FOR  BINARY  AS  #1 
OPEN  62file$  FOR  BINARY  AS  #2 

'  remove  the  IFA  gain  and  offset 
gain  =  101;  offset  =  1! 

FOR i  =  1  TO  N 
El  (n  =  El  fn  /  gain  +  offset 
E2(i)  =  E2(i)  /  gain  +  offset 
PUT#1,,E1(0 
PUT#2,,E2(i) 

NEXTi 
CLOSE  #1 
CLOSE  #2 

’  print  the  first  5  and  last  5  entries  for  El  and  E2 
CLS 

LOCATE  3, 7:  PRINT'  El  E2';  PRINT 

FOR  I  =  1  TO  5 

PRINT  USING '#####  ##.####  #####  ##.####■;  i;  El  (i):i:E2(i) 
NEXTi 

FORi  =  N.4TON 

PRINT  USING  ■#####  ##.####  #####  ##.####■;  i:  El  (i):i:E2(i) 
NEXTi 

LOCATE  20, 7;  INPUT 'press  ENTER  to  return  to  the  Main  Menu';  ans$ 

END  SUB 
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Subroutine  VSET 


DEFSNG  l-M 

SUB  VSET  (dvm%,  Ttunnel) 

’  This  program  displays  the  sum  of  the  X-wire  voltages  and  the  U  and  V 
'  components  of  velocity.  Maximizing  the  sum  of  the  voltages  or  zeroing 
'  the  V  component  of  the  velocity  ensures  that  the  bisector  is  aligned 
■  with  the  mean  velocity  vector. 

DIM  B0(2),B1(2) 

'  initialize  the  screen 

COLOR  3, 0 
CLS 

LOCATE  3, 20:  PRIhTT  %WIRE  AUGNMEhfT  PROGRAM’ 

PRINT :  PRINT  The  X-WIRE  bisector  is  aligned  with  the  mean  velocity* 
PRINT  “vector  when  the  sum  of  tho  x-wire  voHages  (below)  is  maximized' 
VIEW  PRINT  8  TO  24 

‘  read  in  the  hot-wire  calibration  constants 


OPEN  ’c:\steve\caI\xwirecal.cof  FOR  INPUT  AS  #1 
FOR  j  =  1  TO  2 
INPUT  #1.B0®,  81® 

NEXT) 

'  compute  current  air  density,  rho! 

rholow!  =  1.3347;  drho!  =  .4666:  Tlow!  =  23.15:  dTemp!  =  100! 
rho!  =  rholow!  -  drho!  •  (Ttunnel  -t-  Tlow!)  /  dTemp! 

'  initialize  the  HP3852 

CAaiBWRT{dvm%.  RST;CLfl’) 

CALL  IBWRT(dvm%.  ’OISP  OFF! 

DO 

'  take  the  x-wire  voltage  readings  (100  readings  per  channeQ 

CALL  IBWRT(dvm%.  ’RST  400*) 

CALL  IBWRT(dvm%,  ’USE  400*) 

CALL  l8WRT(dvm%.  ’REAL  WAVE(1 99),  CHAN1  (99).  CHAN2(99)’) 

CALL  IBWRT(dvm%,  ’REAL  LI  .12.H1  .H2,AV1  >tV2.Sl  ,S2^ 

CAU  IBWRT(dvm%,  ’SCANMODE  ON:CONF  OCV;ARMODE  BEFORE*) 

CALL  IBWRT(dvm%.  NERM  RIBBON;RANGE  9;RDGSMODE  COMPLETE*) 

CAa  IBWRT(dvm%,  ’PRESCAN  100;  POSTSCAN  O’) 

CALL  IBWRT(dvm%.  *SPER  0.001*) 

CALL  IBWRT(dvm%.  ’STTRIG INT) 

CALL  IBWRT(dvm%.  ’CLWRfTE  SENSE,321-322:ASCAN  ON:AZERO  ONCE;SCTRIG  INT) 
CALL  IBWRT(dvm%.  *XR0GS  400, 200  INTO  WAVE*) 

■  separate  the  WAVE  array  into  2  arrays 

CALL  IBWRT(dvm%.  *SUB  SEPARAT*) 

CAa  IBWRT(dvm%.  ’INTEGER  I.J*) 

CAa  IBWRT(dvm%.  *3  =  O’) 

CAa  IBWRT(dvm%.  ’FOR  1  =  0  TO  199  STEP  2*) 

CAa  IBWRT(dvm%.  ’CHANI  (J)  =  WAVE(0’) 

CAa  IBWRT(dvm%.  ’CHAN2(J)  =  WAVE(I+1)’) 

CAa  IBWRT(dvm%.  ’J  =  J+l*) 

CAa  IBWRT(dvm%.  ’NEXT  !■) 

CAa  IBWRT(dvm%,  ’SUBENO-) 

CAa  IBWRT(dvm%.  ’CAa  SEPARAP) 

'  take  the  mean  of  the  100  readings 

CAa  IBWRT(dvm%.  'stat  LI  .HI  3^V1  .SI  .CHANI  *) 

CAa  IBWRT(d-^m%.  ’stat  12,H2.AV2.S2.CHAN2’) 
rd$  =  SPACE$(16) 

CAa  IBWRT(dvm%.  Vread  AVI’) 

CAaiBRD(dvm%.rd$) 

El  =  VAL(rd$) 
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CALL  IBWRT(dvm%,  Vread  AV2*) 
OLLIBRD(dvm%.  fd$) 

E2  =  VAL(rd*) 

CALL  IBWRT(dvm%,  ’SCTRIG  HOLD') 
CALL  IBWRT(dvm%,  'RST;  CLR') 

'  remove  the  IFA100  gain  and  offset 


gain  =  101:  offset  =  1! 
tempEI  El 
tempE2  =  E2 
El  =  El  /  gain  +  off'^^''. 

E2  =  E2/gain  +  oILot 

'  calculate  wire  effective  velocities 

El  =  ((El  2-B0(1))/B1(1))  2.22 

E2  =  ((E2  2  -  B0(2))  /  B1  (2))  2.22 

'  calculate  velocities  in  wire  coordinates,  then  lab  coordinates 

k  =  .15 

PI  =  3.1416 


mu  =  PI/ 4 


X  =  SQR(XX) 

YY  =  'F2  2-  (k  *  El)  ''  2)  /  (1!  -  k  ''  4) 

IFYY<0ITHENYY  =  Y''2 


Y  =■  SQR(YY) 

U  =  COS(mu)  *  X  +  SIN(mu)  *  Y 

V  =  -SIN(mu)  *  X  +  COS(mu)  *  Y 


PRINT  USING 'El +E2  =  ##.####  U  = +##.###  V  =+##.###';  (El  +  E2);  U;  V 


LOOP  UNTIL  INKEYS  <>•• 


CALL  IBWRT(dvm%,  'DISP  ON;  RST;  CLR') 
VIEW  PRINT:  CLS 
END  SUB 


F-20 


3.  Data  Reduction  Software 


PREDUCE.BAS 


■  PREDUCE.BAS 

'  This  program  reads  in  the  pressure  voltage  readings  from  the  data  files 
’  created  by  ACQUIRE.BAS  and  converts  them  into  pressure  values  using  the 
'  calibration  constants  stored  in  PCOEFF.OAT.  This  program  also  creates 
'  a  data  file  for  GRAPHER. 

DIM  pvolts(a8),  pressure(28),  chord(23),  Cp(23t 
'  initialize  screen  and  input  the  pressure  data  fife  to  be  reduced 

CLS 

LOCATE  1.7:  PRINT  'PRESSURE  DATA  REDUCTION  PROGRAM* 

LOCATE  3,  7: 

INPUT  'Enter  the  4  letter  file  designator*;  file$ 

'  read  in  the  pressure  calibration  coefficients 
LOCATE  5,  7:  PRINT  'reading  in  the  pressure  calibration  coefficients' 

OPEN  'o:\steve\cal\PCOEFF.dal*  FOR  INPUT  AS  #1 
INPLrT#1.m.  b 

PRINT  USING  the  calibration  curve  is  P  =  ##.#)!f##V  +  ##.#####';  m;  b 
CLOSE  #1 

INPUT  *  press  ENTER  to  continue';  ans$ 

'  read  in  the  pressure  voltage  readings 
PRINT 

PRINT  USING  'reading  in  the  voHage  readings  from  &pv.daf ;  file$ 

PRINT 

OPEN  *o:\stove\ht)rfeAprassure\*  +  file$  +  'pv.dat*  FOR  INPUT  AS  #1 
FORi  =  0TO28 
INPUT  #1.pvolt8(0 

PRINT  USING  'port  ##  voHage  =  ###.######  volts';  i;  pvolts(i) 

NEXT! 

INPUT  #1,Ttunnel 

PRINT :  PRINT  USING  Ttunnel  =  ###.####  degC;  Ttunnel 
CLOSE  #1 

INPUT  *  press  ENTER  to  continue';  ans$ 

'  convert  the  voltage  readings  into  kPa  pressures 
PRINT :  PRINT  'converting  the  voltages  to  pressures' 

PRINT 

FOR  i  =  0  TO  28 
pressure(i)  =  (m  *  pvotts(i))  +  b 

PRINT  USING '  port  ##  P  =  ###.#####  Kpa';  i;  pressurefi) 

NEXT  i 

INPUT '  press  ENTER  to  continue';  ans$ 

’  write  the  pressures  to  a  data  file  for  permanent  storage  a 

PRINT :  PRINT  *writing  the  pressures  to  &p.dat*  * 

OPEN  'e:\steve\htxfer1pressure\*  +  file$  +  'p.daf  FOR  OUTPUT  AS  #1 

PRINT  #1 ,  USING '  the  calibration  for  this  data  is  P(kPa)  =  ##.#####V(volts)  +  ##.#####)';  m;  b 

FORi  =  0TO28  ? 

PRINT  #1,  USING*  ##  ##.#######  ##.######';  i;  pvortsfO;  pressure^  1 

NEXT!  :! 

'  calculate  the  tunnel  q,  deltah,  and  velocity 
q  =  pre88ure(26)  -  pressure(28) 

rho  =  (101 .3  +  pres8ure(28))  /  (.28708  *  CTtunnol  +  273.1 5)) 
deltah  =  q*1000l/ 2491 
veil  =  SQR(2000I  *  q  /  rho) 

'  print  the  velocity  information  to  a  file  ' 

PRINT  #1," 

PRINT 

PRINT  #1,  USING*  the  tunnel  q  is  ###.#######  kPa';  q 
PRINT  USING*  the  tunnel  q  is  ###.#######  kPa';  q 
PRINT  #1,  USING*  the  tunnel  deltah  is  ###.#######  in  H20';  deltah 
PRINT  USING*  the  tunnel  deltah  is  ###.#######  in  H20';  deltah 
PRINT  #1 ,  USING '  the  tunnel  velocity  is  ###.#######  m/s';  veil 
PRINT  USING*  thetunnel  velocity  is  ###.#######  m/s';  vel! 

CLOSE  #1 

INPUT '  press  ENTER  to  continue';  ans$ 
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'  'collapse  the  vector*  by  eliminating  atmospheric  pressure  readings 
PRINT 

FORi=  1  TO  10 

PRINT  USING  •  tap  ##  P  =  ###.####  Kpa';  i;  pressure^ 

NEXT! 

FOR  i  =  1 1  TO  20 
pressure(i)  =  pressure(i  +  1) 

PRINT  USING  •  tap  ##  P  =  ###.####  Kpa';  i;  pressure(0 
NEXT! 

FOR  I  =  21  TO  23 
pressure(i)  =  prossure(i  +  2) 

PRINT  USING  ■  tap  ##  P  =  ###.####  Kpa';  i;  pressure(i) 

NEXT! 

INPUT  ■  press  ENTER  to  continue';  ans$ 

'  calculate  the  blade  surface  pressure  coefficients  (Cp) 

PRINT :  PRINT  'calculating  Ops* 

FORi  =  1  TO  23 

Cp(i)  =  (pressureO')  -  pressure(28))  /  q 

PRINT  USING  tap  ##  Cp  =  ###.#####*;  i;  Cp(i) 

NEXT  I 

'  wrHing  the  chord  locations  and  the  Ops  to  the  GRAPHER  file 

LOCATE  15, 7:  PRINT  “writing  the  chord  locations  and  Ops  to  &pgr.dat*;  file$ 

OPEN  'o:\8teve\chord.dar  FOR  INPUT  AS  #1 

OPEN  *e;\8teve\htxfeApro88uro\*  +  file$  +  'pgr.dat*  FOR  OUTPUT  AS  #2 
FOR  I  =  1  TO  23 
INPUT  #1,chord(i) 

PRINT  #2,  USING '#.########  ,  ###.########';  chord(i):  Cp(i) 
NEXT! 

CLOSE  #1 
CLOSE  #2 
END 
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TREDUCE.BAS 


’TREDUCE.BAS 

’  This  program  reads  in  the  temperature  readings  from  the  data  files 
'  created  by  ACaUIRE.BAS  and  converts  them  into  chord  representative  values 
DIM  iemp(30),  chord(23) 

'  initialize  the  screen 

CLS 

LOCATE  1,7:  PRINT  TEMPERATURE  DATA  REDUCTION  PROGRAM’ 

LOCATE  3, 7: 

10  INPUT  ’Enter  the  4  letter  file  designator*;  file$ 

'  read  in  the  temperature  readings 
PRINT 

PRINT  USING  'reading  in  the  temperature  readings  from  &&t.dat’;  file$ 

PRINT 

OPEN  ’e:\8teve\htxfeAtemp\’  +  file$  +  ^dat’  FOR  INPUT  AS  #1 
FOR  i  =  1  TO  30 
INPUT  #1,temp(0 

PRINT  USING  tc  ##  temp  =  ###.######  degC’;  i;  temp(i) 

NEXT! 

CLOSE  #1 

INPUT  ’  press  ENTER  to  continue’;  ans$ 

’  writing  the  chord  locations  and  the  temperatures  to  the  GRAPHER  file 
LOCATE  15, 7:  PRINT  “wrrting  the  chord  locations  and  temps  to  &pgr.dat*;  file$ 
OPEN  ’e;\steve\chord.dar  FOR  INPUT  AS  #1 
OPEN  ’e:\steve\htxfeAtemp\’  +  file$  +  tgr.dat*  FOR  OUTPUT  AS  #2 
FOR  i  =  1  TO  23 
INPUT  #1,chord(i) 

PRINT  #2,  USING  •#.########  ,  ###.########’;  chord(i);  tempfi) 
NEXT! 

CLOSE  #1 
CLOSE  #2 
GOTO  10 
END 
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HTXFER.BAS 


HTXFER.BAS 

This  program  calculates  the  blade  heat  transfer 

DECLARE  SUB  PREDUCE  (file$,  Uinf,  CpQ.  pressureQ,  q) 

DIM  visc(23)  AS  SINGLE 
DIM  Tscan(30)  AS  SINGLE 
DIM  detUI(23)  AS  SINGLE 
DIM  deltaT(23)  AS  SINGLE 
DIM  resi8tivity(23)  AS  SINGLE 
DIM  Cp(23)  AS  SINGLE 
DIM  pressure(28)  AS  SINGLE 
DIM  rho(23)  AS  SINGLE 
DIM  mu(23)  AS  SINGLE 
DIM  vel(23)  AS  SINGLE 
DIM  Re(23)  AS  SINGLE 
DIM  8(23)  AS  SINGLE 
DIM  h(23)  AS  SINGLE 
DIM  Nu(23)  AS  SINGLE 
DIM  Nutheorl(23)  AS  SINGLE 
DIM  Nutheort(23)  AS  SINGLE 
DIM  St(23)  AS  SINGLE 
DIM  Sttheorl(23)  AS  SINGLE 
DIM  Sttheort(23)  AS  SINGLE 
DIM  chord(23)  AS  SINGLE 
DIM  K(23)  AS  SINGLE 
DIM  fttpltl(23)  AS  SINGLE 
DIM  fltpttt(23)  AS  SINGLE 

DEFINTI-J 

’  read  in  thermocouple  distances  along  the  blade  surface  SQ  for  local  calculations 
OPEN  'eAsteveWahSNEW.daf  FOR  INPUT  AS  #1 
FOR  I  =  1  TO  23;  INPUT  #1 , 8(1):  NEXT  I 
CLOSE  #1 

'  read  in  the  chord  positions  of  the  thermocouples  and  pressure  taps 

OPEN  •e:\steve\dat\chord.dat*  FOR  INPUT  AS  #1 
FOR  I  =  1  TO  23:  INPUT  #1,  chord(i):  NEXT  i 
CLOSE  #1 

'  read  In  the  deKal  distances  between  the  surface  thermocouples  and  the 
'  internal  thermocouples  for  the  conduction  calculations 
OPEN  •e;\8teve\dat\DELTAL.DAT’  FOR  INPUT  AS  #1 
FORI  =  1  TO  23 

INPUT  #1,dettal(0 

NEXT! 

CLOSE  #1 

'  main  loop 

FOR  alpha  =  1  TO  4 
FOR  tost  =  1  TO  6 

files  =  ‘hx*  +  LTRIMS(STR$(alpha))  +  LTRIM$(STR$(tost)) 

CLS  :  PRINT  USING  'processing  data  file  fileS 
’  read  the  blade  temperatures  and  current  from  storage 

OPEN  ■e:\8tovo\htxfer\temp\'  +  fileS  +  Tdat*  FOR  INPUT  AS  #1 
FORI  =  1  T030 

INPUT  #1.Tscan(i) 

NEXT! 

INPUT  #1.Curr! 

CLOSE  #1 

Tint!  =  T8can(30)  'Use  tunnel  temperature  as  Tfreestream 
'  print  the  blade  temperatures  and  power  to  the  screen 
PRINT  USING  tunnel  temp  =  ###.###  degC;  Tint! 

PRINT:  PRINT'  I  degC  I  dogC 
FORI  =  1  TO  15 

PRINT  USING'  ##  ###.###  ##  ###.###';i:Tscan(i):i+ 15:Tscan(i  + 15) 

NEXT! 

PRINT  :  PRINT  USING  the  blade  current  is  ##.##  amps';  Curr! 

’  INPUT 'press  enter  to  continue';  ansS 
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'  read  in  the  Pressure  Coefficients  on  the  blade  to  be  used  in  computation  of 

CALL  PREDUCE(file$,  Vinfl,  CpQ.  pressureQ.  q) 

'  ConstrucTdettaTli),  the  internal  temp  differences 

FOR  i  =  1  TO  3 

deltaT(i)  =  Tscan(i)  -  T8can(24) 


FOR  i  =  22  TO  23 

deltaT(i)  =  Tscan(i)  -  T8can(24) 

NEXT! 

deltaTM  =  T8can(4)  -  Tscan(25) 
deltaT(5)  =  TscanfS)  -Tscan(26) 

dettaTf^  =  Tscanfffl  -  n'scan(2^  +  Tscan(26)  +  Tscan(15))  /  3! 

deitaTm  =  Tscanm  -  tscan(27) 

dettaT(8)  =  T8can(8)  -  T8can(28) 

dettaTra)  =  T8can(9)  -  T8can(12) 

dettaT(IO)  =  TscanflO)  -T8can(12) 

dettaT(ll)  =  Tscanlll)  -Tscan(28) 

deltaT(12)  =  T8can(12i  -Tscan(9) 

deltaT(13)  =  Tscan(l3|  •Tscan(28) 

deltaT{14)  =  T8can|l4)  -  (T8can(27)  +  T8can(28)  +  Tscan{7))  /3I 
deltaTfIS)  =  TscanfIS)  -Tscan(27) 

deltaTflg  =  Tscan(16)  -  n’scan(26)  +  Tscan{27)  +  Tscan(6))  / 3! 

deltaT(17)  =  TscanflT)  -tscan(26) 

deltaT(18)  =  TscanilS)  -Tscan(26) 

deltaT(19)  =  T8can(l9)  -  .5  *  (rscan(25)  +  Tscan(26)) 

deltaT(20)  =  Tacan(20)  -  .5  *  ^8can(25)  +  Tscan(26)) 

deltaT(21]  =  Tacan(21)  -Tscan(25) 

'  compute  the  local  viscosities  visc(i)  on  the  blade 
'  the  density  rho(i)  determined  from  the  ideal  gas  law  and  viscosity  mu(i) 

'  is  determined  from  its  linearity  over  a  wide  range  of  pressures 
'  all  calculations  are  based  on  Tavg  =  (TIocal  +  Tinf)/2 
FOR  i  =  1  TO  23 

rho(0  =  (101.3  +  pressure(O)  /  (.28708  •  ((rscan{0  +  Tinl)/2  +  273.15)) 
mu(0  =  4.136E-06  +  ((T8can(0  +  Tint)  /  2!  +  273.15)  *  4.76E-08 
vi8c(i)  =  mu(i)  /  rho(i) 

NEXT! 

'  compute  local  velocities  vel(i)  on  the  blade  (m/s) 
f6r  i  =  1  TO  23 

vel(0  =  SQR((2  *  1000  *  q  *  ABS(1!  -  Cp(i)))  /  rho(i)) 

NEXT  i 

'  compute  local  Reynolds  numbers  on  the  blade  (based  on  distance  from  leading  edge). 

FOR  i  =  1  TO  2  'Use  Stagnation  on  Round-nose  Blunt  body 

rad  =  .005  'radius  of  blunt  nose,  m 

Re(i)  =  2  *  Vinf  •  rad  /  viso(i) 

NEXT! 

F0Ri  =  3T0  23 

Ro(i)  =  vel(i)  *  s(i)  /  visc(i) 

NEXT  i 

'  compute  local  resistivity  of  the  foil  in  ohm-m 
FOR  i  =  1  TO  23 

resistivity(i)  =  2.223  *  10  ''  -10  *  (Tscanfi)  •  1.8  +  32)  +  8.446  *  10  -7 

NEXT  i 

'  conduction  constants 

ukl  =  .026  'urethane  foam  thermal  conductivity,  (W/m*K) 
ak!  =  .0^  'air  thermal  conductivity  (W/m*K)  atSO  C 
'  convection  constants 

Cp!  -  1005.7  'for  air  at  30  C,  (J/kg*K) 

PrI  =  .707 

'  radiation  constants 

epsi  =  .17  'emissivity  for  polished  stainless  steel  at  300  K 
sigma!  =  5.67E-03  'Stefan-Bohzmann  constant.  (W/m''2*K^4) 

'  compute  the  Nusselt  and  Stanton  numbers  for  the  blade,  using 
’  steady-state  2-0  heat  model  (ignore  spanwise  and  chordwise  heat 
'  conduction  through  the  steel  foil  as  <  <  heat  convection). 
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ATot!  =  .013387  'area  of  foil  =  0.26352  m  x  0.0508  m) 


denomi  =  0 
denomr =  0 

FORi  =  1TO10  'Pressure  side 

Qcl  =  ukl*deltaT(i)/deital(i) 

Qrl  =  epsl  *  sigma!  *  (Tscan(i)  ^  4  -Tinfl  4) 
dem  =  (TscanO)  -  Tinfl) 

Qheatl  =  CurrI  ''  2  •  resistivity (0  •  7628  *  10  3  'W/m''2 

h{0  =  (Qheaf!-Qc!-Qr!)/dem 
Nu(i)  =  ABS(h(i))  •s(i)/akl 
St(i)  =  Nu(i)  /  (Pr!  *  Re(i)) 

’  also  calculate  the  theoretical  St#  and  Nu#  by  Kays  equs  9-50  and  12-32 


IF  i  =  1  OR  i  =  2THEN  GOTO  10 

denomi  =  denomi  -h  ((rho(i)  •  vel(i))  ^  1.87  *  (s(i)  -  s{i  - 1))) 

Sttheorl(i)  =  (.418  *  mu(i)  .5  •  (rho(i)  *  vel(0)  .435)  /  denomi  ^  .5 
NutheorlfO  =  Sttheorl(0  *  Pr  *  Ro(i) 

denomt  =  denomt  (((T8can(i)  -  Tinf)  1.25)  •  vel(i)  •  (s(i)  -  $(i  - 1))) 

Sttheort(i)  =  (.0287  *  Pr  -.4  *  ((Tscan(i)  -  Tint)  ^  .25)  *  visc(i)  ^  .2)  /  denomt  ^  .2 
Nutheort(i)  =  Sttheort(i)  •  Pr  *  Ro(i) 

K{0  =  (vi8c(i)  /  val(i)  ''2)*  ((vel(0  -  vol(i  - 1))  /  (s(0  -  s(i  - 1))) 

10  NEXT! 

'  the  theoretical  calculations  for  node  1 

denomi  =  (rho(1)  *  vol(1))  1.87  •  (s(1)  -s(2)) 

Sttheorl(l)  =  (.418  *  mu(1)  .5  *  (rho(1)  »  vel(1))  .435)  /  denomi  .5 

Nutheorl(l)  =  Stfheorl(i)  ^  Pr  *  Re  1) 

denomt  =  ((Tscan(l)  -  unf)  ^  1.25)  *  vel(1)  *  (s(1)  -  s(2)) 

Stfhoort(l)  =  C0287  •  Pr  -.4  •  ((rscan(l)  -  TinQ  .25)  •  visc(1)  ^  .2)  /  denomt  .2 
Nutheortfl)  =  Sttheort(l)  *  Pr*  Re(1) 

K(1)  =  (visc(1)  /  vel(1)  2)  *  ((vel(1)  -  vel(2))  /  (s(1)  -  s(2))) 

K(2)  =  0! 


FOR  i  =  1  TO  2  'Use  Stagnation  on  Round-nose  Blunt  body 
Nu(i)  =  .81  *  ((Re(i))  .5)  *  (Pr!  ''  .4) 

St(i)  =  Nu(i)/(Pr!*Re(i)) 

Nu(0  =  ABS(h(i))  *  8(0  /  ak! 

St(i)  =  Nu(0  /  (Pr!  *  Re(0) 

next) 

FORi  =  1T02 

Re(i)  =  vel(i)  *  s(i)  /visc(i) 

NEXTi 


’  Node  1 1  has  no  heat  generation  term  since  it  is  off  the  foil. 

Qc!  =  uk!  *  deltaT(ll)  /doltal(ll) 

Qr!  =  epsl  *  sigma!  •  ^8ean(11)  4 -Tinf!  ''  4) 

defT!  =  (Tscanfll)  -  Tinf!) 
h(11)  =  (-Qcl-Qr!)/derT! 

Nu(11)  =-ABS(h(11))*s(11)/ak! 

St(11)  =  Nu(11)/(Pr!*Re(11)) 

K(11)  =  (yi8c(11)  /vel(11)  ''  2)  *  ((vel(11)  .vel{10))/(s(11)  -  s(10))) 

'  the  calculations  for  Node  23 

Qc!  =  uk!  •  deltaT(23)  /deHal(23) 

Qr!  =  epsl  •  sigma!  *  (Tscan(23)  ''  4 -Tinf!  ''  4) 
defT!  =  (T8can(23)  -  Tinf!) 

Qheatl  =  Curr!  ''  2  *  rosi8tivity(23)  »  7628  *  10  ''  3  W/m ''  2 
h(23)  =  (Qheatl  -  Qc!  -  Qr!)  /  delT! 

Nu(23)  =  ABS(h(23))  •  8(23)  /aki 
St(23)  =  Nu(23)  /  (Pr!  •  Re(23)) 

denomi  =  denomi  +  ((rho(23)  •  vel(23))  1.87  •  (s(23)  -s(1))) 

Sttheorl(23)  =  (.418  •  mu(23)  •'  .5  *  (rho(23)  *  vel(23))  ''  .435)  /  denomi  .5 
Nutheorl(23)  =  Sttheori(23)  »  Pr  *  Re(23) 

denomt  =  denomt  +  (((Tscan(23)  -  Tinf)  "  1 .25j  •  vcl(23)  •  (s(23)  -  s(1))) 

Sttheort(23)  =  (.0287  *  Pr ''  -.4  *  ((rscan(23)  -  tinf)  .25)  *  vise(23)  ''  .2)  /  denomt  .2 
Nutheort(23)  =  Sttheort(23)  *  Pr  •  Re(23) 

K(23)  =  (visc(23)  /  vel(23)  ''  2)  •  ((vel(23)  -  vol(1))  /  (s(23)  -  s(1))) 

FOR  i  =  22  TO  12  STEP -1  'Suction  side 

Qc!  =  uk!  *  dettaT(i)  /deltal(i) 
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Qf!  =  eps!  •  »igma!  •  (Tseanfi)  ^  4-Tinf!  4) 
delT!  =  (T8can(i)  -  Tinf!) 

Qheat!  =  Curr!  ''2*  resistivftym  »  7628  *  10  3 

h(i)-=  (Qheat!  -  Qol  -  Qrl)  /  delTI 
Nu(i)  =ABS(h(i))  *#(!)/ aki 
St(0  =  Nu(0  /  (PrI  •  Re(i)) 


denomi  =  denomi  +  ((rho(i)  *  vel(i))  ^  1.87  *  (s{i)  -  s(i  +  1))) 

Stthoorl(i)  =  (.418  •  mu(i)  ^  .5  •  (rho(0  *  vel(0)  ''  .43^  /  denomi  .5 

Nutheorl(i)  =  Snheorl(i)  •  Pr  *  Re(0 

denomt  =  denomt  +  (((Tscan(0  -  Tinf)  1.25)  *  vel(i)  •  (8(i)  -  8(i  +  1))) 

Sttheort(i)  =  (.0287  *  Pr  -.4  *  ((Tscan(i)  -  Tinf)  ^  .25)  *  visc(i)  ^  .2)  /  denomt  ^  .2 
NiJtheort(i)  =  Sttheort(i)  *  Pr  •  Re(0 

K(i)  =  (visc(i)  /  vel(i)  ^2)*  ((vel(i)  -  vel(i  +  1))  /  (s(0  -  s(i  +  1))) 

NEXTi 

’  calculate  the  flat  plate  Nu#  for  comparison  purposes 
FOR  i  =  1  TO  23 

fttpltl(i)  =  .453  •  Re(0  ^  .5  *  .88917 
fHpltt(i)  =  .030784  •  Re(i)  ^  .8  *  .88917 
NEXTi 

'  save  data  on  disk 


OPEN  ■e:\sfeve\htxfef\temp\r6d'  +  file$  +  '.daf  FOR  OLTTPLrT  AS  #1 

PRINT  #1 ,  "Heat  transfer  results  for  run";  SPC(1):  fileS 

PRINT  #1,  USING  ‘Electrical  energy  input  =  ####.#  watts';  Qheat! 

PRINT  #1 ,  USING  ‘Uinf  =  ##.#  m/s‘:  Vinfl 
PRINT#1,“ 

PRINT  #1.  Tap  8/c  Re#  K  h  Nu#  NutI  Nutt  Nufpl  Nufpt* 

FOR  I  =  1  TO  23  H  H 

####.#  ####.#  ####.##  ####.#  ####.# 
R8(0:  K(i):  h(0:  Nu(i):  Nutheorl(i):  Nutheort(i):  fltpltl(i):  flfpltt(i) 

**■**"  ******  ####.#  ####.##  ####.#  ####.# 
****'*  '  '■  *(0  /  •1143;  Re(i):  K(0:  h(i):  Nu(0:  Nulheorl(i):  Nutheort(i):  mpltl(i) ;  fHpltt(i) 

CLOSE  #1 

'  write  the  nusslet  plot  files  ■nsttHX_.dar 

files  =  ‘NsltHX*  +  LTRIM$(STR$(alpha))  +  LTRIM$(STR$(test)) 

OPEN  dirS  +  files  + ‘.dat*  FOR  OUTPUT  AS  #1 
F0Ri  =  8T0  2STEP  -1 

PRINT  #1 .  USING  •##.####  ,  ##.## ^  ^  ' 

1143;  Reg^u(i);  Nutheorl(i):  Nutheort(0;  flfpltl(0:  Wpltt(0 


####.#,  ####.#,  ####.#,  ####.#  ,  ####.#•;  .s(i)  / 


print  the  suction  side  numbers  1 . 23  - 1 3 


FOR  i  =  1  TO  1 

PRIhTT  #1.  USING'##.####. 

;1143;  R8^.u(0:  Nutheorl(i):  Nutheort(i):  fftpHI(0;  fffpltt(i) 


####.#,  ####.#,  ####.#,  ####.# .  ####.#':  s(i)  / 


F0Ri  =  23T013STEP-1 
PRINT #1.  USING'##.####  ,  ##.##''  ^ 
.1143;  Reg^u(i):  Nutheorl(i):  Nutheort(i);  ftfptll(i);  fltpIttO) 

CLOSE  #1 

'  write  the  h  plot  files  ‘HHX_.dat' 


####.#,  ####.#,  ####.#,  ####.#  ,  ####.#';  s(i)  / 


files  =  'KH>C  +  LTRIMS(STRS(alpha))  +  LTRIMS(STR$(test)) 

OPEN  dirS  +  fileS  +  '.da#  FOR  OUTPUT  AS  #1 
F0Ri  =  8T0  2STEP.1 

PRINT #1,  USING '##.###,##.## - '; -8(1) /.1 143:  K(i) 

NEXTi 


print  the  suction  side  numbers  1, 23  - 13 
FOR  i  =  1  TO  1 

PRINT  #1.  USING'##.###  ,  ##.##''  ^  ^  s(i)  /  .1143;  K(0 

NEXTi  ” 


FOR  I  =  23  TO  13  STEP  -1 


PRINT  #1 ,  USING  •##.###  - - •;  8(i)  /  .1 1 43:  K(i) 

NEXT! 

CLOSE  #1 

INPUT  'press  ENTER  to  continue';  ansS 
NEXT  test 
NEXT  alpha 
END 


DEFINT  K-L 

SUB  PREDUCE  (file$.  Uinf,  CpQ.  pressureQ.  q) 

DIM  pvott8(28) 

OPEN  'e:\steve\eal\PCOEFF.dat*  FOR  INPUT  AS  #1 
INPUT  #1.m.b 
CLOSE  #1 


OPEN  'e:\3teve\ht)rfei\pras$ura\'  +  file$  +  'pv.daf  FOR  INPUT  AS  #1 
FORi  =  0TO28 

INPL/T#1.pvolts(i) 

NEXT! 

INPUT  #1,Ttun 
CLOSE  #1 


'  convert  the  voltage  readings  into  kPa  pressures 
FORi  =  0TO28 

^^res8ure(0  =  (m  *  pvotts(i))  +  b 

'  calculate  the  tunnel  q,  deltah,  and  velocity 
q  =  presaure(26)  •  pressure(28) 

density  =  (101.3  +  pre8sure(28))  /  (.28708  •  (Ttun  +  273.15)) 

deltah  =  q*  10001/2491 

Uinf  =  SQR(2000!  *  q  /  density) 

'  'collapse  the  vector*  by  eliminating  atmospheric  pressure  readings 
FOR  I  =  11  TO  20 


^^^_re88uro(0  =  pressure(i  +  1) 


FOR  I  =  21  TO  23 


^^^_res8ure(i)  =  pressure^  +  2) 


'  calculate  the  blade  surface  pressure  coefficients  (Cp) 
PRINT :  PRINT  'calculating  Cps' 

FOR  i  =  1  TO  23 


Cp(i)  =  (pressurelO  -  pressure(28))  /  q 
PRINT  USING  'pressure##  =  ##.###  kPa 
NEXT! 

PRINT  USING  ■Vinf  =  ##.#  m/s':  Uinf 
INPUT  'press  Enter  to  continue';  ans$ 

END  SUB 


CP##  -  ##.###';  I;  pressure(i);  i;  Cp(i) 
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VELREDUC.BAS 


•VELREDUC.BAS 

'This  program  reduces  the  binary-stored  fluctuation  data 

DECLARE  SUB  PLOT  (nplot!,  xmini,  xmaxi.ymin!,  ymax!,  cOlQ.  Cl  IQ.  C2!0) 

DECLARE  SUB  STATS  (N!.  El  IQ.  EZiQ.  MEANifl.  STDDEVfft) 

DECLARE  SUB  REDUCE  (N,  El  0.  E20.  ufile$,  vfile$,  caHileS) 

DECLARE  SUB  RealFT  (datQ.  N,  ISIGN%) 

DECLARE  SUB  Fourl  (datQ,  NN.  1SIGN%) 

DECLARE  SUB  SCALES  (EIQ.  E20.  N.  MEANIQ.  STDDEV'O.  MICROIQ,  IMTEGRAUQ,  df) 
DEFINTI-J 
'Main  Program  Block 
'SDYNAMIC 

'  outfileS  =  ’eAsteveLhtxfeAturbVhxturbb.dat* 

caKilaS  =  'e:\steve\ca(\044  84b.cof 
'  OPEN  outfiloS  FOR  APPEND  AS  #5 

N  =  8192  'samples 

fsample  =  20000!  'hertz 

dt  =  1  /(sample  'seconds 

df  =  fsample  /  N  'hertz 

nbins  =  1024 

DIM  El  (16384),  E2(16384),  frequ(16384) 

DIM  MEAN(2),  STDDEV(2),  MICRO(2).  INrEGRAL!(2) 

DIM  uspec(nbins),  vspec(nbins),  freq(nbins) 

DIM  cO(nbins),  Cl  (nbins),  C2(nbins) 

'  FORtest  =  1T0  6 

runs  =  'hxM' '  -t-  LTRIM$(STR$(test)) 

ufileS  =  ‘a;\turb\*  +  runS  +  'el  .daf  'these  files  are  the  binary  raw 
vfileS  =  'aAturbV  +  runS  +  •e2.daf  Xwire  voltage  files  from  VELACQ. 

'  read  in  the  x-wire  voltages  and  converted  to  velocities 

CALL  REDUCE(N,  El  0.  E20.  ufileS.  vfileS,  caHileS) 

CALL  STATS(N,  E10,  E20.  MEANQ.  STDDEVO) 

'  plot  the  first  1 00  u'  and  v' 


nplot  =100 

xmin  =  0;  xmax  =  nplot 
ymin  =  -.25  •  MEAN(1):  ymax  =  -ymin 
FOR  i  =  1  TO  nplot 
c0(i)  =  i 
C1(i)  =  E1(i) 


„  C2(.-)  =  E2(i) 

NEXT  I 


CALL  PLOT(nplot.  xmin,  xmax,  ymin.  ymax,  cOO.  Cl  0.  C20) 


'  obtain  the  FFT  and  energy  spectrum  and  stored  the  FFT  values  in  El  and  E2 


PRINT :  PRINT  *  performing  FFT  on  u1 '  and  ■; 

halfn  =  N/2 

CALLRealFT(E10.halfn,  1) 


'  calculate  the  right  side  of  Parseval's  theorem. 


sumcheck  =  2  •  E1(1)  2  +  2  •  El(2)  ''  2 

FORi  =  3TON-lSfEP2 

sumcheck  =  sumcheck  +  2  *  (E1(0  ''  2  +  E1(i  +  1)  ''2) 

NEXT! 

sumcheck  =  sumcheck  /  (N  +  1) 

PRINT  USING  “the  right  side  of  Parseval's  =  ##.### ^  sumcheck 


'  use  Parseval's  Theorem  to  obtain  E  (n)  from  H(u‘) 

'  u' 

freq(0)  =  0 

E1(0)  =  (2*(E1(1)  ''2/N)  ''2)/df 
FOR  i  =  3  TO  hatfn  +  1  'froqu(i)  arc  the  discrete  Froqs 

j  =  2  *  (i  -  2)  +  1  'El  (i)  are  the  spectral  energies 

frequ(i-2)  =  {i-2)  *  df 

E1(t-2)  =  2*  ((SQR(EIO)  ''  2  +  E1(j  +  1)  ''  2)/N)  2)/df 


NEXT! 

frequ(haHn)  =  10000! 

E1{hatfn)  =  (2?(E1(2)  ^  2/N)  ''  2)/cif 
’  PRINT'*  performing  FFT  on  u2" 

;  CALLReain'(E20.halfn,1) 

■  use  Parseval's  Theorem  to  obtain  E  (n)  from  H(v') 

’  v' 

’  FORi  =  3TOhalfn  +  1 

’  j  =  2  *  (i  -  2)  +  1  'E2(i)  are  the  spectral  energies 

;  ^'.2)  =  2 *  ((SQR(E20')  ^  2  +  E2G  -H)  ''  2)/N)  2)/cff 

;  E2(hatfn)  =  E2{haHn-1) 

’  calculate  the  integral  of  energy  over  the  frequency  for  comparison  to  u’ ''  2 

sumcheck  =  0 

FORi  =  3TOhalfn  +  1 

sumcheck  =  sumcheck  +  El(i-2) 

NEXT! 

sumcheck  =  sumcheck  *  (fsample/  N) 

PRINT 

PRINT  USING  *  The  energy  integration  =  ##.########';  sumcheck 

'  calculate  the  integral  scale  and  the  microscale 

CALLSCALES(E10,  E20,  N.  MEANQ.  STDDEVO,  MICROO.  INTEGRALIft,  df) 

TUINTENS  =  ((((STDD^(I)  ''2  +  STDDEV(2)  '' 2)  /  2) ''  .5)  /  MEAN(1))  *  100 

•  PRINT  #5.  USING  '& .  ###.##  .  ##.###  .  #.##### .  #.#####';  run$:  MEAN(1):  TUINTENS:  MICRO!(1); 
INTEGRAL!(1) 

‘  NEXT  test 

■  CLOSE  #5 

'  reduce  the  ener^  spectrum  file  from  n  points  to  256  points 
'  nperbin  =  N/nbins 

•  k  =  0 

'  FORi  =  1  TO  nbins/2 

'  uspec(0  =  0! 

'  vspec(i)  =  0! 

‘  freqfi)  =  0! 

'  FOR  j  =  1  TO  nperbin 

•  k=k+1 

’  U8pec(0  =  uspecW  +  E1(k) 

'  vspee{0  =  vspee(i)  +  E2(k) 

'  <feq(i)  =  heqji)  +  frequ(k) 

'  NEXTj 

'  uspec(i)  =  uspecfO  /  nperbin 
'  vspec(i)  =  V8pec(0  /  nperbin 
'  ^  nperbin 


plot  the  spectrum 

xmin  =  LOG(freq(1)):xmax  =  LOG{freq(nbins/2)) 

ymin  =  LOG(uspec(1)  / 100001):  ymax  =  LOG(10!  *  uspec(l)) 

FOR  i  =  1  TO  nbins/2 
cO(i)  =  LOG(freq(i)) 

C1(i)  =  LOG(uspec(i')) 

C2(i)  =  LOG(vspec(i)) 

NEXT! 

CALL  PLOT(nbins  /  2,  xmin,  xmax,  ymin,  ymax,  cOQ.  Cl  0,  C20) 

specfileS  =  'eitsteveWaftSPEC.DAP 
OPEN  specfileS  FOR  OUTPLTT  AS  #1 
FOR  i  =  1  TO  nbin8/2 

PRINT  #1,  freq(i),  uspec(i),  vspec(i) 

NEXTi 
CLOSE  #1 

OPEN  •e:\8teve\dat\logplotdat'  FOR  OUTPUT  AS  #1 
FOR  i  =  1  TO  256 

uspecm  =  (uspecfO  *  MEAN(1))  /  (STDDEV(I)  ''2*  INTEGRAL!(1)) 
freqm  =  (freq(0  *  INTEGRAL!(1))  /  MEAN(l) 

PRIm  #1,  uspec(i),  freq(i) 

NEXT! 

CLOSE  #1 
END 
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REM$STATIC 

DEFINTK-M 

SUB  Fourl  (datO.  NN.  ISIGN) 

N  =  2  *  NN 
i=1 

FOR  i  =  1  TO  N  STEP  2 
in >i  THEN 
■reMPR  =  dat(j) 

TEMPI  =  dat(j  +  1) 
dat(j)  =  dat(0 
datO  +  1)  =  dat(i  +  1) 
dat(i)  =  TEMPR 
dat  i  +  1)  =  TEMPI 
END  IF 
m  =  N/2 

1  IF  ((m  >  2)  AND  (j  >  m))  THEN 

I  =  j-m 
m  =  m/2 
G0T01 
END  IF 

11  N^i^"’ 

MMAX  =  2 

2  IFN>MMAXTHEN 

ISTEP  =  2*MMAX 

THETA  =  6.28318530717959#  /  (ISIGN  •  MMAX) 
WPR  =  -2#  •  SIN(.5#  •  THETA)  ''  2 
WPI  =  SIN(THETA) 

WR  =  1# 

Wl  =  0# 

FOR  m  =  1  TO  MMAX  STEP  2 
FOR  i  =  m  TO  N  STEP  ISTEP 
i=i  +  MMAX 

TEMPR  =  WR  •  dat®  -  Wl  *  dat(j  +  1) 

TEMPI  =  WR  •  dat(j  +  1)  +  Wl  •  dat® 
dat®  =  dat®  -  TEMPR 
datO +  1)  =  dat(i  +  1). TEMPI 
datiO  «  dat®  +  TEMPR 

+  1)  =  dat(i  +  1)  +  TEMPI 

WTEMP  =  WR 

WR  =  WR  *  WPR  -  Wl  •  WPI  +  WR 
Wl  =  Wl  •  WPR  +  WTEMP  •  WPI  +  Wl 
NEXTm 
MMAX  =  ISTEP 
GOTO  2 
END  IF 
END  SUB 


DEFSNG  l-M 

SUB  PLOT  (nplot,  xmin,  xmax,  ymin,  ymax,  cOQ.  Cl  0.  C20) 
CLS 

xmax  =  xmax  -  xmin 
ymax  =  ymax -ymin 
FOR  I  =  1  TO  nplot 

cO®  =  <5)®  -  xmin 
Cl  (i)  =  Cl  (i)  -  ymin 
C2®  =  02®  -  ymin 

NEXTi 
SCREEN  9 

UNE  (0. 0)T639. 319).  3,  BF 
FORi  =  0TO10 

Y  =  (1  -i/ 10)  *319 
UNE(0.Y)-(639,Y) 

NEXTi 

xO  =  (c0(1)  /  xmax)  *  639 
yO  =  (1-01(1) /ymax)  *319 
y02=  (1  -02(1) /ymax)  *319 
FOR  i  =  1  TO  nplot 

y1  =  (1  -  Cl® /ymax)  *319 
y2  =  (1-02® /ymax)  *31 9 
x1  =  (cO(i)  /  xmax)  *  639 
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NEXTi 
INPUT  ‘corrtinue';  ansj 
SCREEN  0 
END  SUB 


DEFINTI-M 

SUB  RealFT  (datQ.  N.  ISIGN) 

THETA  =  6.28318530717959#  /  2#  /  N 
Cl  =  .5 

IF  ISIGN  =  1  THEN 
C2  =  -.5 

CALLFour1(datO,  N,  1) 

ELSE 
C2  =  .5 


THETA  =  -THETA 
END  IF 

WPR  =  -2#  *  SIN(.5#  *  THETA)  ^  2 
WPI  =  SIN(THETA) 

WR=  1#  +  WPR 
Wl  =  WPI 
N2P3  =  2  •  N  +  3 
FORi  =  2TON/2  +  1 

11  =  2  *  i  - 1 

12  =  II  +  1 

13  =  N2P3-I2 
M  =  13  -H 
WRS  =  WR 
WIS  =  Wl 

H1R  =  Cl  •  {dat(ll)  +  dat(l3)) 

H1I  =  C1  *(dat{12)-dat(l4)) 

H2R  =  -C2  *  (dat(l2)  +  dat(l4)) 

H2l  =  C2*{dat(l1)-dat(l3)) 
dat(ll)  =  HI  R  +  WRS  •  H2R  -  WIS  •  H2I 
dat  la  =  HI  I  +  WRS  *  H2I  +  WIS  *  H2R 
dat{l3)  =  HI  R  •  WRS  •  H2R  +  WIS  *  H2I 
dat[M)  =  -HI  I  +  WRS  *  H2I  +  WIS  •  H2R 
WTEMP  =  WR 

WR  =  WR  •  WPR  -  Wl  *  WPI  +  WR 
Wl  =  Wl  •  WPR  +  WTEMP  *  WPI  +  Wl 
NEXTi 

IF  ISIGN  =  1  THEN 
H1R  =  datm 
dat(1)  =  HIR  +  dat(2) 
dat(2)  =  H1R-dat(2) 

ELSE 

HlR  =  dat(1) 

dat(1)  =  Cl  *  (HIR  +  dal{2)) 
dat(2)  =  C1»(H1R-dat(2)) 

CALL  FourUdatO,  N,  -1) 

END  IF 
END  SUB 


DEFSNG  l-M 

SUB  REDUCE  (N.  El  Q.  E20.  ufile$.  vfile$.  calfileS) 
'  subroL'tine  to  reduce  hot  wire  data 
CLS 

REDIM  B0(2),  B1{2) 

'  read  in  calibration  constants 

OPEN  caHileS  FOR  INPUT  AS  #1 
FOR  i  =  1  TO  2 
INPUT#1.BO(i),B1(0 
NEXTi 


CLOSE  #1 

’  read  in  hotwire  r:tta 

PRINT  I  JG  •  Reading  in  Hot  Wire  Data  from  &*:  ufileS 
i  =  0 

OPEN  ufile$  FOR  BINARY  AS  #1 
OPEN  vfileS  FOR  BINARY  AS  #2 
FORj  =  lTON 
i  =  i  +  1 
GET#1..E10 
GET  #2, .  E20’) 

NEXT! 

CLOSE  #1 
CLOSE  *Z 
NPOINTS  =  i 

I 

'  calculate  wire  effective  velocities 


PRINT:  PRINT* 
FOR i  =  1  TO  N 


E1(n  =  f(E1j5  - 


Calculating  Wire  Effective  Velocities* 

2-B0(1))/B1(1))  ''  2.22 
2-B0(2))/B1(2)j  ""  2.22 


'  calculate  velocities  in  wire  coordinates,  then  lab  coordinates 


PRINT :  PRINT  *  Calculating  Velocities* 
k  =  .15 

pi  =  4!»ATN(1!) 
mu  =  45! 


mu  =  (mu  / 1801)  •  pi 

FOR  i  =  1  TO  NPOINTS 

XX=  (E1(0  2-(k*E2R)  **  2)/(l!-k 

IFXX<0!THENXX  =  X  ''  2 

X  =  SQR(XX) 

YY=(E2(0  ^2-{k»E1(0)  ''2)/(1!-k 
IFYY<o!THENYY  =  Y  2 

Y  =  SQR(YY) 

U  =  COS(mu)  •  X  +  SIN(mu)  •  Y 

V  -Slf^mu)  *  X  +  COS(mu)  *  Y 

irr 

'  U0  is  now  stored  in  El  0,  V0  in  E20 
END  SUB 


4) 


4) 


DEFINTW 

SUB  SCALES  (El  Q.  E20.  N.  MEANIQ,  STDDEViO,  MICROIQ.  INTEGRAUQ.  df) 
sum  =  0 

FORi  =  3TO(N/2)  +  1 

sum  =  sum  +  ((i-2)  ''  2)  •El(i-2) 

NEXT! 

sum  =  sum  *  df  3 

MICRO(I)  =  SQR(1  /  ((sum  •  2 *  3.14  2) / (MEAN!(1)  ''2*  STD0EV!(1)  2))) 

sum  =  0 
F0Ri  =  3T012 
sum  =  sum  +  El  (i  -  21 
NEXTi 

sum  =  sum/ 10 

1NTEGRAL!(1)  =  (MEANI(I)  /  (4!  *  ST0DEV!(1)  2))  *  sum 
PRINT 

PRINT  USING  *  the  microscale  =  #.#########  m*;  MICRO!(1) 

PRINT  USING  *  the  integral  scale  =  #.#########  m*;  INTEGRAL'(I) 

■  INPUT  'press  ENTER  to  continue*;  ansS 
END  SUB 
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DEFSNQ  kJ 

SUB  STATS  (N.  EIQ,  E20.  MEANQ,  STDDEVQ) 

'  calculatelha  means  of  u'  and  v* 

PRINT :  PRINT  ■  Calculating  Mean  and  Standard  Deviation* 

MEAN(1)  =  01 
MEAN(2)  =  0! 

FOR i  =  1  TO  N 
MEAN(1)  =  MEANO)  +  El  (I) 

MEAN(2)  =  MEAN(2)  +  E2(i) 

NEXT  i 

MEAN(1)  =  MEAN(1)  /  N 
MEAN(2)  =  MEAN(2)  /  N 

’  calculate  the  standard  deviations  of  u'  and  v'  and  remove  the  mean  for  FFT 
STDDEV(I)  =  01 
STDDEV  2i  =  0! 

FOR  i  1  TO  N 

El  (fl  -  El  |0  -  MEANflf  'removing  mean  for  u’ 

E2(i)  E2(i)  •  MEAN(2)  'removing  mean  for  v' 

STDDEV(I)  =  STDD^O)  +  El(i)  ''  2 
STDDEV{2)  =  STDDEV(2)  +  E2(0  ''  2 
NEXTi 

STDDEVO)  =  SQR(STDDEV(1)  /  N) 

ST0DEV(2)  =  SQR(STDDEV(2)  /  N) 

PRINT:  PRINT*  mean  std dev* 

FOR i  =  1  TO  2 

PRINT  USING  *  +##.###  +##.###*:  MEAN(i):  STDDEV(i) 

NEXTi 

'  calculate  u'  ^  2  mean  to  compare  to  the  energy  integral, 
sumcheck  0 
FOR i  =  1  TO  N 

sumcheck  =  sumcheck  4  El  (i)  ''2 
NEXTi 

sumcheck  =  sumcheck /N 

PRINT :  PRINT  USING  *  u"'2  mean  =  ###.########  and  *;  sumcheck; 

'  calculate  the  left  side  of  Parseval's  theorem 

as  ^  Kl 

PRINT  USING  “the  left  side  of  Parseval's  =  ##.###^'''"'*:  sumcheck 
END  SUB 


F-34 


Appendix  G:  Calibration  Procedures 


Detailed  pressure  transducer  and  x-wire  calibration  procedures  are  listed  below, 
complete  with  example  calibration  results. 

Pressure  Transducer  Calibration.  The  Scanivalve  calibration  setup  is  illustrated  in 
Figure  G1 .  Pressure  is  applied  to  Scanivalve  using  shop  air  and  vacuum  is  applied 
using  a  vacuum  hand  pump.  The  "known"  pressure  is  measured  with  a  King  Instru¬ 
ments  25"  water  manometer  with  1 .5  digits  of  accuracy. 

Once  executed,  the  PRESSCALbas  pressure  transducer  calibration  program 
prompts  the  user  to  trigger  a  voltage  reading  and  input  the  associated  pressure  in 
inches  of  water  (inc/ies// 2  0)-  The  program  converts  inchesHzO  to  kPa,Vne 
units  required  for  data  reduction.  After  all  the  calibration  data  is  collected,  it  must  be  to 
a  straight  line  using  any  curve  fitting  technique.  The  calibration  steps  are  listed  below: 

1 .  connect  the  Scanivalve  as  shown  in  figure  G1 ; 

2.  execute  the  PRESSCALbas  program  and  follow  the  program  prompts; 

3.  starting  with  atmospheric  pressure,  record  the  pressure  transducer 
pressure  and  voltage  using  the  PRESSCAL  program; 

4.  increment  the  pressure  and  record  the  voltages; 

5.  once  the  upper  limit  of  the  pressure  transducer  range  is  reached, 
increment  back  down  to  atmospheric,  recording  the  voltages; 

6.  repeat  steps  3  -  5  for  vacuum; 

7.  after  recording  the  pressures  and  voltages,  use  any  straight  line  curve 
fitting  technique  to  determine  the  intercept  and  slope  of  the  linear 
pressure  transducer  voltage  /  pressure  relationship  and  write  the 
constants  to  a  files. 
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Using  a  +/-  2  psi  pressure  transducer,  the  calibration  curve  generated  for  this 
experiment  is  shown  in  Figure  G2.  The  cun/e  demonstrates  the  linearity  of  the  pressure 
transducer  over  the  operational  range. 

X-wire  Calibration.  The  X-wires  for  this  experiment  are  calibrated  with  known 
velocities  generated  with  a  standard  calibration  tank,  illustrated  in  Figure  G3.  Shop  air 
enters  the  tank,  passes  through  baffles  into  the  stilling  chamber,  then  passes  through  a 
low-loss  nozzle,  exiting  to  atmosphere  in  a  uniform  jet.  The  low-loss/isentropic  nozzle 
insures  that  total  pressure  will  be  conserved  and  the  jet  velocity  V ,  can  be  calculated 
based  on  the  tank  pressure: 


where  P  y  is  the  tank  pressure 

P ,  is  the  pressure  of  the  jet 
p !  is  the  density  of  the  jet 

For  this  calibration  configuration,  P ,  is  atmospheric  pressure  and  p ,  is  calculated 

from  the  ideal  gas  law  based  on  P ,  and  the  tank  air  temperature.  The  tank  air  temper¬ 
ature  is  measured  by  the  data  acquisition  system  and  the  tank  pressure  is  measured  by 
the  user  with  a  20"  Miriam  Micromanometer  with  4  digits  of  accuracy.  The  user  enters 
the  tank  pressure  in  Laches  HzO  and  the  calibration  program  calculates  V 

An  elaborate  calibration  scheme  is  described  in  the  works  of  Galassi  and  Acree  in 
which  the  X-wire  probe  is  1)  examined  microscopically  to  determine  the  bisector  angle; 
2)  eyeball  aligned  in  the  calibration  stand  probe  fixture;  and,  3)  calibrated  through  both 
a  velocity  range  and  an  angle  range  (by  rotating  the  probe).  After  the  calibration  data  is 
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collected,  the  data  is  manipulated  to  determine  the  calibration  constants  a  and  b  and 
the  correction  factor  k  ,  and  then  the  probe  is  removed  from  the  calibration  stand  and 
eyeball  aligned  in  the  tunnel  probe  fixture. 

For  this  experiment,  a  new  tunnel  probe  fixture  was  fabricated  which  also  serves 
as  the  calibration  stand  probe  fixture,  eliminating  a  realignment  of  the  probe; 
additionally,  a  new  probe  alignment  procedure  was  developed.  Instead  of  measuring 
the  bisector  angle,  software  was  incorporated  which  uses  the  data  acquisition  system  to 
measure  the  wire1  and  wire2  voltages  for  a  fixed  jet  velocity.  For  a  bisector  angle  of  90 
degrees  and  identical  wires,  the  voltages  will  be  identical  and  the  bisector  will  be 
aligned  with  the  jet  when  the  sum  of  the  voltages  is  maximized.  Although  no  X-wires  are 
"perfect",  this  method  was  determined  to  be  as  good  as  eyeballing  the  alignment. 

Another  change  to  the  X-wire  calibration  procedure  was  the  elimination  of  the 
angle  calibration.  As  mentioned,  studies  (Champagne,  1 965)  have  found  that,  for  a 
limited  velocity  and  angle  range,  k  varies  linearly  with  21/d  of  the  probe  from  k  =  0.2 
for  21/d  =  200  to  A:  =  0.0  for  21/d  =  600  to  800.  With  2l/d  =  300,  our  TSi  T1 .5  probes 
have  k  =0.15. 

The  new  X-wire  calibration  procedure  is: 

1 .  configure  the  calibration  stand  as  illustrated  in  Figure  G4; 

2.  install  the  probe  in  the  calibration  stand,  paying  attention  to  the  flow 
geometry  and  wire  positions; 

3.  execute  XWIRECALbas  and  follow  the  instructions  and  prompts; 

4.  align  the  probe  using  the  "Probe  Alignment"  menu  selection  and 
maximizing  the  voltage  for  a  fixed  jet  velocity; 

5.  calibrate  the  probe  by  selecting  "Calibrate  the  Probe"  and  following  the 
instructions; 

6.  calculate  the  probe  calibration  coefficients  by  selecting  "Calculate  Coeffi¬ 
cients"  and  following  the  instructions; 
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7.  remove  the  probe/probe  fixture  and  install  in  TCTF. 

A  typical  calibration  curve  is  illustrated  in  Figure  G5.  Since  no  temperature 
compensation  is  applied  in  this  calibration  configuration,  a  new  calibration  must  be 
performed  for  a  change  in  temperature  of  more  than  5  degrees  F. 


Sconivolve 


monometer 


Figure  G1  Pressure  Calibration  Setup 


Figure  G2  Typical  Pressure  Calibration  Curve 
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Figure  G3  Hot  Wire  Calibration  Tank 
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Figure  G4  Hot  Wire  Calibration  Setup 
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